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4 THÉORÈMES GÈiSÉRÀUX SUR LE POTESTIEL NEWTONIEN 

Théorème de Qreen. — Si U et V sont deux fonctions con- 
tinues à rintérieur d'un volume T limité par une surface S, on a 
Tégalité 

(0 / U-^t/cr= / UAWt 



1 \]~-d^= l UAVi 

t/ 8 t/ * 



f \ àx àx ày dy dz dz J 



Dans cette égalité, on peut permuter U et V, et en retranchant 
les deux égalités obtenues, on a 



(^) / (u^-^-S-)'^' = / (U^V-VAU)rf.. 





Cas particuliers de Végalité (i) : faisons U ^ i. 
(3) I ^rfcr= 1 AVrfT. 




an 




Faisons : 



U = V 




(4) / V^cfT:= / VAVrfT 



ÔV 








)-(-^r-(4^)']- 



Si V est une fonction satisfaisant à Téquation de Laplace, 
AV=o, 



(^) / v^-=/[(-^)V(-^)V(4^)>, 



le second membre étant positif, le premier Test aussi. 

Si U et V sont des potentiels, sur une sphère de rayon suffi- 
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samment grand, U et V sont de l'ordre de — , leurs dérivées 

premières de Tordre de -^ . On peut alors appliquer la for- 
mule (i) à l'espace extérieur a une surface S et intérieur a une 
sphère de rayon très grand. Si ce rayon est suffisamment grand, 
l'intégrale prise sur la surface de la sphère © est négligeable : 

c'est-h-dirc est de Tordre de -5- ; mais on a la relation 






car les éléments de l'intégrale se détruisent deux a deux. 

Quant aux intégrales des seconds membres, leur somme donne 
des intégrales étendues à l'espace entier, et on a, en ajoutant les 
deux égalités obtenues : 

m C^^K^TJ , C I ^'^ ^^ , ÔV dU , dV dU\ , 

(6) j \51Yd. + j ^__+__+— _jrf.=o. 

les intégrales étant étendues à tout l'espace. 
De même, la formule (2) donnera 



(UAV— VAU)rfT = o; 

et, en se rappelant que 

— AV = 4*^^? 

— AV,= 4Tcp^, 

on a 

(7) / (V,p-VMrfT=o 

l'intégrale étant étendue à tout l'espace. 
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Cette formule est analogue à celle que Ton trouve dans la 
théorie de réleclrîcitc, 



2j'^^^'~Aj"^'^^^' 



Dans la formule (7), faisons 

on en conclut la formule 

(8) / (Vrf? — crfV)rfT = o. 



Travail de F attraction. — Supposons un système de masses 
m'^m"^ ..., etc. attirées par un autre /w'^, m'\^ ...., etc. 

Soit V'p V'j,... le potentiel aux points ni' ^ m"^,,,, etc., dû aux 
masses attirantes m'j, ?n'\^ ..., etc. 

V, V..., le potentiel aux points m\^ m'\,,,^ dû aux points 



m\ m" y 



Les masses attirées se déplaçant, on a comme travail 
en appelant ox, 01/, o.Zy le déplacement du point de masse m ; 



Si je pose 



j aurai 






t = 511. 



Si les masses attirantes se déplacent aussi, mais seules, le 
potentiel \^ deviendra 



# « 
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et on a 



ni=\mrjy[. 



Si les deux déplacements ont lieu simultanément, on a 
Supposons que les masses attirées forment un volume T, on a 




11= I pV^rf-, 



cette intégrale pouvant d'ailleurs s'étendre à l'espace entier, 
puisque p est nul en dehors de T. On a 

C . 

oll = I V^SprfT 

étendue de même h l'espace entier. 

Si la masse attirante est la même que la masse attirée, ou 
a V = V., 

0ll= / \XorfT, 

ce qui peut s'écrire en tenant compte de (8) 

l\ 

oH = I poVc/T, 



ou encore 



311=. /i!X±I!Lrf,. 



Posons maintenant 



w=J^.,. 
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rîntégrale étant toujours étendue a tout l'espace ; W est Ténergie 
du système. 
On a 

On peut d'ailleurs écrire, puisque AV = — 4^^?» 



W = -g^ I VAVrft 




et, en tenant compte de (6), 



(9) 



rintégrale étant toujours étendue à l'espace entier. 

Conditions d'équilibre. — Ceci étant rappelé, nous allons com- 
mencer Tétude du problème que nous nous sommes proposé. 

Considérons une masse fluide, isolée de toute influence exté- 
rieure, tournant autour d'un axe fixe que nous prendrons comme 
axe Oz^ d'un mouvement uniforme de vitesse angulaire (o; prenons 
les axes Ox et Oy, invariablement liés h la masse fluide. 

Soit/; la pression en un point (or, ?/, z); p ne dépend que de 
.r, y, Zy la force s'exerçant sur l'élément de volume dt a pour 
composantes 

ox Oy oz 

Ecrivons les conditions d'équilibre relatif, en remarquant que 
l'accélération de Coriolis est nulle. Si on appelle X, Y, Z les 
composantes suivant les axes de coordonnées de la force agis- 
sant sur la molécule (.r, ?/, z)y on a 

ôV 
V ^V , , 
„ OV 



et, en posant 
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U=V+~(a:^+r), 



on a 



^n dU 

— 



ôa; ' (Sx ' 
dp dU 






d'où Ton tire 



dx-\ — r^rfy +-^rfc==:p (-r^ — dx -{- - — dij -\- -—- dz] , 

dp = pd\J, 



Cette relation montre que p est uniquement fonction de p ; U est 
par suite, aussi uniquement fonction de p et, par suite, de s. 
On arrive alors au théorème suivant: 

Théorème. — Quand une masse fluide, animée d'un mouve- 
ment de rotation autour d'un axe fixe^ est en équilibre relatifs les 
surfaces de niveau sont les surfaces d'égale pression^ et aussi 
d'* égale densités 

A la surface de la masse, on a /> = o. Donc U est une cons- 
tante. La surface libre est donc une surface de niveau. 

S'il n'y a pas de rotation, on a U = V. 

Dans le cas général on a AV:=: — 4^p. AV est donc simple- 
ment fonction de U. 

p 1 r ♦ . 1 ^U ÔU OU ^ , , . 

nemarquons que la force totale — r — ►, -^r — > -^r— est en gênerai 

normale aux surfaces de niveau, U = C. 

Ces conditions d*équilibre sont nécessaires. 

Une condition nécessaire et suffisante d'équilibre est que le 
travail résultant d'un déplacement virtuel soit nul. Ce travail 
comprend le travail de l'attraction, plus le travail à(\ h la force 
centrifuge. 



• « 



'.à 



|t 
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Le premier a pour expression 8W, que nous avons déjà 
étudié, le second a pour expression 



pd'z {{i>^xox + w^yôy) 




w 




û— -. I p (.r2 + y2)tfT; 



et, si J est le moment d'inertie par rapport a Oz, le travail dû à 
la force centrifuge est 



J. 



La condition d'équilibre est donc 



5W+— oJ=o 



pour tout déplacement compatible avec les liaisons. 
Cette condition n'implique pas nécessairement que 






W + — J 

2 



soit maximum ; il faut encore d*autres conditions ; mais je dis 
que si cette fonction est maximum, l'équilibre est stable. 
:[ Soit T la force vive relative, on a l'égalité 



_/w+-^j)=c". 



to*" 



Posons W + — J =U, et supposons que pour le maximum 

U soit égale à U^. 
On a 



fc T-U=T,-U,. 



RELATION ENTRE LA MASSE, LE VOLUME ET LA VITESSE n 

Si on donne maintenant à U une valeur U<Uq et à T,, une 
valeur a, on aura 

T-U = T,-U„ 

T = a-(U.-U). 

La force vive dans les instants suivants sera inférieure à Ty, 
l'équilibre est donc stable. 

Relation entre la niasse, le volume et la vitesse de rotation. 
— Une condition nécessaire de l'équilibre est, qu'en tout point 
de la surface libre, la force soit dirigée vers l'intérieur, autrement 

■VT ▼ 

une partie se détacherait. Il faut donc -^ — <o et par suite 




La formule (3) transforme cette inégalité en la suivante : 




AUrfT<o, 



Donc il faut que Ton ait 








c'est-a-dire 




2tù~(h<i0, 



En désignant le volume par T et la masse par M, il vient 

— 4^M + ato'^T < o. 

Dans le cas particulier oii la densité est constante et égale à p, 
on a 

20)" <47^p. 
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Dans rapplîcatîon de cette formule, il faut se rappeler que les 
unités sont telles que la force d'attraction newtonienne s'ex- 
prime par la formule 

^ mm' 

r=-—-r- 
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LA MASSE EST SANS MOUVEMENT DE ROTATION 

Simpliûcation des formules générales. — Etudions le cas 
particulier d'une masse homogène fluide et sans mouvement de 
rotation. 

En appliquant les théorèmes généraux, on trouve que la sur- 
face libre doit être une surface équipotentielle V = V,,. 

Si on prend la densité de la masse pour unité de densité, le 
potentiel en un point a pour expression 



V = 




J/énergie est 




pVrfT / VrfT 




V est continu dans tout l'espace ainsi qu'à la surface libre ; il est 
nul h Tinfini ; a l'extérieur il n'a ni maximum nî minimum et est 
compris entre o et V^ ; a l'intérieur il est supérieur a V„ et a un 
maximum. 



L'intégrale / r — di étendue h une surface fermée cntou- 
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raiit le corps a pour valeur / AV^/t, étendue au volume limité 

A\Wt = — 4^ / ?ch=^ — 4'ï^^I-» 

Ce théorème est encore vrai pour la surface limite, et il s'ap- 
plique ([ue le corps soit creux ou plein, et môme si Ton a affaire 
à une distribution superficielle. 

Comme enfin nous supposons le corps homogène et de den- 
sité I, on a 





d7 ==: I AV<fT= — I 4'^p^'z= — 4'^T 



en désignant par T son volume, et par S une surface qui le 
renferme complètement, en particulier sa propre surface. 

Si on porte sur la normale en M h la surface extérieure un 
segment MM' = rf/îp(fig. 2), le potentiel en M'a pour expression : 

Or, on a 




on en conclut que 



OV 

<o. 



dn 



t 



On a vu que cette condition était nécessaire pour l'équilibre : 
ici elle est toujours remplie. 

Soit Vq le potentiel de la masse fluide sur sa surface, et suppo- 
sons maintenant une couche électrique de masse M, en équilibre 
à la surface du corps, et produisant dans l'espace un potentiel 
V'tel que sur la surface et à son intérieur le potentiel soit y=y^, 

A l'extérieur, on a \' =V. En ellet, sur la surface V' = Vo et 
à l'infini Y'= o. Or V et V sont deux fonctions satisfaisant à 
l'équation de Laplace, ayant les m(>mcs valeurs sur la surface et 
sur une sphère de rayon très grand. Ces deux fonctions coïnci- 
dent donc : c'est le principe de Dlrichlct. 
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On a, d'après des formules connues, 



Ù/l 



On a d'ailleurs 




^"^ ' ^'^''=-^'^' 



Donc la masse M répandue à la surface a pour valeur T et réci- 
proquement si 

M = T, V=:V. 

T , 

Le rapport —:r7— est ce que Ton appelle la capacité électrique 

' 

du corps. 

Cette notion va nous ôtre utile pour résoudre le problème 
proposé. 

Théorie de Liapounoff, — Parmi les figures d'équilibre, il y 
a une solution évidente qui est la sphère : mais peut-il y avoir 
d'autres figures d'équilibre ? Nous n'en savons rien : mais si Ton 
se propose de chercher les figures d'équilibre stable, le pro- 
blème est différent. 

Le maximum absolu de W est bien atteint pour la sphère, 
mais a priori nous ne pouvons affirmer que W n'ait pas d'autres 
maximums relatifs. Cela a été démontré par Liapounoff; la 
démonstration est assez longue et on va procéder par étapes. 

L'énergie a un maximum. — Le potentiel d'un volume 
quelconque par rapport à un point quelconque, est inférieur au 
potentiel d'une sphère de même volume par rapport à son centre. 

Du point M, comme centre, je décris la sphère i] ayant le même 
volume que le corps S (fig. 3). 

Le corps S est partagé en deux volumes : le volume K inté- 
rieur à S, et le volume A extérieur. 

Soit B le volume de S extérieur à S. 



i6 
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On a par hypothèse B = A. 

Si U est le rayon do hi splière, on a 




Potentiel de A < 



Potentiel de B 



11 
B 



> 



U 



A^ 
R 



Fig. 3. 



D'où 

Potentiel de B > Potentiel 
de A 

Potentiel de S = Potentiel 
de K-j- Potentiel de A 

Potentiel de il = Potentiel 
de K -j- Potentiel de B. 



Donc 



Potentiel de 1' > Potentiel de S. 
Si V est le potentiel au point M, et R le rayon de la sphère S, 
on a donc 



V<2-R-. 



L'énergie qui est égale h 



\(h 



2 



est donc inférieure à 



-izW I (k 



et il vient 



W<7:1VT 



en appelant T le volume de la niasse. 
On a d'ailleurs 

T=-|-7;R\ 



d'après la définition de R. . 
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Donc 

o 

L'énergie W due h un volume T est donc limitée supérieure- 
ment ; W a donc un maximum. Je dis que ce maximum n*est 
atteint que si la surface extérieure est une sphère. 

Nous allons d'abord démontrer le théorème suivant : 

La capacité électrostatique a un minimum. — De tous les 
corps ayant même volume, la sphère est celui qui a la moindre 
capacité électrostatique. * 

En effet, la formule (9) (p. 8) donne 

-=^/[(-|^)"+(f)V(-^)*]-. 

l'intégrale étant étendue a tout l'espace. 

Supposons que sur un corps conducteur de volume T il y ait 
une masse électrique T en équilibre. Si V est le potentiel corres- 
pondant, sur la surface du corps on a 

CV' = T 

par définition de la capacité C du conducteur. 

A l'intérieur du corps, V est constant ; à l'extérieur, V satisfait 
à l'équation AV' = o et est nul à l'infini. 

L'énergie électrique est, d'après une formule connue de l'élec- 
tricité : 

1 I T^ 

W' = — CV'-= — 4t-- 

2 2 C 

Elle est d'ailleurs égale à 



l'intégrale étant étendue, soit h tout l'espace, soit au volume 

PoiNCAR^. — Figures d'équilibre. a 




r8 
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extérieur seulement, puisque V est constant h Tintérieur du 
volume. 

Nous allons voir qu'elle est inférieure à Ténergie W duc au 
potentiel newtonien du corps supposé homogène et de den- 
sité égale à V unité. 

Soit V = V + U ce potentiel : je pose 



W = 




)[^^^]'+[^^]-[^^^]'i-. 



l'intégrale étant étendue à tout respacc. 



W 




d\ 



-/ \ i 



d.v 






d\' ^ '- 



dz 



V- 




ùv OU ov ou ov ou 



dx dx dij Oy 0; 0; 



rfT 




ou 

0.r 



-'f)'+ 



mh- 



La première intégrale est W. 

La deuxième se décompose en deux : Tune prise h l'intérieur 

dV OV OV 
du corps, qui est nulle, puisque -:- — == - — = — — =: o. 

* ^ i^X 01J 03 

L'autre, prise à l'extérieur, se transforme par la formule de 
Green en 




,, ou 



0« 



d'J — 



V'AUiT. 



c 



La seconde de ces deux intégrales est nulle puisque AU est nul 
a l'extérieur. 

La première est également nulle : en effet, sur la surface, V 
est constant et égal à Vy. 
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L'intégrale est donc 



ï9 




or 



d7 = ^7zT, 



0/2,. 






dn. 



cette intégrale est donc nulle. 
On a par suite 






La troisième n'est pas nulle, car U n'est pas nul à l'intérieur 
du corps, puisque V y est constant, et que V y varie. 
La troisième intégrale est alors positive, et on a bien 

On a donc, en réunissant les deux résultats précédents : 

7rR^T>W> ^' 



2C ' 

d'oii 

T 
C> 



27îR^ ' 

et enfin 



!àO 
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C a donc un minimum. Je dis que ce minimum ne peut être 
atteint que pour la sphère. 

Le minimum est atteint pour la sphère. — Supposons qu'un 
corps conducteur ait la forme pour laquelle ce minimum est 
atteint : on a 



^ c 4 



kfmH^ïHw- 




-rr — désignant la dérivée suivant la normale à la surface équî- 
potentielle passant par le point (.r, ?/, ^), 



2W = 



Déformons le conducteur de façon que l'élément de surface 

rfT devienne d^ (fig. 4). Appelons 
^ la projection sur la normale à la 
surface du déplacement dv compté 
avec son signe. 

L'augmentation du volume est 
alors 




Fit?. 4. 



L'augmentation de l'intégrale est 



j;./t. 



INl 



s^=/t[(^)"-(f)"+(^)'] 




Ceci suppose que le potentiel reste constant; mais, par hypo- 
thèse, la figure est en équilibre. Donc il y a un minimum pour 
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-p- ; la variation de Tintégrale due à la variation du potentiel 

est donc nulle, puisque cette intégrale est minimum. 
Or, si on appelle [x la densité superficielle, on a 

ÔV . 

= — 47r[A. 



ô/i 



On a donc 



-1-— — = I 47ra^ÇrfT 



Or, il faut que dC =o : donc 

lK%d<T = o 

et comme on doit avoir en même temps / ÇrfT:=:o, quel que 

soit Ç, il faut que Ton ait [jl = C'^ 

On a d'ailleurs M = u.S en appelant S la surface du conduc- 
teur : on a donc 

Supposons maintenant que le conducteur se déforme en restant 
semblable à lui-même : la capacité est proportionnelle à la 
racine cubique du volume, et on a 

rfC I dT 



C ~ i T ' 
d'où on conclut 

S^= i27wTC. 
Tous les corps dont la capacité est minimum sont donc tels 
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que cette relation existe entre leur capacité, leur volume et leur 
surface. 

Mais alors, pour tous les corps de même volume, la capacité est 
minimum si la surface est minimum ; or, parmi tous les corps 
de même volume, la sphère est celui qui a la surface minimum. 
C'est donc celui dont la capacité électrostatique est mini- 
mum. 

Nous arrivons maintenant au théorème que nous nous étions 
proposé de démontrer. 

La sphère est la seule ûgure d'équilibre. — Je dis que W est 
maximum quand le corps de volume T a la forme d'une sphère. 

Soit T le volume correspondant, donnons-lui un accroisse- 
ment (fig. 4) 



(n = 




Comme on a 



W= 





2 



V dp (k. 



Le corps étant supposé de densité égale à i , dp est égal h — i , 
G ou I ; dp n'est différent de zéro que dans le voisinage de la 
surface, qui est une surface équipotentîelle, le corps étant sup- 
posé en équilibre ; c'est-a-dire que V est constant. 



d\\=\\ I dpd^ 
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et rintégralc est égale à rfT. On a donc 







On peut établir une autre relation entre rfW et dT, 
Supposons que le corps varie en restant semblable à lui-même : 

Ténergie varie en restant proportionnelle à la puissance — du 

o 

volume. On a donc 

W ~ 3 T ■ 

Eliminant nîW et dT entre ces deux dernières relations, on 
trouve 

W = 4-VJ. 
;> 

D'après la définition de W, on voit que Ton peut appeler 

2W 



le potentiel moyen. Il est égal à 



5 « 



Or W est aussi donné par l'intégrale 




mH^Mm-- 



étendue i\ l'espace entier. 

I.a partie étendue à l'espace intérieur est, d'après la formule 
de Green (4) » 




j4 masse homogène fluide 

Sur la surface, V est constant. L'intégrale a donc pour 
valeur 



VAVrfT, 



ce que Ton transforme immédiatement, puisque — AV = 4"^^ 
= 4?^, en 



AVrfT — -TT- I VAVrfT 








rfT+ 



VT 

" +W 



2 

5 



o o 

La partie de l'intégrale étendue au volume intérieur est donc 
h la partie étendue au volume extérieur, dans le rapport cons- 
tant — , en supposant évidemment le corps en équilibre. 

Mais l'énergie de ce corps est en raison inverse de sa capa- 
cité électrostatique. 

Nous voyons donc bien que la sphère est le volume pour lequel 
W est maximum. 



LA MASSE EST ANIMEE d'uN MOUVEMENT DE ROTATION. 



Formules générales. — Supposons maintenant le corps, 
homogène, soumis à une rotation w autour d'un axe fixe. 

Appelons J le moment d'inertie autour de cet axe. 



J= / rfT(a.'' + y'). 



MASSE EN ROTATION aS 

Si le corps change de forme, de façon que la projection du 
déplacement de l'élément de surface drj sur la normale soit Ç, 
J varie et on a 



rfj= / ÇrfT(j;2 4-y'). 



^2 

Nous avons appelé U la fonction \ -\ (-î^^ + y^)» 

W représente toujours l'énergie due au potentiel new- 
touien. 
On a 




= / UÇrfT. 



A la surface d'équilibre U est constant et égal à U^ 
On a donc 



(O^ 




iWH rfJ=U, I !:rfcr = U„rfT. 



D'ailleurs si le corps se déforme en restant semblable h lui-même, 

5 
W et J varient comme la puissance -rr- de T. 



On a donc 



rfW+- 


1 


dJ 


w+ 


■2 

0) 


J 



5 <fr 



3 T ' 



et par suite 



|-nT = W + — J=U. 
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On a 

AU = 20)2 _|_ ^Y^ 

== 2(0" 4^?> 

= 2(0^ 4'^» 

On a déjà vu (page ii) qu'une condition nécessaire pour 
l'équilibre est que 

2W' ^ 4^5 

c'est-h-dire que AU soit négatif. 

Donc U, à l'intérieur du corps, ne peut avoir de minimum, et 
comme sur la surface il est constant et égal il U^, à l'intérieur on 
a U > U,. 

Si AU est nul, U est constant et égal à U^. 

Enfin si AU était positif, U serait dans tout le corps inférieur 
à U,. 

Considérons maintenant l'intégrale 



Ué;T=2W-+--^j. 



Suivant que AU est négatif, nul ou positif, le premier membre 
est supérieur, égal ou inférieur à 



V,di = U 



A^ = uj=4-(w+4j). 



On a donc dans les trois cas : 



AU<o ^^W + -^jW 2W + -^J W >(oM 
AU = o A(w + -^j) = 2W+— J W=a>'-J 

O \ 2 / 2 



.2 



AU>o -|-(w + — j)>2W+ — J W<w^J 



MASSE EN ROTATION 2^ 

Limite de la vitesse de rotation, — Supposons maintenant 
que co varie d'une façon continue, La figure se déforme d'une 
manière continue, et on a 

2 ô 

l.a figure étant supposée en équilibre, W -] J est maximum 

ou minimum, en tous cas 

.2 



11 reste 



rfW + -^rfJ = o. 
2 



3 
5 



Donc 



d\] 



dtù 



' >o. 



Uç croit quand <o augmente. 

V étant le potentiel newtonlen, on a sur la surface 

Si donc la surface rencontre l'axe de rotation, U^ est le poten- 
tiel newtonîen au pôle. Si cj croît, ce potentiel croît également. 
Uq ne peut d'ailleurs dépasser 27:R^ (p^g^ ^6), R étant le rayon de 
la sphère ayant même volume que la masse fluide. 

En divisant membre à membre, on a 



w+-^j ^'o 



Si on fait croître w indéfiniment, w^ finira par dépasser ?:. 
W sera inférieur à w^J ; 

w'^J 2 



> 



2 
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et par suite 

2 diù d\i^ 



< 



3 co ^ U, 

Si oi croît indéfiniment, U^ devra croître indéfiniment ; mais 
Uq ne peut dépasser iizlA} ; il est donc nécessaire, ou bien que w 
cesse de croître, ou bien que la surface d'équilibre ne rencontre 
plus Taxe ; on a alors une figure annulaire. 

Nous verrons que parmi les figures d'équilibre, il y a des 
ellipsoïdes de révolution ou des ellipsoïdes a trois axes inégaux. 
Pour les premiers, to est inférieur h 4^ X 0,112. Pour les 
seconds, w < 4^ X 0,093. 

Il est aussi possible qu'il existe une succession de figures 
d'équilibre, pour lesquelles w passerait par une infinité de maxi- 
mums ou de minimums. Ici, le raisonnement ne s'appliquerait 
plus et w pourrait augmenter indéfiniment. 

Permanence de l'axe de rotation, — Jusqu'ici nous ne nous 
sommes préoccupés que d'un mouvement de rotation uni- 
forme. 

On. pourrait se demander s'il peut y avoir des figures d'équi- 
libre relatif pour une masse fluide animée d'un mouvement de 
rotation non uniforme. 

Nous allons prouver que c'est impossible. 

Nous pouvons d'abord supposer le centre de gravité du corps 
immobile; puisqu'il n'y a aucune force extérieure, son mouve- 
ment est un mouvement rectiligne et uniforme. 

L'équilibre ne sera pas troublé si on ajoute une nouvelle liai- 
son au corps, si on le solidifie. Alors le mouvement du solide 
autour de son centre de gravité est un mouvement a la Poîn- 
sot. 

D'après le principe de Dalembert, le travail virtuel de la force 
totale, pour tout déplacement compatible avec les liaisons, est 
donc nul. Dans le cas qui nous occupe, il n'y a qu'une liaison, 
c'est l'incompressibilité du fluide, s'il s'agît d'un liquide. Si l'on 
a affaire à un gaz en équilibre, l'équilibre subsiste a fortiori en 
supposant le corps incompressible. 

Si ox, 0?/, o.z représentent les déplacements virtuels d'un point 
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a', y, ::, lu condition d'incompressibilité s'exprime parla rela- 
tion 



d ^ , d ^ d 



11 peut y avoir trois sortes de déplacements virtuels : 

I® Les déplacements de toute la masse, considérés comme les 
déplacements d'un corps solide. L'application du principe de 
Dalembert montrerait que le corps doit avoir un mouvement à 
la Poinsot autour du centre de gravité, comme on l'a déjà vu. 

2** Il y a les déformations du corps. 

3** On peut supposer que la masse soit animée d'un mouvement 
tel que les molécules se déplaceraient sur des surfaces à densité 
constantes. 

Les surfaces d'égale densité étant aussi des surfaces équipo- 
tentielles, on aurait 

Op -. dp , dp , 



d.f Oy de 

Les surfaces d'égale densité restant les mêmes, il en est de 
même de la forme extérieure. Tout se réduit donc a des dépla- 
cements internes dont nous allons étudier le travail virtuel. 

Soit un point M de masse /?2, de coordonnées ^r, y, z, animé 
d'une vitesse x\y\d ^ l'accélération du mouvement étant x" ^ y'\^" * 
La force d'inertie a pour composantes mx" , ^^^l/' -> niz" , 

Appelons w^, w^, Wg les composantes de la rotation instantanée 
suivant les axes liés aux corps O.r, Oy, Oz, 

D'après les formules connues, on a 



x' 




"J3.'/ 




^2-. 


y' 




w,r 





(03^, 




— 


LÙ^X 




^iZ/- 



L'accélération a pour composantes 



f = [(^^z' — o)3.r') + {o}\z — coVr), 



3q masse HOMOGENE FLUIDE 

L'équilibre relatif ayant lieu, il y a équilibre entre le travail de 
rattraclion et la force d'inertie mx^' , mif , mz'' . 

Cette force se compose de deux parties : Tune est ce que Ton 

trouverait si Ton supposait le mouvement 
uniforme ; l'autre est due à l'accélération 
du mouvement de rotation. 

A un instant /, Taxe de rotation est OP 
(w^, (o^, 0)3) ; à l'instant t-\-dl Taxe instan- 
tané de rotation est OP' (fig. 5). Le vec- 

PF 

teur OQ équipoUent au segment 

tend, lorsque dt tend vers zéro, vers 
une position limite, et ses projections sur les trois axes sont 

{iù\, (0^2, W'3). 

Si on prend cet axe OR lui-même comme axe 0^ 




a){ = 10.^ = o. 



La force d'inertie due à l'accélération angulaire se réduit à 



0)3?/, — oy^z, o. 



Appliquons maintenant le principe de Dalembert aux déplace- 
ments virtuels que nous considérons. 

Le travail de l'attraction est nul, puisque la forme du corps n'a 



9 



pas changé. Le travail de la force centrifuge est égal à — oJ, en 

posant (0^ = Wi + io| -f- co^. 

Ce travail est nul puisque J ne varie pas. 
Le travail de la troisième partie est 



^'3 /. m{i)lx — xoy). 



Puisqu'il y a équilibre, il doit être nul. Or, on peut toujours 
choisir un déplacement virtuel tel que 



2_^m{ylx — xhj) 



ne soît pas nul; il suffit de supposer un courant autour de 0^. 
La somme aurait pour mesure l'aire du plan des xoy^ limitée 
par la projection du courant sur ce plan. 
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Il faut donc que (O3 soit nul, c'est-à-dire que Ton ait un mou- 
vement de rotation uniforme. C. Q. F. D. 

Il faut de plus que Taxe de rotation soit un des axes principaux 
d'inertie; cela résulte de l'étude du mouvement d'un corps 
solide. 

Stabilité de l'équilibre. — Jusqu'ici on ne s'est pas préoccupé 
des conditions de stabilité de l'équilibre. 

Lejeune-Dirichlet a démontré que la condition nécessaire et 



9 



suffisante de la stabilité de l'équilibre est que W + J soit 

maximum. 

Mais il faut faire la distinction entre la stabilité temporaire et 
la stabilité séculaire. Lord Kelvin qui le premier a fait cette dis- 
tinction, appelle stabilité séculaire ceEe qui a lieu même en 
tenant compte du frottement; tandis que la stabilité temporaire, 
qui subsiste quand on ne tient pas compte du frottement, ne 
subsiste plus si on en tient compte. Lord Kelvin a démontré que : 

La condition nécessaire et suffisante de la stabilité séculaire 
de l'équilibre relatif est la condition de Dirichlet. 

Evidemment cette condition est toujours suffisante ; mais sup- 

posons que W H J ne soit pas maximum, il pourra, d'après 

le théorème de L. Kelvin, y avoir équilibre stable, à condition 
qu'il n'y ait pas de frottement. Mais s'il y a un frottement, si petit 
qu'il soit, l'équilibre ne sera plus stable. 

On ne peut appliquer immédiatement le théorème de L. Kelvin, 
car il suppose que tout mouvement détermine un frottement, ce 
qui n'est pas le cas pour cette masse fluide : en effet, si la masse 
isolée dans l'espace se déplace d'un bloc, h la façon d'un corps 
solide, il n'y aura pas de frottement. 

Solide équivalent. — Avant d'aller plus loin, il est nécessaire 
de définir ce que nous appellerons le corps solide équivalent 
à la masse fluide. C'est un solide où, à l'instant considéré, 
les molécules ont la même position que dans le système fluide. 

La vitesse de son centre de gravité est la même que pour le 
iluide. 
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Les axes principaux d^inertie ayant évidemment la même posi- 
tion, les trois moments de rotation autour de ces axes sont les 
mêmes que pour la masse fluide. 

Son mouvement est donc bien défini a l'instant considéré ; 
mais il faut remarquer que le solide équivalent à Tinstant t, n'est 
pas le solide équivalent à l'instant t\ 

Théorème. — La force vive du fluide est équivalente à la force 
vive du solide équivalent^ augmentée de la force vive du fluide 
dans le mouvement relatif par rapport à des axes invariablement 
liés au solide équivalent. 

En effet, appelons Ç, t), Ç, les composantes selon ces trois axes 
de la vitesse absolue de la molécule m ; \', 't\\ Ç' les composantes 
de la vitesse de la molécule correspondante du solide équiva- 
lent ; Ç'', r/', 'Q' les composantes de la vitesse relative de cette 
molécule par rapport au solide équivalent. On a 

? = ?' 4- ?", 

■'i = -fi' + vi", 
Les trois forces vives considérées ont pour expression : 



T"=2-^(i"^+-/l"'+r), 



et il faut démontrer que 

T = T' + T". 

On a 



T = 2 ^[(?' + \"f + (r/ + y,")^ + (Ç' + 1:")^] 
= T' + T"+ 2 ;„(Ç'^'/ +.Y.^" + !;'Ç"). 
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Je vais démontrer que 



2 m{lV + ,S" + Ç'C) = o. 



Cette expression représente le travail d'un système de force S'', 
tel que la force agissant sur la molécule m ait pour composantes 
/??$'', wr/', mtl'y mais ce système est la différence de deux autres S 
et S', dont les forces sont m\, m't\y wÇ, et /wi\ mr\\ mX^ ; ces deux 
systèmes de forces sont équivalents. La résultante générale est la 
même et les moments résultants des deux systèmes sont les 
mêmes par rapport aux trois axes, par la définition même du 
solide équivalent. 

On a en effet 

puisque la vitesse du centre de gravité du solide équivalent est la 
même que celle de la masse fluide. 
De même on a 

puisque les moments de rotation sont les mêmes autour des trois 
axes de coordonnées. Donc S = S' et S" est nul; son travail l'est 



aussi. 

Par suite 



^7«(i1" + r.V + Ç'Ç")=o. 



Condition de stabilité de Lejeune-Diricblet. — Revenons 
maintenant h notre masse fluide. 
On a l'égalité 

T — W = C\ 

D'où 

T' + T^' — W = C\ 

Cette égalité suppose qu'il n'y a pas de frottement : autre- 
ment T'-}~T'' — W ira constamment en diminuant. 

PoiKCABÉ. — Figures d'équilibre. 3 
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Si on pose 

T' + T' — W = £, 

—7— est nul s'il n'y a pas de frottement. 

Mais alors c'est qu'il n'y a pas de mouvement relatif, et par 
suite T'' = o. 

Si T">o, -$-<o. 

Dans ce mouvement de rotation le théorème des aires s'applique, 
et on peut considérer comme des données les trois moments de 
rotation. 

Je vais démontrer que la condition nécessaire et suffisante de 
stabilité est que T' — W soit minimum. 

Pour une position d'équilibre relatif, on aura 

np; . 'p/ 

1 1 Q 

W = Wo 

T" = o. 

Si l'on part de cette position d'équilibre, et que l'on fasse varier 
très peu les éléments des mouvements des molécules, mais de 
façon que les moments de rotation ne soient pas altérés (c'est 
une restriction sur laquelle nous reviendrons), la masse fluide 
ne s'écartera pas beaucoup de la position d'équilibre. 

En effet, soient T' et W, les valeurs de la force vive et de 
l'énergie, on a 

T' - W > T'„ - W,. 
On a a fortiori 

T+r'_W = £ est d'ailleurs voisin de T^ — Wo- 

Cette quantité ne peut jamais croître, donc elle restera toujours 
voisine de T^ — Wo- 

Les valeurs des variables ne différeront que très peu de leurs 
valeurs pour Tq et Wo, donc il y a équilibre stable : la condi- 
tion est suffisante. 



DÉTERMINATION DE VAXE DE ROTATION 
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Elle est aussi nécessaire. En effet, si T^ — W'o n'est pas mini- 
mum, je puis choisir T' et W de manière que 



et même que. 



T'_W<T',— Wo 



T' + T"_W<T' — W,, 



T" étant suffisamment petit. 

Mais alors T'tJ- T'' — W ne peut que décroître ; T" lui, tendra 
vers zéro; T' et W ne pourront alors reprendre leurs valeurs pri- 
mitives. 

Il y a une restriction à lever. Supposons que l'on écarte le 
corps de sa position d'équilibre sans conserver la valeur des 
moments de rotation autour des axes principaux ; on sait alors 
que le solide équivalent aurait un mouvement où les moments de 
rotation auraient des valeurs peu différentes de celles qu'ils 
avaient précédemment et l'on est ramené au premier cas. 



L'axe de rotation est le petit axe de Vellipsoïde d'inertie. 
— Je vais démontrer maintenant que l'équilibve n'est stable que 
si l'axe de rotation est le 
plus petit axe de l'ellip- 
soïde d'inertie relatif à 
la masse fluide. 

Soient (fig. 6) Oxyz un 
trlèdre fixe de coordon- 
nées ; OA,OB, OÇ, les axes 
principaux d'inertie; w la 
vitesse instantanée de rota- 
tion dont les composantes 
sont ci)^, -Wj, Wg suivant les 
trois axes fixes,/;, q^r sui- 
vant les axes OA, OB, 
OC ; on a Fig- 6. 

Soit OR l'axe instantané de rotation, J le moment d'inertie 
par rapport à cet axe. 
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On a 

Soient u^, I^-ail^s» ^^^ moments de rotation par rapport aux trois 
axes de coordonnées. Le segment OM, dont les projections sont 
([jLp [x^, [X3), est fixe dans l'espace (théorème des aires), et Ton a 

Les projections de OM sur les axes OABC sont A/?, By, Cr, 
et l'on a encore 

jjL^ =Ay+BY+cv^ 

La condition de stabilité, comme nous l'avons dit, est que 
T — W soit minimum. Or, si on change l'orientation sans 
changer la forme du corps, W ne variera pas. Donc il faudra que 
T' soit minimum. OM est fixe dans l'espace. 

Si on pose 

rj v/B"= Y, 
/• y/c"= Z. 
On a 

T = — (X^ + Y^ + Z^) 

jji^ = (AX^ +BY^ + CZ^) = C^ 

Donc, le point X, Y, Z doit être sur un ellipsoïde donné et sur 
une sphère de rayon minimum. Donc Taxe instantané de rota- 
tion est le plus petit axe de l'ellipsoïde donné et par suite, de 
r ellipsoïde d'inertie. 

Si A est la plus grande des quantités A, B, C, l'axe de rota- 
lion sera l'axe des x. Alors on ap = to, q = o^ r=o. 

A = J, 

T' = ""' 



J • 



DÉTERMINATION DE LAXE DE ROTATION 3; 

Il faut donc que 

(i) w — r=w— ^^ 



2J 

soit maximum ; donc que 



8(w-4)=o 



Il sudTit aussi que 



(,) (w+ ^ ,) 

soit maximum (p. lo), condition qui n'est pas nécessaire. 
Donc 

,2 



«s 
Ù 



W + ^J=o. 



Les deux conditions donnent : 



►W + -^-p~ = o, 



8W+ — ôJ = o. 



On doit donc prévoir que si la condition (2) est remplie, la 
première le sera également. 

Supposons qu'il en soit autrement, c'est-à-dire que Ton puisse 
avoir en même temps 

/ 2 2 

' _-- CO _ ___ (1) _ 

w + — J<w, + — J. 

2 2 

Comme on a [x = tjJo, on conclurait facilement 

w„— ^<w-'''' ^" 



2 J • 
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En additionnant membre a membre, on aurait 



2 2 



(•j-j,)j<^(jj,-V), 



d'où 

J'— aJJ. + V<o, 

(J-J.)'<o, 

ce qui est absurde ; mais la réciproque n'est pas vraie ; la 
deuxième condition n'entraîne pas la première. 



CHAPITRE III 



FONCTIONS SPHÉRIQUES 



Expression de AV dans un système triple orthogonal, — 
Soient 

^(•^7 y, z) = Y, 

trois familles de surfaces, constituant un système triple ortho- 
gonal. 

Considérons le petit parallélépipède rectangle dont les faces 
sont les 6 surfaces. 

Les arêtes de ce parallélépipède ont pour longueurs 

ad^, bdr\^ cdX^, 
La diagonale est 

ds^ = a}d^' + Pdr^' + c^d'Q. 

Supposons une fonction V (1, v), 'Qj et cherchons à déterminer 
la valeur de T'intégrale 



/ 



AVrfT 
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étendue au volume du parallélépipède; elle est égale a 




étendue a sa surface. 

Le volume étant infiniment petit, la première intégrale a pour 
expression 

La seconde est la somme de trois intégrales étendues chacune 
aux deux faces opposées du volume. 
Pour la face $ = 5^, on a 

r=— hcdr.d'C, 

Pour la face opposée Ç = Sq + ^5> ^^ ^ 

et par suite pour l'ensemble des deux faces, ou a 



d / hc dV 
On a donc 



r ("v^f) ''■'"^''^^- 



\Nahc.d\ dr\ dX, 



,e, ,„rd ic ôV d ca ôV , d ab dV 1 



et finalement on a 

hc dV 



abc ^J ^l \ a ôÇ / 



Application aux coordonnées polaires. — Dans ce système 

on a 

a: == /• sin 9 cos o, 

1/ =r sin sin cp, 

c = /• cos 9, 
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et, en posant cos 9 = [jl, on a 



I — \L^ cos cp, 



y =^r y I — ^^ sin o, 

L'élément linéaire dans ce système de coordonnées a pour 
expression 

ds^ = dx^ + dxf + rf;:'^ = dr'' + /--^rfe^ + /-^ sin^ôrfcp^ 

= rf;.2 + ;.2 _^!^+ ,.2(,_^2W^2^ 

I JJL ^ 4/4 

Au lieu de Ç, y^, J^, nous avons /', [jl, cp et dans ce système 



a 
h 








v/i- 


■[^^' 


c 
abc 


— r\ 


-I-% 



AV 



On calcule alors racilement 
d^V . a dV . 



I d r 


■^^ da J 

I d^V 




' r\i u*) d'^' 



Polynômes aphèriques. Fondions sphèriques. — On appelle 
polynôme sphérique un polynôme homogène en x^ y, z satisfaisant 
a Téquation de Laplace, A P = o. 

S'il est de degré n en x^ y, z il contient ^ — coeffi- 

cients arbitraires. Mais il y a un certain nombre de relations 
entre ces coefficients puisque l'on a AP„ = o. Cette relation est 

de degré n — a. Il y a donc n • relations entre les coef- 
ficients. 
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Il y a donc 

(/i+i)(/^+2) _ n{n—i) ^ 2 ^_^_ j 

2 2 

polynômes sphériques, indépendants du degré n. 

On peut exprimer a;, y, ^ en coordonnées polaires : on aura 
alors 

p ,.n Y 

Y„ étant uniquement une fonction de 'f et de [x, que nous 
appelons fonction sphérique d'ordre n. 
On a évidemment la relation 

Considérons maintenant l'égalité (2) démontrée page 4 

Supposons que la surface S soit une sphère de rayon i, que U 
et V soient deux polynômes sphériques d'ordre m et n, le second 
membre sera nul, et on aura 





7 " ""ôT^ "^ 






Donc, en tenant compte de l'égalité (i), on a 




[m — n) I P„j P,j rfo" := o. 



Si m est différent de n, on a 



PmPn^^= O, 
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1-intégrale étant étendue a une sphère de rayon i ; et on peut 
encore écrire ceci : 

Y„jYnrf'fi[X = 0. 



t/-l Jq 



Si /w = 72, ce résultat n'est plus vrai. 

Fonctions sphériques fondamentales. — On a vu que Ton a 



x = r \J \ — iji^ cos '^ = r y I — 



1 e'+e-' 



tS — Ï3 

y = /• y I — [JL^ sin '^ = 7' y I — ^ — 



11 



Si on remplace :r et y par ces valeurs dans P„, on aura 

P„ = V ;•» ey'?Xï, 



p=_n 



XJJ étant une fonction de [jl seulement ; ce sera un polynôme 

en jx et en \J \ — [I^. Les polynômes oùy; est pair, ne contiendront 

y I — |JL* qu'à des puissances paires. Les polynômes où p est 

impair, ne contiendront y/ i — jx^ qu'a des puissances impaires. 
Les premiers seront donc des polynômes en jjl, les seconds 

seront des polynômes en [x multipliés par y \ — [x^. 
Substituons la valeur de P dans l'équation A P„ = o. 
On obtient 

- I — o. 
I — tx J 

En annulant le coefficient de /" - ^ e'^'?, on obtient l'équation 
à laquelle doit satisfaire XJJ. C'est 

rf[X L à^x \ I — [X' 
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L'intégrale générale de cette équation est une fonction trans- 
cendante, mais il y a une intégrale particulière qui est de l'une 
des deux formes exigées par la définition de XJ, et qui est 
déterminée à un facteur constant près. 

On voit que Xg et X~^ satisfont à la même équation. On a 
donc 

A étant une constante quelconque. On a ainsi la forme des 
2/? -|- I fonctions sphériques cherchées. 

On peut écrire autrement 2/i + i fonctions indépendantes en 

posant 

XJ, Xicoscp, Xf,cos2îf ,..., Xjcos/iy, 

X^sincp, X^sin 23 ,..., Xj sin /2'j. 

Deux fonctions sphériques diffèrent donc soit par /i, soit par p^ 
soit que n et p soient égaux, mais que l'une contienne le sinus, 
l'autre le cosinus de p'^. 

On a déjà démontré que 




est nul si m^n; nous allons achever la démonstration et démon- 
trer que l'intégrale 

YTdz> 




est nulle si Y et Y' sont deux fonctions différentes. 

Il suffit de le démontrer dans le cas de deux fonctions de 
même ordre. 

On a alors 



Xj; XJ a[A I cos p'f cos çrcp rf'p 



1 .7 
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OU bien 



41 /**« 

Xj X^ d]k I sin /7'^ sin qrf rfy, 

-1 .7 

Xj^ ^/[x I sin/;'^ cosp'fdo. 

-1 •/ 




Dans les trois cas, la seconde intégrale est nulle. 
Enfin il est évident que 



rf(x / Y'df 





est différent de zéro. 

Détermination des polynômes X};. — Proposons-nous maîn- 
Icnant de déterminer XJ, et pour cela, commençons par le cas 
où p = o. 

L'équation est alors 

C'est l'équation qui définit les polynômes de Legendre, et 
on a 

X^„ = -f— (l— UL^K 

" du/" ^ ^ ^ 

Xj (Cos 9) est donc une fonction sphérique où 9 représente 
l'angle d'une direction avec l'axe 0;:;. 

Par raison de symétrie, Xj^ (Cos e), où t représente l'angle de 
la direction considérée avec une autre direction quelconque, mais 
fixe, sera aussi une fonction sphérique ; mais on a 

cos e = cos 9 cos 9^+ sin 9 sin 9' cos (*^ — '^j 
et cela, quels que soient 9 et 0'. 
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Donc 

X^ I cos cos 0' + sin sin 6' cos (© — cp') 1 

sera encore une fonction sphérique. 
Si on pose 

sin 9' e-'î' = 2$, 

sin ^' e^'i' =:. 2rf\^ 

d'où on conclut, en ajoutant membre à membre après avoir mul- 
tiplié ces deux équations respectivement par e*? et e~ '?, 

sin 6' cos (cp — o') = Çe'î +7,e~'? 

et en multipliant membre à membre, 

sin2 6'==4ÇT, ; 

d'où 

cos-9' = 1 — 4?''! > 

la fonction 

est donc encore une fonction sphérique quels que scient \ et y. 
En particulier, si on fait 

est une fonction sphérique. 

Si je développe suivant la formule de Taylor, j'ai un poly 
nome : 

X« [[x + v/7^^e'>1= X„(a) + e-VT^Il?"-^ + ... 

1,2... JJ duJ' '" I . 2. . ./l Û^JJt," 






m # 
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Chacun des termes de la somme est une fonclion sphérique. 
On a donc défini les fonctions X,^ à un facteur près. 



Propriétés des fonctions spbériques. — Soit U une fonction 
de /•, et P„ un polynôme sphérique : 

Proposons-nous de calculer A (U P„). On a 



mais 



" ' do; clr ' àx dr ' /• 

A (UP„) = P,|-^ + - ^J 

2 rfU / dP„ OP,. dP,. \ 

t, parce que P'* est une fonction homogène, 

_ r rfnj , 2(«+i) rfU -] 

Proposons-nous maintenant de calculer 

a(uy„)=a[-^xp,.] 

puisque P„=:/"Y„. 
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Couche sphéroïdale. — Nous allons appliquer ce qui précède a 
Tétude d'une masse fluide sphéroïdale, c'est-à-dire d'une masse 
fluide dont la figure d'équilibre est peu différente d'une sphère, 
ou, d'une manière plus précise, dont la surface est comprise entre 
deux surfaces sphériques de rayons voisins. Une couche sphéroï- 
dale sera une couche comprise entre deux surfaces sphéroïdales. 

Action d'une couche sphéroïdale, — Soit(fig. y) rfo* l'élément 
de surface d'une couche de rayon i, s l'épaisseur de la couche 

sphéroïdale, d'z l'élément de volume 
de la couche sphéroïdale ; on a 

On peut supposer que la masse 
contenue dans rfr est sur l'élément 
rfo-, et est représentée par ikd^r 

Soit un point M dont les coordon- 
nées cartésiennes sont x, y, z^ les 
coordonnées sphériques 

/•, [JL=:COSO, cp. 




Fig. 7. 



Soit un point M' dont les éléments analogues sont x' ^ y\ z! \ 
/•' = I , (jl' = cos 6', ç^ 

La distance MM' est égale a i + /-^ — 2/- cos (OM, 0M% 

Cos Y = cos (OM, OM') = cos cos 9' + siu sin 6' cos (3) — d) . 

Développement en série de- ^..., . — Le potentiel en M, du à 
la couche sphéroïdale, est 



V = 



e'rfd' 



MM 



/ t 



s' est une fonction de [jl et de cp ; /• est plus petit que i , 
Or 



MAP :^ I 27- cos Y+ /'^ = (l /'C'Y ) (l /'C-^V), 



w=(-"")' ' ( 



I — /-e-'Y 
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On a 



(. \ _L I 1.3 o a 
I — re'ï)" 2 = iH re'T^ ^/.V'v+... 
/ 2 2.4 

1.3.3... (2/1 — i) 

2.4.0... 2/1 

C'est une série absolument convergente pour r < i . 
De même 

/ A--!- . ï .-1 . 1.3... (2/1 — 

( I — re-«T ) 2 = I H /•e-'ï+ ... H 7^^ ^ /-"e-^'ï + . . . 

\ / 2 2.4... 2/1 

Le produît de ces deux séries est une série encore conver- 
gente, si /• < I ; le coefficient du terme en z*** est 

. 1.3... (2/1 l) . , 1.3... (2/2 3) I . ax; 

2.4 .. 2/1 2.4... (2/1 2j 2 

, 1.3... (2/1 — 5) 1.3 , ... , , 1.3.5... (2/1 — i) 

2. 4**- (2^/2 — 4j ^'4 2.4.0... 2/1, 

Le coefficient du terme en e'"ï égale celui du terme en e"^"^ 
On a donc un polynôme h coefficients réels du cosinus de 
l'angle y et de ses multiples jusqu'à n ; comme les coefficients de 
ce polynôme sont tous |)ositifs, ce polynôme atteindra sa valeur 
maximum pour y== o. Mais alors : 



- = = I + /• + ... + / 

' I — /• 



M. 



MM 

Donc les coefficients du développement sont inférieurs ou égaux 
h l'unité. 

On peut donc écrire 



[M' ""2j^"''"' 



MM' 
A„ étant une fonction de [x, de îp, de \k' et de cp'. 

Développement en série du potentiel d'une couche spbéroï- 
dale, — Le potentiel V sera alors égal à 



^z-" / A„£Wt. 



PoiNCARÉ. — Figures d'équilibre. 
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Les termes de V sont moindres que ; 



L'erreur commise quand on prend pour V la somme des n pre- 



•» I t'dfj. 



miers termes, est moindre que 






La série sera donc convergente si 7*<i. 

* I 
rrjrrjpeut ètrc dévcloppé suivant les puissances croissantes 

des variables x, y, Zy 



MM' 



y^ (x' — xY+(jy' -yy + [z' — zf 

et, si on revieat aux coordonnées polaires, on aura un développe- 
ment suivant les puissances de /•". 

Si j'appelle V„ l'ensemble des termes de degré n dans V, j'au- 
rai 

v=v,+v.+...+v„ + ... 

Or V satisfait à l'équation de Laplace : on a donc 

AVo + AV,+ ... + AV„ + ...=o. 

Cela exige que chacun des polynômes V^ soit un polynôme 
sphérique; et, si on remplace dans V„, x, ?/, -c, en fonction de r,[jL 
et îp, on aura un résultat de la forme z" Y„ où Y„ est une fonction 
sphérique. On a donc 



=y'-v.. 



Cette formule est valable si /• est plus petit que i . 
Supposons maintenant r supérieur k i, et considérons le 
point Mj image du point M, tel que OM^ X OM = i (fig. 8) ; 



DEVELOPPEMENT EN SÉRIE DU POTENTIEL 



5i 



j'ai OM =i r, OM^ =— . Si M' est un point de la sphère, on a 
les égalités 

MlVr _ OM _ OM' \/0M X OM' 



M'M, OM' OM, 
Or le potentiel en M, est 



y/OM'xOM 



M' _./0M 

m; ~ V OM, - '• 



On a donc 

* V 
/' 

Or le potentiel en M,, point 

situé à une distance — du 
centre, est 






2 







M 



Fig. 8. 



Le potentiel en M est donc 



V Yn 



Supposons maintenant que M et Mj se rapprochent de la 
surface de la sphère en restant inverses Tun de Tautre. Les 
valeurs du potentiel M et M, sont voisines, mais les dérivées 
ont des valeurs différentes de part et d'autre de la surface et 
on a 

dV dV, 

O/' Q1\ 



mais 






V. = Vr, 



V;-2_ 



dr 



,.3 
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On a donc lorsque /* tend vers i 

Dans les mêmes conditions, puisque — :r — est fini, on a 

lim -:r (/• I) =0, 

enfin 

lim V (i — r^) = o. 

Si on ajoute ces quatre égalités^ on voit que l'on a 



Posons 



liml 2r -^+ V J= 4-6. 






Si /• est inférieur à i, <[> peut se représenter par une série de 
la forme 



*= 2 (2«+i) /•»¥„; 



^ est donc une fonction satisfaisant h l'équation de Laplace, 
puisque c'est une somme de polynômes sphériques. 

Si r tend vers i, ^ tend vers ^izt. Si on démontre que ce 
développement est encore valable pour /•= i, on arrive au théo- 
rème de Laplace : 

Théorème. — Une fonction quelconque de deux variables est 
développable en une somme de fonctions sphériques. 

Potentiel en un point de la couche spbèroïdale. — Soit e la 
fonction donnée sur une sphère de rayon i. 
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Elle définît à Tîntérieur de la sphère un potentiel V, prenant 
sur la sphère la valeur ^ize. 

Ce potentiel est développable en une sérîe entière 



(0 2^-'"" 






convergente si /• est inférieur a i ; et, si on regarde r comme une 
variable complexe, cette sérîe définit une fonction $ convergente 
à l'intérieur du cercle de rayon i, et qui sera obtenue en faisant 
dans V, r=pe*ï. 

Dans V, faisons r= e*»; on a une fonction ^' définie sur le 
cercle de rayon i, sauf peut-être en un certain nombre de points 
singuliers. Si la série (i) converge pour p = i, elle représentera 
<ï>' d'après le théorème d'Abel. 

Je dis qu'il en est ainsi dans le cas où ^' est développable en 
série de Fourier, 

Nous avons alors 

Je dis que l'on a 

A,;=A„ B,;==o. 

En effet, d'après le théorème de Gauchy, la suite des valeurs 
^' définit à l'intérieur du cercle de rayon i une fonction qui 
n'est autre que la fonction V (pe'-). 

Cette fonction est développable suivant la série (*i), et on a 
d'autre part 

2m J r''^^ 2.nz J ^ 

les deux intégrales étant prises le long du cercle de rayon i. 
Mais la première intégrale se réduit à A'„, la seconde à B'„. On 
a donc 

a;=a„, 

B;=o. 

Donc, SI 4>' est développable en sérîe de Fourîer, la série sera 
encore convergente pour v=ze'^ et en particulier pour /•= i. 
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Tout revient donc h démontrer que '^' est développable en 
série de Fourier, c'est-à-dire que $' est une fonction de cp admet- 
tant une dérivée, sauf en quelques points singuliers, et qu'en ces 

points singuliers l'intégrale / | ^' | rfy a un sens. 

Nous allons maintenant démontrer que e est développable en 
une somme de fonctions sphériques. Le potentiel en un point M 
de coordonnées (/•, 8, y) est 



/ MM' ' 



e' étant la densité de la couche sphérique, au point (6',<p'). 
On a 






d/- / MM'* dr 



/ 



/ e' (cos Y — ;•) d's' 



Gomme cela a été démontré, nous avons 



e = *r9.ol = 



Supposons que OM soit une direction fixe que nous allons 
prendre comme axe Oz \ 

6 = 0, [/.= !, 
MM' =1 — 27* cos Y + /'^ = I — 2'xV -f- ''^• 

L'élément de surface est d's' =d^'d'o\ et on a pour valeur en M 
de la fonction ^ 

— — r / ''''?'' 

1 (,_2r!.' + r')Vj, 
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car e' seul dépend de 'y. Si on pose 
on a 




_i (i — 2r[x' + 7'2)' 



Supposons que la sphère soit partagée en un nombre fini de 
polygones sphériques dont les côtés soient des arcs de courbes 
analytiques, et que dans chacun de ces polygones e' soit une fonc- 
tion continue de \»! et de cp'. 

Dans ces conditions, F([jl') sera une fonction analytique de [x', 

excepté pour les valeurs de \k' correspondantes aux sommets des 

polygones, ou aux parallèles de la sphère tangents aux côtés des 

F (uO 
polygones. La fonction '3 ■ sera alors intégrable entre — i et 

Si on pose /• = e*t', on aura 

MM'=2r — lr-\ ;-] — \k' =2r| cosij/ — ^'\ 

et si on pose ensuite 

cos^ — 1^'= 2^> 
on aura 

^"' F [^') d^! 






z^ 



Nous démontrerons plus loin que $ est une fonction de i^ 
admettant une dérivée, et est par suite développable en une série 
de Fourier. 

Enfin, si nous développons suivant les puissances de r, nous 
aurons une série convergente pour)7*|<i. Elle sera sûrement 
convergente même si | /• | = i , dans le cas où s ne comprend qu'un 
nombre fini de termes. 
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Voyons quelques exemples qui nous iservîront par la suite. 
Si on fait c' = , on aura 

271 



F (fx') = 



/2k 



si on fait e.' =^ — , on aura 

71 



On pourra prendre pour fonctions sphériques correspondantes 
Y.=-^ouYo = — . 

27Î 7t 

u! 
De même, si e' = —^—, on a en appliquant la formule de la page Sa 



Dans le cas qui nous occupe 

Yo = Y2 = ...==:o 

et 

$Q=2/* = 2e*t'. 

Théorème de Laplace, — On a toujours le droit d'écrire 

a et i étant des constantes quelconques ; <P sera alors de la forme 

2a+2ie*^' + ^^ 
et si on pose 



2 /-g- \/r 



on aura 
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Ceci posé, pour établir que O est développable en série de 
Fourîer, il faut démontrer : i** que ,. - existe et est fini sauf en 
un nombre limité de points singuliers de la fonction $ ; 2** que 



|<^| (2^ est fini en ces points singuliers. Il y a une difficulté, 



car Z s'annule si [x' = cos ^. Il faut donc que F (jjl') s'annule pour 
[x' = cosij/ et alors on doit étudier vers quelle limite tend la 
fonction lorsque [x tend vers cos ^. 

I® — 77— existe. En effet on a 

Or, je puis disposer de a et A, de façon que l'on ait 

Fj (cos ij/) :=: O, 

F( (cos l{^) = O. 

Si donc cos ^ — y/ est un infiniment petit du premier ordre, 
F(jx'), qui est analytique au voisinage de cos ^, sera un infini- 
ment petit du second ordre, et, sous le signe d'intégration, les 

3 5 
dénominateurs sont respectivement d'ordre — et — ; le nunié- 

rateur est d'ordre 2. L'intégrale est donc finie, sauf si cos ']^ est 
une valeur singulière de F([ji.'). Il faut qu'en ces points singuliers, 
où ^j; = <{^ç, l'intégrale 




soit finie. 
Supposons 

a = F (cos ^ J 
b =0. 
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Oh a alors 

Le second membre doit rester fini. La première intégrale est 
finie. 11 faut démontrer qu'il en est de même de la seconde. 

C'est une intégrale double qui, aux points singuliers, doit être 
un infiniment petit d'ordre inférieur h 2, et aux lignes singulières 
d'ordre inférieur à i. 

Or, nous avons ici une ligne singulière, la ligne [jl' = cos i|/, 

3 
mais pour ^ = cos i^^ le dénominateur est d'ordre — > le numéra- 

teur est d'ordre i, sauf aux points singuliers de F([x'). Donc 

l'élément dç l'intégrale devient infini de l'ordre — , l'intégrale 

a un sens. 

Aux points singuliers F(ijl') z=I=. F(cos ^\ la fonction est d'ordre 

3 

— , l'intégrale a encore un sens. 

La fonction ^ remplit donc les conditions nécessaires pour 
êtr^ développable en série de Fourier. 

La formule de Laplace est valable pour /•= i, et l'on a 



^2< 



2W+l) Y„. 



Le théorème de Laplace est ainsi démontré. 



CHAPITRE IV 



MASSE FLUIDE HETEROGENE. PROBLEME DE GLAIRAUT 



Préliminaires. — Supposons le mouvement de rotation de. la 
masse sufSsamment lent : la surface différera peu de celle d*une 
sphère, le carré de la vitesse de rotation to^ sera considéré 
comme un infiniment petit du premier ordre. 

On a déjà vu que si le mouvement de rotation est nul, les 
surfaces équipotentielles sont également des surfaces d'égale 
densité. Le potentiel V^ en un point est uniquement fonction 
de r et l'équation de Laplace se réduit à 



(0 






dV 



dr 



f = - 4«?- 



La masse contenue entre deux sphères infiniment Toisines de 
rayon r et r-\' dr est 



et la masse totale contenue à l'intérieur de la sphère de rayon r 

est 




Tzr^pdr. 



Si on appelle D la densité moyenne de cette sphère 




4 



^Tzr'pdr = -5- Tt/'^D. 



9 » 
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D est une fonction de r dont la dérivée D' est négative dans 
Thypothèse d'un équilibre stable. 

En différentiant cette dernière équation, on a 

En supprimant le facteur /\Tzr^di\ on a 

(2) 3p=3D + rD'. 

En prenant encore une fois la dérivée, on obtient la rela- 
tion 

(3) 3p' = 4D' + /D" 

qui nous sera utile par la suite. 

Faisons encore la remarque suivante. On a évidemment : 



3D ^ 3D + /D' ^ o. 

La première inégalité est vérifiée puisque D' est négatif; la 
seconde Test, car p est positif. On a donc 

o^-^^-3. 

On a -^— =0 si la masse est homogène. 

On a -r— = — 3 si toute la masse est concentrée au centre. 

Le potentiel de la sphère de rayon /• en un point de sa surface 
est 

La force qui s'exerce sur une molécule de la surface est 

Puisque p et V^j sont uniquement fonction de /*, on peut les sup- 
poser fonctions l'un de l'autre, et, en vertu de la relation (i), 



DEVELOPPEMENT DU POTENTIEL EN SÉRIE 6i 

AVq est fonction de V^j, et on peut écrire AV^=/][Vq); en diffé- 
rentiant 

^V, ^V. Ç, _4„,D 

En tenant compte delà relation (3), on a 

4D' 4- rD" 



(5) r (vj = 



/•D 



Développement du potentiel en série, — Supposons main- 
tenant un mouvement de rotation lent. La surface se déformera, 
mais peu. Le potentiel en un point d'une sphère sera dévelop- 
pable en une somme de fonctions sphériques, et Ton aura 

Y représentant une fonction sphérique et H étant un coeffi- 
cient qui est une fonction de /• de Tordre de w^, puisque V dif- 
fère peu de Vq. 

L'équation de Laplace peut s'écrire 

(6) AV=AV,+^A(HY) = -4î:p. 

On a d'ailleurs, d'après une formule démontrée plus haut 

(P- 47) • 

en supposant que n soit l'ordre de la fonction sphérique Y, 
Si on pose 

U = V + — (x-^ + i^'), 



la condition d'équilibre est, comme on Ta démontré plus haut, 
que les surfaces U = constante coïncident avec les surfaces 
d'égale densité. 
On a donc 

AV = — 4^? = fonction de U, 



» » 
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iù^ étant très petit, V diffère peu de V^ ; 

AV diffère peu de AV,,. 

La relation entre AV et V diffère peu de la relation entre AV,, 
et Vq, et on peut par suite écrire 

AV=/-(U) + ?(U), 

'j étant une fonction dont la valeur est constamment de Tordre 
de a)'^ 

Si Ton développe le second membre suivant la formule de 
Taylor, on a 

AV = /(V,) + (U-V.)/-'(V,)+... 
+ 'f(Vo)+(U-V,)-/(V,) 



En vertu des hypothèses faites, on peut négliger les termes 
infiniment petits du deuxième ordre et écrire 

AV =/-(v,) + (U - V,) f (VJ + ? (V,) 

= AV, + (U-V„)/'(V„) + <p(V„). 
Enfin, en tenant compte de (6), on peut écrire 

2a(HY) = (U-V.)/'(V„) +«p(V„). 

On a d'ailleurs 



w^ 



{x^-^y^) est lui-même décomposable en une somme de 
lonctions sphériques 



w^ 



, (^^ + |/')=2CY. 



9 9 O 



{x^+y') = — [x' + f + z^ + :^ (^» + y» _ 23») 
w^;-'*' (ù^i^ x^ -+- y^ — 23^ 



3 ' 6 /•=' 
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Le développement se réduit à la somme des deux fonctions 



sphérîques Y = i et Y= =^^2 ~~ affectées respective- 

ment des coemcients Cj ^ — 5 — et C^ ^ 



3 ^ 6 ■ 

Les autres coefficients sont nuls. 
On peut dans ces conditions écrire 



U 



-v.=2("+<^)^' 



^A (nY)==2(H + C)Y/'{VJ + «(V.). 

Détermination des coefûcients du développement — Pour 
trouver les valeurs de II, on peut égaler terme à terme les coef- 
cients d'une même fonction Y. 

A(HY) 
contient Y en facteur ; on a donc 

i(IIÏ)=Y(IH-C)/'(V.) 

Si on prend Y = i , on a 
(S) ,n=(H+^rf^^+,m. 

H est donc déterminé par une équation linéaire du second 
ordre. 

Si on prend 



Y = 



r' 



H est déterminé par l'équation : 

wV'\^ /•D"+4D' 



(9) A(HY)=(lI+ ^)y 



rD 
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Pour les autres, on a simplement Téquation 

f 

(.0) A (HY) = HY -:5^^±^. 

Ces dernières équations ont une intégrale évidente, H = o; 
nous allons démontrer que Ton ne doit pas en considérer d'autres. 

Dans le cas où Y est une fonction du premier ordre, on peut 
trouver la valeur correspondante de H. En effet, il y a trois fonc- 
tions indépendantes du premier ordre : 

X y z 



> 1 

r r r 

Considérons Tune d'elles, — par exemple, on a vu que Ton 

avait 

AV.=/'(V„). 

Si je prends les dérivées des deux membres par rapport à .r, 

• » • 
jai: 

^ (^v„)=-^/(\;), 



dx ^ "' dx 
dV„\ dV 



K-d^)=/^(^«) 



U 



àx ' 



A 



\ dr ' àj;)~' ^ '^ dr ' àa:' 



j7 

On voit que si Y est égal a —, l'équation 



MHY) = .n'i£:+i2: 



est satisfaite pour II = a — -j-^ , 11 en est de même si Y = — 

* ar r 



ou 1 = — . 



LES SURFACES DE NIVEAU SONT DES ELLIPSOÏDES 

On a alors pour les termes du premier ordre 
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rfV„ I ax-\~ hy + cz 



dr \ 



r 



). 



a^ b^ c étant des constantes quelconques. — Nous verrons pro- 
chainement qu'on peut les supposer nulles, à condition de 
prendre le centre de gravité de la masse comme origine des 
coordonnées. 

Les surfaces de niveau sont des ellipsoïdes. — Considérons 
une surface de niveau primitivement sphérique, lorsque la rota- 
tion est nulle. 
Pendant la rotation (fig. 9), 
' le point M est venu en M' sur la 
nouvelle surface de niveau. Le 
déplacement normal est 

Ç = MM'cos(/-, MM'). 

Avant la déformation, on avait 
en M 




t'ïg- 9- 



AVo=— 47îp. 



Après la déformation, la densité en M' est encore p. La dcn- 

AV 
site en M est 0^= — -- — - , en appelant V le potentiel en M, et 

on a la relation 



AV = AV.-!; 



dr 



On a de plus 



rfAV, rfAV„ rfV„ 



dr 



dV, • dr ' 



et, en se servant des relations (4) et (5), 



AV = AV.+ -^!:(4D'+;-D"). 



PoiNCARÉ. — Figures d'équilibre. 



5 



tii',. 



• # 
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Nous avons crailleurs déjà écrit 

av = av,h-Vaiiy. 

On a donc 

Revenons maintenant aux équations déterminant les coeffi- 
cients II. 

L'augmentation de volume entre la figure d'équilibre au repos, 
et la figure d'équilibre dans le mouTement de rotation, est nulle. 
On a donc 



ou 




AHYr/cr = o 



Ceci peut s'écrire 



2'f('-) / Vrf^=o. 



Chaque terme s'annule séparément, sauf toutefois le premier, 
pour lequel Y est égal à l'unité. On a donc, en désignant le coef- 
cient de Y^ = i , par 11^^, 

A(HJ = o, 
mais IIq est uniquement fonction de R si bien que 

dt^ ^ r dr 
I/intégrale générale de cette équation est 

II„ = A+-, 

A et B étant deux constantes. 



APLATISSEMENT * 6; 

Maïs remarquons que V = Vq+ SIIY; à l'infini, V doit être 
nul comme V^^. Donc A est nul. De même B est nul ; autrement 
le potentiel au centre serait infini, ce qui est absurde. 

i/équation (8) qui définit H se réduit dans ce cas à 

/•'(V„)-^ + ?(V„) = o. 

Cette équation définit o (VJ. 

Il nous reste encore à voir que pour les autres valeurs de Y, les 
seules valeurs possibles de H sont les valeurs II = o. Si nous 
admettons ce résultat, nous voyons que î^ se réduit à un seul 
terme, et que Ton a 

t = — ^^^^ • =YiV> 



On a donc 



^_ ^'+r^ ''% ('•). 



Ceci montre que les surfaces de niveau sont des ellipsoïdes. 
En effet, on a 

{^' +f+ z'] = (/• + î:)^ = ,^ + 2,-:; + ç^ ; 

mais nous pouvons négliger JJ^, et il reste 

C'est Téquation d'un ellipsoïde de révolution autour de Taxe 0". 

Aplatissement, — L'aplatissement est le rapport de lu 
différence entre le rayon équatorial et le rayon polaire, au 
rîiyon moyen, que nous pouvons supposer égal à r. 

Pour avoir ces rayons nous pouvons calculer !J pour les valeurs 
(.r = y = o, 5 = /•) et pour {y = z = o, :r = /•), on trouve les 
deux valeurs respectives 



^-=^ 2/-'' 

V 2 f 

V = /y/-. 



'■^ ('•) 



» m 
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La différence est 

et raplatissement est 

e = J-^-L- = 3rJ; (r) 

__^_M__ — il 

en désignant par Y la fonction sphériqué 

X- ■+- If 2?- 

On a vu que 

3 SAIIY 

^""47: /•D" + 4D'' 

La somme se réduit ici à un seul terme. L'équation (g) donne 

MHÏ) = (H + 4i)Y/-W = (H + 41)Yi2:+!£: 

Donc, 

47t/-^DeY = (H+-Î^W, 

Ae/-^DY= -|^[aHY + -^A/'^y1. 

On a d'ailleurs aussi 

A/-^Y = o, 

puisque Y est une fonction sphériqué du second ordre. 
On a aussi 

X AHY = (/D" + 4D') e/ Y. 
4^: ^ ' 

L'équation devient donc 

Ae/-=DY = (;D"+ 4D') e/T. 
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On a démontré que si P est un polynôme spliérique d'ordre /?, 



on a 

U' 



A(up)=pru"+2(/i+i)^i. 



SI nous posons pour /i = 2 

U=eD, 

p=/n, 

On aura Téquation 



éPieD) . 6 rf.eD ,^„ \\)'e 

I 

en développant les calculs on trouve, 

e"ly + ae'D' + -1 eD' + — ^''D = o. 

C'est une équation linéaire du second ordre que Ton peut 
ramener à une équation du premier par le changement de 
variables. 

re' 

^=— ' 

, re" re'^ , e' 

e e e 

'^ / , 2 



L^équation devient 

(12) D(/V4--/l' + 5-/i)+2/-D'(i+70 = o. 

^^'^5^ V équation de Clair a ut. 

Limites de tj. — Nous avons vu (p. 60) que Ton avait 

— 3<^^<o. 

Dans Téquation (12), substituons ces valeurs extrêmes de rl)\ 
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et considérons les deux paraboles (fîg. lo], oâ l'oa a porté en 
iibscieses 1^ et en ordonnées jt,' 




Le point dont les coordonnées sont /-r/et r, ne peut se trouver 
qu'entre les deux paraboles (dans la région hachurée sur la figure.) 

Je vais maintenant démontrer que r, est compris entre o et 3. 

r, ne peut être négatif. En eilet, s'il est négatif pour une 
valeur i\, on pourra trouver deux nombres positifs x et p tels que 
l'on ait 



c et ^ pourront même satisfaire ii riuégalité 
0<a <2 <5 <p. 
Je dis que t^ ne peut dépasser — a. 



UMITES DE T^ -X 

■ En effet, si r, atteignait le valeur — a pour une valeur /• infé- 
rieure à t\^ —T-!- serait négatif, /•— 7-^ serait négatif. 

Le point dont les coordonnées sont /-r/ et r^ = a ne se trou- 
verait pas alors dans la région hachurée. 
De même Tj ne peut dépasser — p. En effet, si pour une valeur 

d'i . . , vd'f 

de /• inférieure ii/' on avait r < /• , -~ serait positif, • — 7-^ serait 

" dr ^ dr 

positif aussi et le point représentatif serait dans une région 
dii plan où le point (/t/, t^) ne peut se trouver. 

Donc si pour une valeur de r,,, v^ est négatif, il est, pour toute 
valeur de /'inférieur à r^, compris entre deux nombres — a et — [4, 
on a donc 



re^ 



< — a; 



e' a 

— < 

e r 



En intégrant de r^ à r </o, on a 



L^.>L(^)-, 



^0 



e > A/'"" «. 

Or e croîtrait plus que lorsque /• tend vers zéro, ce qui 

est impossible pour /• suffisamment petit, car e serait très 
grand. Donc y) ne peut être négatif. 

Je dis maintenant que */; est inférieur à 3. Le raisonnement est 
le même. 

En effet, supposons que pour ;•:=: r^ on ait r^ = r^^ > 3, je dis 
que pour /• inférieur à r^, on aurait encore r, > 3. 

En effet, l'inégalité ne cesserait que si 7; devenait inférieur 
à 3, ou si r, devenait négatif en passant par l'infini. Cette dernière 
hypothèse est à écarter d'après ce que Ton a vu plus haut. Il 
reste donc a montrer que r, ne peut devenir égal à 3. En effet, 
pour r, = 3, r/ serait positif, rr/ serait positif, ce qui est absurde. 
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le poînt (/-/i', r\) ne serait pas compris entre les deux courbes. 
On a donc toujours 

y; <3. 
On a d'ailleurs 

/•y/ -j- v;^ — 71 — 6 < o pour r, > 3. 
A fortiori on a 

/•r/ + 7i^ — -a — 6 — 5 (vi — 3)<o, 

/V +(•/!— 3)2 <o 

et en multipliant par dr qui est néf^atif si r décroit de r^ à zéro, 
on a 

dr , rfyi 

r -7 Tvr>0 5 



c'est-à-dire 



/• 7) — 



/• 7i — 3 /•, rjo — 3 



Le second membre est une quantité finie ; le premier est une 
somme de deux termes, dont le premier seulement devient infini 
lorsque r tend vers o, Tinégalité est donc impossible. 

On a donc bien 



o<Yi <o 

Les limites peuvent-elles être atteintes? 
Si on a- ' == o, Téquation se réduit à 

/V+yi^+5vi = o, 
et admet la solution 

r^ =0. 



C. Q. F. D. 



LIMITES DE r, "/i 

Si - = — 3, Téquatlon se réduit a 

'V + ("^ — 3) (-^1 + a) = o 

qui «idmct la solution yj = 3. 

Le premier cas correspond h un point intérieur de la masse 
supposée homogène; le second correspond aux points extérieurs 
à la masse, et en particulier au cas où toute la masse est concen- 
trée au centre. 

Je dis que dans ces cas particuliers, les seules intégrales 
admissibles sont les solutions yj = o et y^ = 3. 

En effet : 

Supposons - = o. L'équation se réduit à 

/V+'^i^ + 5yi=o, 
dr dr^ 



r y;" + J7i 

dont la solution est • 



5 7, + ô 



Pour r = o, on aurait 



^i 



7.4-5 



00 ; 



donc 71 = — 5 ; or pour r quelconque, 7, est toujours positif : 
donc la constante d'intégration ne peut être que — oo , et l'inté- 
grale se réduit à 7, = o. Il en est donc de môme pour 7j = 3, 
dans le cas de la masse concentrée au centre. A l'extérieur de la 
masse, on a l'équation 

fW+'fi^ 71—6= o, 

dont l'intégrale générale est 

L/•5=L-i:f^+c^ 

7. + 2 

La constante d'intégration ne peut être que o et on a 'f^ = 3. 
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Ce sont les seuls cas où y, puisse atteindre ses limites. En effet, 
supposons d'abord que pour une valeur i\y on ait r, = o, pour /* 
différent de zéro, r, est supérieur à zéro, donc r/ [r^ = o, et en 
portant ces valeurs dans Téquation de Clairaut, il vient 

D^ = o, 
mais comme 

/•D'+3D = 3?, 

on a D = p. La densité en ce point est donc égale à la den- 
sité moyenne, et comme par hypothèse la densité ne peut 
diminuer vers le centre, on a un noyau homogène. 

Supposons de même que pour r = r^ on ait r^ = 3, de 
même que précédemment il faudra que 7/ = o ; donc 



D 



p =0 en ce point, et pour /• >ry on a encore p ==0; p ne peut 
augmenter puisqu'il y a équilibre, et par suite r, = 3 pour /• > i\. 

Il reste à voir ce qui se passe pour /• < t\. Je dis que Ton a 
encore r, = 3. S'il en était autrement, r^ serait plus petit que 3, 
r^ (7, — 3) serait une quantité qui diminuerait avec r, sa dérivée 
serait positive. On aurait donc 

7V+5(7, — 3)>o. 

Or, Téquatlon de Clairaut peut s'écrire 

,V + 3 (r. - 3) + (r.- 3)' + [^ + ô] ["i +yi] =o. 

Les deux derniers termes étant essentiellement positifs, il en 
résulte que l'on devra avoir 

'V+5(t,— 3)<o 

ce qui est contraire à ce que l'on a établi plus haut. 
Pour /• < /'y, on a donc encore r^ = 3, et par suite 



[^+a]=„. 



FORME DES COURBES ISTÉGRÀLES 



7*> 



/ 



Donc p = o pour r < r^. 

La masse de la planète est donc entièrement concentrée au 
centre. 

Forme des courbes intégrales, — Etudions maintenant la 
Torme des courbes représentées par l'équation 

/V+V + 5r. + ^' (i +r.) = o. 

Traçons les axes de coordonnées Or. et Or, 
11 suffit évidemment d'étudier les courbes comprises dans les 
régions où Ton a 

o<Y, <3. 

Traçons la droite D d'équation r, = 3. 

Par un point du plan passe en général une courbe et une seule, 
41 moins que r/ ne soit pas une fonction holomorphe. 
Or, on a 

V+5Y.+ — iy-(l+T.) 



v=- 



r 



r/ n'est pas holomorphe si r == o ; la droite r = o est une courbe 
intégrale de l'équation ; il y a sur cette droite deux points remar- 
quables r, = o et Tj = — 5, pour lesquels r/ est indéterminé 

/•D' 
à condition que -yr— tende en eDFet vers zéro avec /*. Or^ on a 



rW 
"IT 



VD 



Le second membre est essentiellement négatif, la densité d'une 
couche étant inférieure à la densité moyenne de la masse inté- 
rieure à cette couche, puisque par hypothèse il y a équilibre. 

Lorsque r tend vers zéro, p tend vers une limite p^ et l'on a 

?<l^<?o. 

? P 

?0 ^ 
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9 1 . 1 '^ 

mais I tend vers zéro avec /*, donc -p- 

Po ^ 

zéro et par suite — p-- tend vers zéro avec r. 



I tend aussi vers 



Aux deux points 



( T. == o 
/ =^0 ^ 



\ fi = 3 



les deux courbes intégrales ne sont pas définies d'une façon 
unique. 

Faisons maintenant r suffisamment grand, alors 



D 



= — 3. 



L'équation devient 



/V + Y,^ — r, — 6=0. 



Elle s'intègre immédiatement, et son intégrale générale est 



,.5 /.- 



3) 



7. + 2 



=:C•^ 



Les courbes représentatives sont toutes asymptotes à la droite 
7, = 3, mais elles ne sont pas nécessairement asymptotes a la droite 

/•= o, car on a — r— ^ — 3 pour 

/• = o : le raisonnement ne s'ap- 
plique plus. 

Considérons le rectangle OABC 
(fig. Il) formé par les droites 
Y. = o, /• = o, Y, = 3, r = i\, i\ 
étant une valeur quelconque. Par 
un point P, pris sur le côté OA, 
passe une courbe intégrale et une 
seule. Pour cette valeur 




Fig. II. 



r.=--5->o; 



la portion de la courbe qui pénètre dans le rectangle est donc 
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dirigée vers le côté AB, quel que soit le point P.. Une courbe 

intégrale ne peut donc sortir du rectangle par le côté OA, car si 

elle ressortait en un point P', elle aurait une forme analogue h 

la courbe PMP' (fig. 12), ce qui est impos- 

dr 
sible puisque --z — ne peut être nul dans le 

rectangle. 

Donc cette courbe ne peut sortir du rec- 
tangle que par un point M de la droite AB. 
De même, si l'on considère la courbe inté- pjg, ,2. 

grale passant par un point Q du segment 
GB, la tangente en ce point a comme coefficient angulaire 

dv r 




d'fi 



^b^^] 



quantité qui est négative puisque Ton a 

2/'D' 



D 



— 6. 



La branche de courbe qui pénètre dans le rectangle est donc 
dirigée vers l'arc AB, le même raisonnement que précédemment 
montre que cette courbe ne peut sortir du rectangle que par un 
point N de la droite AB. Evidemment, d'ailleurs, le point M est 
à gauche du point N puisque deux courbes intégrales ne peuvent 
pas se couper. 

Les points M ont donc une limite M' lorsque le point P se rap- 
proche de C, de même les points N ont une limite N' lorsque le 
point Q se rapproche de C. 

Le lieu des points M' sera une courbe complètement comprise 
entre les droites 7^ = 0, yj =3. Il en est de même pour les lieux 
des points N'. 

Est- il possible que M' et N' coïncident ? Supposons, ce qui 
sera démontré un peu plus loin que cela arrive. Alors il y 
aura une courbe intégrale et une seule qui passe par le point 0. 
Cette courbe sera évidemment la courbe intégrale convenable, 
puisqu'elle restera constamment comprise entre les droites '/\ = o, 



p « 
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Y, = 3. lien résulte que la fonction r, (/•) est parfaitement déter- 

mince, est donc parfaitement déterminé ; e est alors, en 

e 

désignant par E la base des logarithmes népériens 






La constante e^ sera connue si on se sert de Tune des équations 
précédentes, par exemple 



e/-^ 



'>=f ("+-^)- 



Une seule courbe est acceptable. — Démontrons maintenant 
qu'il y a une seule courbe intégrale acceptable. Nous allons, pour 
cela, démontrer qu'il y en a une et une seule qui passe par le point 
r-.=^o^ 7j = o, et que toutes les autres passent par le point r=:o, 
Tj = — 5 . 

Démontrons d*abord qu'il y a au moins une courbe passant 
par le point yj = o, et une passant par v) = — 5. 

Nous avons déjà vu quel était le signe de rr/ pour diverses 
valeurs de vj ; /V est négatif si r^ est inférieur a — 5, ou supé- 
rieur à + 3 et positif si t) est compris entre o et — 2. 

Dans les autres intervalles, les deux signes sont possibles. 

Pour Ti == — 5, 7*7]' peut être nul ou négatif, mais non positif. 

Pour 71 = — 2, frri' ne peut être que positif ou nul. 

/•D' 
Enfin, nous avons vu que — ^:- — est une fonction tendant vers o 

avec /', de façon, qu'étant donné un nombre positif a voisin de 
zéro, on peut déterminer un nombre j\ tel que l'inégalité 

/• < f\ 
entraîne 

o>-^>-a. 

L'équation de Clairaut peut s'écrire 

'W + (•/! — ?) (r. + y) = o, 



UNE SEULE COURBE EST ACCEPTABLE ^9 

^ et Y étant des nombres positifs, tendant Tun vers zéro, Taiitre 
vers 5 lorsque a tend vers zéro. 

Ceci va nous donner des indications sur le signe de rf\\ lorsque v 
est inférieur a r,. On a le tableau suivant: 



'I 



00 



y o p. -|-ao 



rr' 



I '^ i 



•j 



Il y a deux intervalles pour lequel le signe de /V est douteux^ 
mais ces intervalles sont moindres que ceux précédemment 
déterminés. 

Ceci posé, je vais démontrer que pour '•</',,*ir, est nécessaire- 
ment inférieur a p. Supposons que pour une valeur i\<ii\ on ait 
Tj = T,^> ^. Je dis que si r < i\ on aura encore r > p. En effet, 
rinégalité ne pourrait cesser que si 7/ = ^, mais alors dans Tin- 
tervalle, V^^^^ positif^ ce qui est contraire à la règle des signes. 

Or l'équation peut s'écrire : 

'V+(^-p)(^ + ï)+2(-^ + aj (i +ri) = o; 



par hypothèse 



/•D' 



D 



-> — «, 



on a donc 



'^/+(^.— P)(^.+r)<o 



A fortiori on a 



■7 — ôhT" H — < ^; 



et, si on intègre de r, h r<r^, 



L ' ' > L ' 1— 

L 1 7r<L 



r T.— fJ /•, Vil—? 
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Or 7, — ^ serait positif : on aurait a fortiori 

L_<L h 



Le second membre est fixe, le premier peut grandir au delà de 
toute limite, Thypothèse faite est donc absurde. Par suite 

T, a pour limite zéro lorsque ;• tend vers zéro. 
Je dis maintenant que l'équation de Clairaut admet au moins 
une intégrale tendant vers — 5. 

En effet, étant donnée une valeur t,^ comprise entre — y et — 5 

— 5<vii<— Y, 

considérons une courbe intégrale passant par le point (z*,,'^»,), je 
dis que si r^ < z*^, r,^ < v. Autrement r/ serait négatif pour Yj > — - y^ 
ce qui est contraire à la règle des signes. 
On a donc pour,/* suffisamment petit 

— 5<ri< — y; 

mais Y pourra être choisi aussi voisin de 5 que l'on voudra ; il en 
résulte quer, tend vers — 5. 

Soient t,^ et t,j les deux intégrales dont nous venons de démon- 
trer l'existence. Il leur correspond des intégrales e^ et e, de 
Téquation linéaire du second ordre 

46>"D + le'D -f -1 D' + A e'D = o. 
L'intégrale générale de cette équation sera 

A et B étant deux constantes. L'intégrale générale de l'équa- 
tion de Clairault sera 

re' kre^ -\-^re[ 



P 
Ti 



RELATION ENTRE V APLATISSEMENT, LA PESANTEUR, ETC. 8c 

Or, quand r tend vers zéro, 7\^ tend vers — 5, r^^ tend vers zéro. 
Il en résulte que Tj^ — Tj^ pourra être inférieur à — 4- On aura 
donc 

re/ re! 



i~^^<_4. 



c c 






e, e, r 



.-4 . 



c'est-à-dire que la limite de —^ est zéro quand r tend vers 

zéro. Donc quand r tend vers zéro, la limite der, est la môme que 
celle de tj^ tant que B est différent de zéro. Si B = o la limite 
est fjQ. Il en résulte qu'il y a une intégrale et une seule de Téquà- 
lion de Claîraut tendant vers zéro avec /• : c'est celle qui corres- 
pond a B = o. 

Relation entre l'aplatissement, la pesanteur et la force 
centrifuge à l'équateur. — A l'extérieur de la planète, comme 
cela a déjà été remarqué, on ne peut parler de surfaces d'égale 
densité, mais on peut encore parler de surfaces de niveau. 

Pour ces surfaces, on a toujours l'équation 

AV=AVo + A(IIY)=o; 

maïs AVq = — 4*^^^ est nul à l'extérieur de la masse fluide. On a 
donc 

A (II Y) ==o 
A(HY) =y| H'' +4 11^ — "-^^^0^ Hl= o. 

Faisons /i= 2, on a 

H''+— H'— -4-11 = 0. 

L'intégrale générale de cette équation est 

II = A/-'- 4- B/-3 ; 

PoiNGARÉ. — Figures d'équilibre. 
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mais à l'infini II doit être nul, il faut donc qae A == o. On a donc 

H = B/-3 ; 

si M représente la masse fluide, V le volume enfermé par la sur- 
face de niveau considérée, on a 

" V ~ /-^ ' 

en posant L = — - — . 
L'équation 



ilevient 






D'où 



On a donc 



9 r-^B 3.a)V -| 

471 L C/- ^ 6C J 



re -\-2.e ^^ ' 



i5 

SttC 

i5w* 



2 3 



87iD * 



Comme on a posé 



©n a finalement 



re'=er^ 
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Cette relation établie pour Textérieur de la masse, est encore 
valable a la surface. On peut lui donner une autre forme. 

La force centrifuge a Téquateur est égale à wV ; la pesanteur 

à la surface est -5- tt/'D. Si on désigne par o, le rapport de la 

première à la seconde de ces quantités, on a 

e(ri+2) = -^3. 

Or r, est compris entre o et 3 ; il en résulte la double inégalité 

^ >e>-^; 



4^2 

Tj = 3 seulement si la planète est homogène. L'aplatissement 
ne peut atteindre la limite -^que si la planète est homogène. 

Détermination complète des coefficients. — Nous allons 
maintenant démontrer que les équations 

AHY^-P'+/P" IIY 

ru 

qui déterminent les coefficients II du développement 



v=v„+2hy, 



n'admettent que la solution II = o. 

En effet, soit n Tordre de la fonction sphérique, posons 

L'équation devient 

A(Z/'»D^ ) == -^ — Zr»Y ; 

/•"Y est un polynôme sphérique ; en appliquant la formule 



(PU)=ru'' + 2^idli. U' "Ip, 



# w 
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on ol)tient finalement 

(ZD)" + ,^!1±1 (ZD)' = ZD''+ i^ , 
ce qui, développé, s'écrit : 

Z"D + 2Z' (d'+ ^^-y^A + -^^^ ZD' = o. 

Si l'on fait le même changement de variables 

r/J 

on a 

;•£ +e'* + (2W+ iJîH jT — (s + /i — i) = o. 

Celte équation se discute de la même manière que l'équation de 
Clairaut ; — — — est compris entre les limites o et — 6. 
Ou est donc amené à considérer les deux équations : 

/•£'+ £^ +(2/1+ l) £=0, 

/•£' -f- ^^-\- (2/1 — 5) £ — 6 (n — i) =0. 

On construit les courbes en portant £ en abscisses et te' en 
ordonnées (fig. i3) ; on obtient deux paraboles, et comme pré- 
cédemment on démontre que £ ne peut être compris qu'entre 
o et 3. 

A l'extérieur de la planète on a A (IIY) =0, ce qui peut s'écrire 






L'intégrale générale de cette équation est 
Il devant s'annuler à l'infini, on a a = o, et par suite 



DÉTERHUNÀTION PRÉCISE DE L'APLATISSEMEKT 
MaisE doit être compris entre o et 3 ; doac b doit être nul. 




Fig. i3. 

Il^o est donc la seule solution possible, à l'extérieur de la 
masse. Il en est donc de même à l'intérieur. 

Ce résultat ne subsiste pas si n ^ i ou 3, mais nous avons 
Jéjhe 



Détermin&tiOD précise de raplatlsaement. — Nous avons 
trouvé 

iV=iV,+4i.p'C, 

i(HY)=4'=?'i:, 



3Ç 



Nous avons donc l'équation 



A(HY)=^™fVY, 
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et, en se rappelant la formule 

A(IIY)=Y[n"+^-<^H], 

dans laquelle /z = a, on arrive pour déterminer H à l'équation 
linéaire avec second membre : 

, 2H' 6H 4 , 

L'intégrale générale de l'équation sans second membre 
étant 

où a et [3 sont des constantes d'intégration, on déterminera l'in- 
tégrale de l'équation avec second membre, en regardant a et jï 
comme des fonctions de /*. La méthode de variation des cons- 
tantes donne 

(i) aV2+^V-« = o, 

par suite 

et, portant ces valeurs dans l'équation primitive, on a 
(2) 2a'/- — 3^3V-^ = 4" ^^?''-. 

Les équations (i) et (2) déterminent a' et ^' 

«=-75-^?. 



3' — — i?i P^'r^ 



On a donc 



II) 1 * i5 f ' 



//• 
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Cette fonction comprend deux constantes arbitraires a et b. 
Comment faut-il les choisir ? Soit r, le rayon de la planète ; pour r 
supérieur à r^, p = o, f ' = o, par suite les intégrales se réduisent 
à des constantes 




et 




Pour ;• =00, II doit s'annuler. Donc le coefficient de /• est 
nul; on a donc a=f\. 

De même pour r = O5H doit être fini; donc le coefficient de 
r~^ doit être nul, si r=o, on a donc b = 0. 

II est ainsi complètement déterminé. 

Portons cette valeur de H dans l'équation 



•2D=-#- [11+ 



<ùh'^ 



4'^ \_ 6 

Il vient 

ep'dr p-rrr f ^P 

^ 5/'D f ' 



8 7 1^*^ 

/•^ dr + 



8tt D 





Que devient cette expression à la surface pour r = ;j?D devient 
1), la densité moyenne de la planète, et on a 



3 / , - , 3 to^ 



(3) 




îp étant comme on l'a déjà indiqué le rapport de la force cen- 
trifuge à l'équateur, à la pesanteur. 



9 • 
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Moments d'inertie de Vellipsoïde, — Calculons l'intégrale 
/ Y^rfcr, Y étant la fonction 



r' 



On a 

X = r\/i — [x^ cos <p, 

y = r\/i—it} sin 'f , 

On a donc à calculer 

I I — 3a^ I rf[xrfy = 27r I I I — 6[Ji^+9[x* Irfjji 

—1 




/v. 



T= — zr— 



Calculons les moments d'inertie A, B, C de la masse 



On a 



Or on a 



AV = — 4^p, 

^-B = C^Y{z'-\-xyx, 

27:(A + B + C) = rAV(^*+y»+z^)rfT. 
AV = AV„4-A(HY) 
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et de plus 

rA(HY)/-VT = o ; 

en effet, on a d-z^f^drdyA'^-=drdv. 
Il vient 



f\W{i^dx = fii{mydrdi 



= I ^y r'dr.Ydv^ I f{r)dr j Y<f<7, 



ce qui est nul, car le second facteur est nul. On a donc. 

27r(A + B 4-C) = |'^Y„l■^d■z 

= fAYydrdv, 

et comme AV, ne dépend que de r, on a 

27t(AH- B + C) = 47t f^Yydr. 
Or 

inuis on a 

dy„ 4toD 

dr 3 ' 

donc 

A4-B + C=-Ç. / ,»J-{Ti,^)dr. 







Si on intègre par parties et que Ton remarque que A = B, on 
voit que l'on a 
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En négligeant des infiniment petits on a A = C ; et par suite 



=i[,,o, 




Calculons maintenant C — A : on a 

47r(C — A) = rAV(.s2— x2)rfT, 

47,(C— B) =f\Y{z^—f)d'. 

Donc 

8t:(C — A) ==—j\Yr'Y(k 



en désignant par Y la fonction sphérique 
On a, comme on a vu précédemment : 



r' 



Stz[c—a)=— fiY.r'Yd^ — f^myxck, 

La première intégrale est nulle, car c'est le produit de deux 
intégrales : 



C^^\r'dr j 



Xdi, 



dont la seconde est nulle. L'intéorrale 

D 



/ A(IIY);nWT= / M^.r^\^d-z 




^ffi-rWr I Wcr; 



et comme on a 



(4) / W=r = 
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il vient 
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8<C-A)=-~^ I ^ 




AIIY) ,, 



Tzeo'r^dr 



i 



ep'r^dr. 



Donc on a 



A = Z- I e^'r^'dr 

ID 




i 



quantité positive. 
On a donc en se reportant à la formule (3) page 87 

. ? __ 9 (C-A) 
' 2 ~ 87: /-JD * 

On aurait pu faire une hypothèse différente de celle de Clal- 
raut, par exemple, supposer un noyau solide non homogène, 
recouvert d'une masse fluide; la relation établie plus haut sub- 
sisterait encore. 



Comparaison de la théorie et des observations. — Posons 

C — A 

— — = J, on a 



* 2 




m # 
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e^ est mesuré par la géodésie, on trouve — :r-p ; y est mesuré 
par la physique, on trouve ■ • ; J est mesuré par la préces- 
sion des équinoxes, on trouve o,oo32j53 = . > ^ — on a donc 

00<J • O I 

l'égalité 

r M3 = ' ^=0,49145. 



1 '"1 



Revenons maintenant à l'équation de Clairaut 

/-/i'D + (ri^+ 5ri)D + 2/-D'(i -\-r) = o ; 
on a 

^'' L "^ J 2v/i+'/i 'fVi+^i V'+^i 



__ 5/-^D r 



i+'^iH ^ + 



/V , r(i+r,) D' 



10 ô 

On peut donc écrire Téquation de Clairaut 

^' L "^ J v^i+T, L 10 2 J 



-1 



J'intègre de o à /'^ et j'ai 




\/i+^ 



L 2 10 J 



On peut itppliqucr le théorème de la moyenne, puisque 5/*^ D 
est une fonction positive ; \ désignant une valeur comprise 
entre o et 3, et correspondant a une valeur de /* comprise entre 
o et i\^ j'ai 
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Voyons les variations de 



ifL=K; 



v/' + i 

lorsque Ç varie de o à 3 on trouve : 

l = o K = I 



% = -^ K = I , oooyS maximum 



2 



K = I , 0002 



Ç = o,5[5 K -= I 

; — V,, = 0,544 K, = 0,99954 

$ = I K = 0,989 

; = 2 K = 0,9 

^ = 3 K = o,8. 

On voit que K varie très peu et reste compris entre i,ooo^5, 
et o, 99954. K est donc très voisin de i ; et si dans la formule on 
remplace yj, par sa valeur 0,544? on trouve 



Cr'Ddr 



='TÂz\/^ + ''h 



/•,^'D, 5K 

= 0,497 — > 0,49663. 

Précédemment on a trouvé pour la môme valeur 

0,49145. 

Comment peut-on expliquer cette différence? 
On peut admettre qu'il y a des erreurs sur les quantités 
qui entrent dans le calcul : c'est-à-dire des erreurs sur J et e,. 

La première équation donnerait alors comme valeur exacte : 



'/'•' 



Ddr 

= 0,49143 — 3o5,3i oe^ + 155,26 oj. 
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La seconde équation donnerait : 



d. kJ?' 



— {ri^-+-2)de 



2K fr'Ddr , ^ . 

Supposons 8J = o ; pour que ces formules soient comparables, 
il faudrait que Se^ soit égal à — 0,000027 et la valeur de e^ serait 

alors TT-TTTT. 

295,05 
Si on suppose que l'erreur est commise sur J, on a 

SJ = 0,000034. 

J'aurai alors pour valeur -r tt . 

* 002,18 

Cela est-il admissible ? Remarquons d'abord que le J qui 
figure dans cette formule n'est qu'approximativement — ~ — . 
Cela revient à remplacer J par 

J 



2 , 



2 
L'erreur commise est 6J =-:j- J^ : on trouve 

SJ = 0,000006. 

Ce n'est pas de là que provient l'erreur. 

D'autre part les mesures de nutation sont assez précises. Si 
on désigne par N la constante de la nutation, on a 

N = — ^J [5,36542]; 
en désignant par [jl le rapport 7^ de la masse lunaire, à la masse 
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terrestre ; le nombre entre parenthèse étant non le coefficient 
mais le logarithme de ce nombre. 

Soit p la précession des équinoxes : c'est la somme de deux 
termes 

celui du à l'attraction de la lune 

,y = 866i35 — '^ J = 34,38 ; 

I +[X. 

celui du à l'attraction du soleil 

//' = 487 1 ,o5 J = 1 5,95. 
Posons 



on aura 



op 





1^ 




w • 




I 


+ !^ 




SX 

N ■ 


8î 




oJ 
J ' 




4 s 


oJ 
+ J 




+P"- 


oJ 

J 




+ P- 


SJ 
J 


'7 





p est mieux connu que N, on peut donc supposer 8y; = o ; et 

oJ 
résolvant les deux équations où Ton regarde— r- comme égal 

o,oooo34 x3o5,3i =o,oio3, on trouve : 



0£ 

e 



0,OIÎ> 



oiN ^ 

— — — = o,oo5. 

n 

Cela revient h prendre pour N la valeur 9"i7 au lieu dc()''2i 

et pour UL la valeur -77 — - au lieu de 77-. Ces chiffres paraissent 
* ^ 8i,j 0:1 

difficiles à admettre. 

On peut alors supposer que l'erreur provient de la somme des 

erreurs commises sur chacun des termes. 
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Une autre explication serait que Thypothèse de Clairaut n'est 
approchée qu'aux termes du second ordre près. M. Callandreau, 
qui a fait le calcul, a montré que considérer les termes du second 
ordre reviendrait à creuser rellipsoïde de g mètres seulement. 
C'est négligeable. 

On pourrait alors se demander jusqu'à quel point l'hypothèse 
de Clairaut serait exacte : par exemple, on peut supposer la terre 
complètement solide sauf la mer. Au moment de la solidifica- 
tion, les couches obéissaient a la loi de Clairaut, mais depuis ce 
moment, la vitesse de rotation aura pu varier (par exemple, les 
marées diminuent cette vitesse de rotation) ; et de plus l'aplatis- 
sement des ellipsoïdes internes aura pu se modifier par suite des 
tassements qui ont produit les continents et les montagnes. 

Enfin on pourrait supposer qu'une portion intérieure de la 
terre est restée fluide : les deux parties fluide et solide, agiraient 
alors séparément dans le phénomène de la précession : mais 
Thypothèse de la fluidité intérieure de la terre n'est guère vrai- 
semblable. 

Jusqu'à présent nous avons considéré les planètes comme 
fluides et nous avons montré quelle devait être la distribution 
des densités à l'intérieur. 

Supposons maintenant que l'on ait une masse solide, de forme 
quelconque, où la distribution des densités est quelconque, maïs 
entièrement recouverte d'une masse fluide de faible épaisseur, 
soumise à la pesanteur et à la force centrifuge due à la rotation. 

La surface extérieure sera une surface de niveau, sur laquelle 
on aura 

mais à l'intérieur les surfaces équipotentielles ne sont pas néces- 
sairement les surfaces d'égale densité. 



CHAPITRE V 
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Théorème de Stokes. — Supposons que Von ait déterminé la 
forme de la surface extérieure S, la masse M et la rotation o), 
on pourra déterminer en tout point extérieur la valeur du poten- 
tiel V. 

Il y a une Infinité de lois de distribution de lu matière dans le 
noyau, qui peuvent donner une même surface extérieure, mais 
quelle que soit cette distribution, le potentiel V, à Textérieur, est 
complètement déterminé par les éléments S, M et o). 

En effet la fonction V satisfait à l'extérieur de la masse à l'équa- 
tion AV = o. A l'infini V = o ; sur la surface on a 



^ (x' +f) = U = K, 



et enfin l'égalité 



d'j = — 4^M. 



dn 



Le problème de la détermination de V est complètement déter- 
miné. En effet, soit V une autre solution. Sur la surface, on a 

V — V'=U — F = K — K'=constante. 

On a d'ailleurs A(V — V') = o. 

PoiNCARÉ. — Figures d'équilibre. 7 
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Donc V — V est le potentiel d'une charge électrique répartie 
sur la surface, puisque c'est une fonction constante sur cette 
surface. 

La charge a pour expression 

an 

l'intégrale étant étendue à la surface extérieure ; mais elle est 
nulle puisque Ton a par hypothèse 



4t:M, 







V — V'.le potentiel d'une charge nulle, est donc nul. C'est ce que 
Ton voulait démontrer. 

L'expérience du pendule qui peut faire connaître les valeurs 
de ^ à diverses altitudes ne peut donc rien nous apprendre sur la 
distribution des masses intérieures, que ne nous fassent connaître 
les éléments considérés dans le théorème de Stokes. De même, 
nous ne pouvons rien apprendre de plus par les perturbations 
produites sur le mouvement des autres astres. 

Développement de ^^, en série de polynômes. — Soient 

O l'origine des coordonnées (fig. i4), M' un point attirant de 
coordonnées x\ y', z', de densité p', M [x, y, z) un point attiré. 
En prenant les coordonnées polaires on a : 

X = /-y/i— [j? cos cp, 
y =^ /y i -^- ji.^ sin cp, 

z = /'[A., 




] i/=r'\/i — ^ Y-'^ sin ç', 



Fig. 14. \ -s J^7 [A j 



V. 



SUR UN CERTAIN DÉVELOPPEMENT EN SÉRIE 99 

MM'* = r* — 2 ri^ cos y -+- r'- 



= x^ + y^ + z^ — i {x^' + y y' + zz') + x" + y" 

I I 



«/2 



MM' 


/ 


■vA 


2 /•' cos y 


+ 






;• 


/•^ 




/•* 


— r' 


' + y' + =S 








/•'■^ = 


-y* 


+ y'^+^'^ 






os Y 


a: 


jj -f- ,. ,y -t- 


r 


/•' 



Développons ».-.., suivant les puissances croissantes de 



MM' ^ ^.....-.....v.. ^ , 

le coefficient de — — sera un polynôme II,, [x, y, z ; x', y', ^') . 
On doit avoir les deux relations 

^Vmm^ 

la première étant obtenue en considérant .r, //, : comme varia- 
bles, x'y'z', comme fixes; c'est l'inverse dans la seconde. 



■^ ( MM' j ~" 
^' \ MÂT/ ^ 



d^ 


I 


MM' 


fl^ 


I 


\J 


MM^ 



+ 


ô- 


I 

Mi\l' 




d,f 


. 1 . 


ô- 


I 
MM' 



ô' ' 



MM' 

- = o, 



d- 



I 



MM' 



On sait que l'on peut écrire 



I „ /•' ...//•' ^ " 



v/ 



= I+P,— +... + 1^ - +• • 



1—2 .-^cosyH r 



où P„ est un polynôme en cos v de la forme 

AqCos'*y + '''^2^^^"~^T H~- 
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si on y remplace cos y par sa valeur, on aura 

Le polynôme Q sera un polynôme homogène et de degré n 
par rapport h chacun des systèmes de variables [x^ y, z^ /•), 
{x\ y' f z'y r'). Il ne contiendra d'ailleurs /• et /•' qu'à des puissances 
paires, on pourra donc remplacer /*^ et /•'^ par leurs valeurs 

^2 _J_ y2 _|_ ^2 gt ^/2 _!_ y/2 _^ ^/2 . 

on aura donc : 

Q!n{^,y,z\x',y\z') 



cosy+ 



;./2 Zà 












MM 

Q' représente un polynôme homogène et de degré n par rapport 
à chacun des systèmes dje variables x^y^z\ x' y y' ^ z' , chacun des 
deux ensembles y entrant de la même façon que l'autre. 
On a d'ailleurs 



Donc on a 



A'- 


I 

MM' 





A' 


Q» . 


— o; 

4 



et, comme /• ne dépend pas de x\ y\ z', on a : 

A'(QQ=o. ■ 

Q' est donc un polynôme sphérique d'ordre n par rapport aux 
variables x',y', z' et par suite par rapport aux variables {x, y y z). 
Nous savons qu'il y a ip + i polynômes sphériques d'ordre p^ 
de forme 

p — ,.p qH'A I M ^1 2 ^^ ^^ r ^ 
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OÙ q prend les valeurs 

— P, — T'+i, —1.0,1, 2, . . . p — i,p. 

On peut donc écrire 

formule où P est un polynôme fondamental en .r, y^ z, Q un 
polynôme en {x\ y, z') ; Q' est d'ailleurs fonction de 



cos 



Y = UL[X' -f" V/ï 1-^' V/ï y-'^ cos (cp '^j 



Q' ne dépend donc que de la diDFérence {o — 'f') : il en résulte 
que le coefficient de P^^ne peut être que 

O =r'^e— «^?/i ul'^^ =î ^^ — ^— 



et on pourra écrire : 



l ^^ ^p, 9 1 /). <? 1 /^ - ?• • 



MM' 

A^^ est un coefficient numérique que Ton peut déterminer faci- 
lement. 

Application, — Supposons que /• soit suffisamment grand par 
rapport au rayon de la planète. On pourra développer le poten- 
tiel suivant les puissances croissantes de — — et on aura 

D'après ce que Ton a démontré plus haut, la connaissance de 
S,Met(o suffit pour déterminer V, donc ces éléments permet- 
tent de déterminer B^, ^ ; or, on a 

V— / .P^ — VA^ / o'F d- 
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il vient donc 



Bk=Ap, / p'PV-,'^^ 



P'^ ç est l'un des polynômes sphériques fondamentaux d'ordre p. 
On doit donc savoir calculer 



P'p;,rfv. 



Considérons d'abord le polynôme d'ordre o : Po,o= i ; on a 



p'rfT' = M. 



Les polynômes d'ordre i sont x' ^ y\ z' : on connaît donc la valeur 
des intégrales 

pVrfr, 1 ç'y'd'z, j p'z'd'z. 



Ce sont, à un facteur près, les coordonnées du centre de gra- 
vité de la masse. Supposons que la surface extérieure possède 
un centre de symétrie, et prenons l'origine des coordonnées en 
ce centre ; les intégrales sont nulles puisque ce sont les coeffi- 
cients des termes en —3-, -4", "x? dans le développement du 

potentiel et que les surfaces équipotentielles doivent admettre 
l'origine comme centre de symétrie. Le centre de gravité de la 
planète coïncide avec le centre de symétrie. 
Les polynômes de second ordre sont 

^Vi y^^\ ^^, 
z'^ — jc'\ z'^ — y'\ 

Si / x' y' ^'d'z est nulle, ainsi que les deux autres intégrales 
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analogues, les axes de coordonnées sont les axes principaux 
d'inertie du système. Si / (x'* — y'^) prfr est connu, on connaî- 
tra la différence entre les moments principaux d'inertie. 

On a trouvé une relation entre e^/'j, D^ et o); les deux premiers 
éléments sont donnés par la connaissance de S ; D^ est connu 
puisque Ton connaît M et S ; cp est connu puisque Ton connaît o) ; 
C — A est aussi connu. On doit en conclure que la relation déduite 
de l'équation de Clairaut est vraie, quand même la distribution 
intérieure de la matière ne serait pas conforme a l'hypothèse 
faite pour l'établir. 

Relation entre l'aplatissement, la pesanteur et la force 
centrifuge à l'èquateur. — Proposons-nous de calculer les pre- 
miers termes du développement de ^,,^ • 

On a d'après la formule du binôme 

I I r i ( iî^ i ^^ \ 



■^TK—'"'''^-—) + • • -J' 



négligeant les termes de l'ordre de — j-; nous avons 

I I , /•'cosY , i r /'^ . 3 2 2 1 



MM 

3 



2/ 



J- (xx' + i/y' 4- sz'Y 



= — + -p- {xx' -\-i/y'+zz') 
j-[(^' + y'+z') {x'''+y"+z'^) — i{xx'-[-yy' + zz'y]. 



2/ 

Le coefficient du terme — r- peut s'écrire 

27* ^ 

— 6{xx' yy' -}-xx' zz' -\- 1/y' zz'), 
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On a alors 



y . P'^T 






dx 



MM' 



- y") ?Wt 



I (.r^ +y^—2z') , ^^, _^ ^„ _ ^_„^ ^,j^ 



4 



M» 



-«Sf 



a'y'^'ch'. 



Si M est la masse, a, ^, y les coordonnées du centre de gravité, 
A, B, C les axes principaux d'inertie, on a 

V= -^ + ^^f+^l+Il +_L^[3(.^_,^) (B- A) 

Supposons que la surface soit de révolution, que Torigine des 
coordonnées soit placée au centre, que les axes de coordonnées 
soient les axes principaux d'inertie : alors 

B = A 

et 

Calculons maintenant la valeur que Ton déduit pour V de la 
forme extérieure de la surface, que nous supposons un ellipsoïde 
d'aplatissement e^. Ce potentiel est de la forme 

M 
V=-^+IIY + SQZ, 
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OÙ Y est la fonction 



,^^_±t 



HZ 



r' 



et Z est une fonction sphérique différente de Y. 

H et Q doivent être d'une forme particulière, puisque le poten- 
tiel doit s'annuler à l'infini, nous avons déjà vu que l'on avait 

H= "" 



r' ' 



_ 9. 



Q= 



jjin + 1 



Ilg et q^ sont des constantes ; on a donc 



M_ H.Y .VI fi.:/^ 



r /•' 






A la surface libre on a 



U = V + — (^^ + y^)= C'«. 

o 



Il vient donc 



M (0^-^ w^-^Y 

\ ^, 1 r ^ 

r o 

Posons 



HqY I V QqZ 



On a vu (p. 69) que l'on avait 



«1 



de plus on a 



/•.Y ' 



'• = '-^[' + ^] 
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Tous les termes étant petits, sauf le premier, on peut y négli- 
ger e\ et on a 

M^^ M M / e,Y 

r 






L'équation U = C® se réduit donc à 
M Me.Y , a)V? , co^/jY , H, Y , Sy„Z _^,, 



,.3 • ,^ + 1 



Comme par hypothèse nous avons un ellipsoïde : ^r^ = o ; donc 
V se réduit à 



M H,Y 



/• /•* 



et on conclut, en comparant les deux valeurs de V, que 

C — A 



Ho== 



On a donc l'équation 

M (Ue, C-AX toVî , toV?Y _^,, 

""VITT ^T^r"* 3 ' 6 —^ ' 

on en conclut que le coefficient de Y est nul, et que l'on a 

.2 ,.2 



Me, C — A (0 Vi 



On a d'ailleurs 



ir^ 2/'î 6 

_ M 

M 



f\ 



{j}^f\ = CO £r^ = cp 



1 T i? T «2 

' i 



Il vient donc 



2rJ ~^^ 3 6 ' 

C— A _ 2r,y, / (p\ 

rî - 3 r~T;' 
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et, en supposant le rayon de la terre égal à l'unité, on a 

C_A=if^(..-X). 

Pesanteur en un point de la surface. — Les composantes 

, , ôU ôU ôU 

de la pesanteur sont -r; — ,-t:— , -r— , et on a 
* ox oy dz 

La composante suivant le rayon OM est, s étant le cosinus de 
l'angle de la verticale avec le rayon OM, 



ÔU 



:j^=g-COS£ ; 



comme e est petit, g cos e est sensiblement égal a g, aux quan- 
tités du second ordre près, on a donc 

M]_ 

^~ dr ' 

.2,. 



_ M (C — A) Y 2ioh^ (o^rY 



Si maintenant on fait /• = i -]- 2^, 

I 



,.2 = ï — <^ 



on trouve 



g^=M — 2M!; + A(C — A)Y 
_£ ,_AL ^ 



3 6 



r 2? ^ I 



2W 


(o^Y, 


3 


V 2(0^ 


o>=Y 


^ 3 


3 


Y 1 =^?A'o 


'f?oY 


^ ' 3 


3 


Y+ ^'-^ 
^ 3 




5?\ 





io8 
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C'est la formule de Claîraut qui donne la pesanteur en un 
point, en fonction de la latitude de ce point. 

Influence de l'altitude. — Si on s'élève en ballon, la cor- 
rection est facile à faire, on n'a, dans la formule, qu'a rem- 
placer /• par I + Ç + A, au lieu de le remplacer par i + Ç, mais 
ce n'est pas dans ce cas que l'on a à faire la correction d'alti- 
tude : le plus fréquemment cette correction se présente pour un 

lieu situé sur un plateau, et il faut 
tenir compte de l'attraction duc à 
ce plateau. 

Il faudra donc tenir compte des 
aspérités positives ou négatives; il 
faudra d'abord choisir une surface 
définissant le niveau général du 
plateau : cette surface sera la sphère 
décrite sur OM comme diamètre 
(fig. i5), O étant le centre du sphé- 
roïde, M le lieu de l'observation. 
On calcule l'attraction due au vo- 
lume compris entre les deux sphères, 
l'erreur est peu considérable et d'autant moins que l'attraction 
des montagnes est très faible. On prend, en général, comme 
densité de la masse la densité moyenne des roches h la surface, 

c'est-à-dire D^, en représentant par Dj la densité moyenne 

du sphéroïde. 

Correction de Bouguer. — Supposons une couche sphérique 
homogène d'épaisseur h et de densité p ; le volume de cette 
couche sera 

47r/'i/i, 
et sa masse 

La correction apportée à la mesure de g sera 4'*^phrl et on a 



Mgr. I.). 



■6 ^^''■' 



D,,- 
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D'un point M situé à une certaine distance au-dessus de la 
couche, menons le cône tangent à la sphère (fig. 16) ; le cercle 
de contact partage la couche en deux portions, et nous savons 
que les attractions de chacune de 
ces deux portions sur le point M 
sont égales. 

Appliquons ici ce théorème. 
Nous avons supposé que la surface 
de la planète était limitée à la 
sphère décrite sur OM comme 
diamètre ; si donc on suppose 
que l'attraction de la masse MRS 
se compose de l'attraction d'un 
certain nombre de couches homo- 
gènes et minces, la partie de ces 
couches limitée à la sphère (OM) 
aura sur M la même action que la 
moitié de la couche totale. Il faut donc calculer l'attraction sur M 
d'une couche entourant le sphéroïde de densité p et d'épaisseur h 
égale a l'altitude du point M. La correction à apporter a la 
pesanteur sera la moitié de cette attraction. 

L'épaisseur h de la couche peut se représenter par un dévelop 
pement de la forme 




Fig. iG. 



le potentiel dû à cette couche, d'après une formule connue sera 

4tîAY 



2 



(2/2+ i)/ 



.n + 1 



l'accroissement de g sera 



6g = 



5SV 



2 



AY47r 



71 + I 
2/2 + I ' 



car on fait ensuite r = i dans la formule, et on a 






3L+yjl+LAY. 
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Si on se borne au terme n = o 



i: D ' 



si on en prend la moitié, c'est la correction de Bouguer. 

La correction totale comprend donc deux termes : i® celui 
qui provient de l'augmentation d'altitude sans tenir compte du 
plateau; 2° celui où l'on tient compte de la masse du plateau, on 
a finalement la formule 

8^ 2/1 3 ph 



go '\ 2 ^i^ 



l' 1 



Si on applique cette formule on trouve une valeur trop faible : 
la valeur fournie par l'observation est, en général, sur les conti- 
nents, inférieure à la valeur calculée ; elle lui est au contraire 
supérieure sur les océans et dans les îles isolées. 

M. Faye a ainsi trouvé qu'il vaut mieux ne pas tenir compte 
de la correction de Bouguer, ce qui correspond h des masses qui 
devraient se trouver sous les océans, et h des vides qui existe- 
raient dans les continents. 

A la surface des mers on a U = Uq. La surface U = Uy est une 
surface qui, se prolongeant sous les continents, doit donner une 
surface de révolution, mais il n'en est pas ainsi. Cette surface 
s'appelle le géoïde. 

Pesanteur en un point de la surface du géoïde. — Considé- 
rons une sphère de rayon voisin du rayon de la terre ; soit Ç l'al- 
titude du continent au-dessus de cette sphère, on peut écrire 

Ç' la distance à l'ellipsoïde, peut s'écrire 

Ç — Ç' est la distance de la sphère a l'ellipsoïde au lieu consi- 
déré ; 
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Soit h l'altitude du point au-dessus du géoïde, 



=E 



ks: 



Zf — h est la surélévation du géoïde au-dessus de l'ellipsoïde : 
on a 



K'-h==^(ll-rj-k)Y, 



Ç — 2^' se réduit a deux termes : celui qui correspond a Y = i , 
et celui qui correspond a Y = -^ * 



A la surface du géoïde U = ^^ = constante ; en un point d'al- 
titude A, on a 

OU 

dr 



U rr _J fi . 



On a d'ailleurs 



dU 






Donc 



U 



V 



_=,+V(II_A.)Y 
- = i+V(lI-A— y)Y. 

^^ 



Les coefficients Y sont connus : ils correspondent au développe- 
ment de 



0)2 






i^' + f); 



si on rétablit l'homogénéité, on a 



V _ ï S(H — A' — y)Y 



a- /• 



et comme on a : 

OV 

O' ^^ 
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on en conclut : 

I , s(«+i)(H-A-y)Y 






Les opérations géodésiques font connaître A, les opérations 
de nivellement font connaître îj, les observations pendulaires 
font connaître^. Théoriquement, deux séries d'observations suf- 
firaient pour connaître 2J, A etg^; mais pratiquement les observa- 
tions d^aucune des trois espèces ne sont assez nombreuses, ni 
assez précises pour qu'on puisse se passer d'une série d'observa- 
tions. 

Jusqu'à présent l'on peut dire que l'on a trouvé la surface du 
géoïde peu différente de celle d'un ellipsoïde de révolution. 

Nous allons, dans ce qui suit, montrer qu'il peut exister des 
figures d'équilibre autres que l'ellipsoïde de révolution. 



CHAPITRE VI 



FONCTIONS DE LAME 



Coordonnées elliptiques, — Considérons les quadriqucs homo- 
focales 

^2 t/' z^ 

+ --1 J- — i=o. 



2 ,.2 



Par un point de l'espace passent trois de ces surfaces : en effet, 
.r, 7/, z étant donnés, on a pour déterminer ")} une équation du 
troisième degré dont les trois racines sont séparées par les 
nombres a^y b^, c^. 

Il en résulte que l'une des racines, la plus grande, supérieure 
à a', correspond à un ellipsoïde; la racine moyenne, comprise 
entre a^ eti", définit un hyperboloïde à une nappe; la troisième, 
comprise entre b^ et c^, définit un hyperboloïde à deux nappes. 

Soit p^, [A^, v' ces racines, on a donc 

p2 > a^ > jjL^ > i- > v^ > c'\ 

Réciproquement si p, [ji, v sont donnés, on a trois surfaces se 
coupant en huit points, placés symétriquement par rapport aux 
divers plans de coordonnées. 

Si on ne considère que les points situés dans un trièdre trirec- 
tangîe donné j les trois nombres p, ul, v définissent un point et un 
seul. 

On peut écrire 



^2 j^2 ^2 
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Si on multiplie par X" — a^, et qu'on fasse ensuite \^ = a 



'/ 



on 'A 



~" (a^ — 0^^ i^a^ — c') 



on obtiendra donc 






et de même pour y et r, 



- ./ {f-c'W-c'){^'-c') 
" "~ V (c^ — o*) {c' — i^") 

Le signe qu'il faut prendre devant le radical correspond au 
trièdre des coordonnées que l'on considère. 
Soit f{x) un polynôme, et posons 

R=/'(?^). 

M==/-(.^^), 

RMN sera un polynôme symétrique des trois quantités p^, [jl^, v'^ 
racines de Téquation 



x'^ y^ 



et par suite une fonction entière de x^^ y^, z^. 
Si on pose 

ou bien 
OU encore 



ÉLÉMENT LINÉAIRE EN COORDONNÉES ELLIPTIQUES ii5 

et que Ton prenne pour M et N les mêmes fonctions de [x et 
de V, RMN sera un polynôme en x^^ y^-i^^j multîpliépara', parj:^ 
ou par xyz. 

Supposons que R soit un polynôme entier en p* ; on peut le 
décomposer en facteurs du premier degré et écrire 

R=n(p*-Xî). 

Alors 

RMN = n (p^ — ).î) ([X* — \\) (v^ — XJ) 

RMN sera donc, dans tous les cas possibles, un produit de fac- 
teurs de la forme (i) multiplié par un certain nombre de quan- 
tités jr, y, z. 

Élément linéaire en coordonnées elliptiques, — Les sur- 
faces p =C^, [JL==C^, V = C'® forment un système triple ortho- 
gonal et l'élément linéaire peut s'écrire : 

ds' = a}d;J + P'diL' + fdy\ 

a, p, Y étant des fonctions de p, j;., v que nous allons déterminer; 
on a 

ds' = dx^ + df/^ + dz\ 

dx pdp [jL^ix vc/v 



X p^ — a^ iK^ — a^ 

dy pû?p )xd^ 



dz pdp I [kd\u 



y' 


«•^' 




v/fv 


V* 


p ' 




vdy 



z — p» _ c^ "^ [x-^ — c" ' v^ — c^ 
On a donc 

[222 

[X^ |,2 ^2 -1 

+ vwf — ^^ + - y' H ^^1 
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et ensuite 

Je vais poser 

_^. (P^-^^)(p^-^')(P^-0 



B' = 



(l-^- 


-«^)(f^ 


p-^ 


i^)(p^- 


-c^) 


(v^- 


- o') (v' 


f*' 


i^) (v' - 


-c») 






v^ 







A^ et C^ sont positifs, B^ est négatif; donc A et C sont réels, 
B est purement imaginaire. 
Avec ces notations on a 



a 



?= 



V 
i 



Q 
Q 



yv^ — p*est purement imaginaire, donc |3 est réel. 

Equation de Laplace. — Voyons ce que devient Téquation 
de Laplace en coordonnées elliptiques. Nous avons déjà établi 
(p. 4o)> qtie si on désigne Télément linéaire par 

on a 






L'équation AV = o devient donc : 



IM^^h"- 
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On a 

BC y i^V^ - 0^) {p^ - a') 

BC 
et l'équation de Laplacc est alors : 






A seul dépend de p. Ou a donc 



en multipliant par ABC on obtient Téquation 

2A(.--v-)i[.^].o 

Les coefficients de cette équation sont réels, car A", B^, C" sont 
réels, et par suite leurs dérivées le sont aussi. 

Fonctions de Lamé. — Parmi les solutions de l'équation de 
l.aplace, il y a des polynômes de la forme RMN : on les appelle 
fonctions de Lamé. 

Remarquons d'abord que nous avons les identités évidentes : 

^(i^' — v') = o; 

2j'{\J? — ^r) = o, 
Il en résulte que l'équation (2) est équivalente à la suivante 
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Supposons que l'on aitV= RMN et que R satisfasse à l'équa- 
tion 

M et N satisferont à des équations analogues, et la fonction V 
satisfera h Téquation de Laplace. Cette condition suffisante est 
aussi nécessaire, c'est-h-dire que si V est de la forme RMN et 
satisfait à l'équation de Laplace, R satisfait à l'équation (3), et par 
suite M et N satisfont aux équations analogues obtenues en rem- 
plaçant p par [L et v, et A par B et C. 

En effet si on fait u=v l'équation doit être encore vérifiée ; 
R se réduit a M en ce cas, et l'équation se réduit à 

_^„ rf^M ^ dB dM . „ ,,, -.-, 

B' -r^ + B --z 3 h Hu^M + KM = o. 

Existe-t-il de telles fonctions? Considérons la fonction/? (//) 
définie par l'équation différentielle 

p'u = 2 \/{pu — e^) {pu—e^) [pu — e^) , 

où 

o. 





^1 -t- ^2 -t- ^3 — < 


Posons 






p- —pu + h, 




a} e^ +/^ 




h- e^ + A, 




c' -.e^-\-h; 


on en conclut 





n= . 

3 

On a donc 

p'u = a \/(p^-a^)(?^-^l(p'-c-^) 

-=2Ap. 

On a d'ailleurs 

, du 



I 

i donc 
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A= '' 



du 



et par suite 



rfR _ rfR rf? _ dR 



L'équation 



devient 



do do du du 



^wbwh^"''+^^''''° 



^+(Hp.+ K)R = ., 



Peut-on choisir H et K de manière que R soit un polynôme 

en p, ou un tel polynôme multiplié soit par yp* — a', soit par 

\/p^ — b\ soit par \/^'^ — c^? 

Ce polynôme, si on y remplace p^ par p (w), sera une fonction 
doublement périodique admettant les deux périodes de p[ii) 
2a>j et 2(i>2. Elle admettra une double infinité de pôles qui seront 
les pôles de p (w), c'est-à-dire u = a/ww^ + 2 /i(o^, m et n étant des 
entiers positifs, négatifs ou nuls. Ces pôles qui sont des pôles 
doubles de/? (m) seront pour le polynôme des pôles d'ordre aw. 

De plus y p* — à^ = \ P^^ — ^u ^^t une fonction doublement 
périodique admettant les périodes 2Wj et 4^^ et ayant comme pôles 
simples les pôles de p (w); donc si on ajoute a w la valeur 2(0p 210, 
ou 2(^3, la fonction R garde sa valeur ou change de signe sui- 
vant qu'elle contient o ou 2 radicaux, ou bien qu'elle est le pro- 
duit d'un polynôme par i ou 3 radicaux. 

Donc R est une fonction paire, ou une fonction impaire, mais 
dans un cas comme dans l'autre on doit avoir le développement 

T^ ^A ^4 ^O 



«» iC-"^ z^"-* 



si on porte cette valeur dans l'équation, en se rappelant que 

p[u)=—^ -i-dll^ 
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on obtient l'équation : 

[n{n-^i)a, , (n — i) {n — 2)a^ , "1 

Pour que cette équation soit vérifiée il faut d'abord que le 
terme en ^ disparaisse, ce qui n'est possible que si l'on a 

Il = — n [n -{- i). 

On a d'ailleurs si q est le degré de R, et qu'il faille lui ajouter p 
radicaux, (p = i , 2 ou 3) : 

n = iq ■\'p. 

Formation des fonctions de Lamé- — Ainsi donc les poly- 
nômes de Lamé, doivent s'ils existent satisfaire aux conditions 
suivantes. 

Regardés comme fonctions de ^, y, z ce sont des produits de 
facteurs quadratiques multipliés par un ou plusieurs des facteurs x^ 
y^z. Si on remplace^,?/,2;par leur valeur en fonction de p'^^ji."^, v^, on 
obtient une fonction de la forme RMN, où R désigne, soit un poly- 
nôme en p^, soit un tel polynôme multiplié par 1,2 ou 3 des 
radicaux 

v/p=^ \/7^^, \/7^^^^ 

et M, N désignent les mêmes fonctions ou p est remplacé par [jl 
et V. 

On a pour Q huit formes possibles, h chacune desquelles cor- 
respond une forme pour R. Si n est le degré de cp, le degré de R 
en p^ sera suivant le cas 



n n — I Jl — 2 n — 3 



, ou 



2 2 



suivant que R contiendra o, i, 2 ou 3 radicaux. 

Enfin si dans R on remplace p'^ par pu^ on obtient une fonction 
qui, lorsqu'on remplace u par a -f- sto^, change son signe ou reste 
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invariable, suivant qu'il y a un nombre impair ou pair des deux 

radicaux y/p^ — i^ety/p'^ — c", car \^pu — e, ne change pas de 
signe pour u-\-:hù^. 

De même quand on augmente u de 2(0^, y /j/^ — e^(il\pn — e^ 
changent de signe et par suite R changera de signe ou non, sui- 
vant qu'il contiendra un nombre impair ou pair de ces radi- 
caux. 

Ces résultats sont résumés dans le tableau suivant, oii la pre- 
mière colonne contient la forme de Q, la seconde la forme corres- 
pondante de R, la troisième le degré du polynôme qui entre 
dans R en fonction du degré de Q, la quatrième et la cinquième 
indiquant si la fonction change de signe pour l'adjonction à u des 
périodes aco^et 2w^; Q^ désigne uniquement un produit de facteurs 
quadratiques. 




R ne contient pas de radical 
RconlicnH/p- — a^ 



» 
» 
» 
» 



n 

2 



/ n — I 

1 

\ 

\ 

I n — 1 




+ 
+ 

+ 



+ 




+ 






1^ 
+ 



/ »+ 



•2 

n — I 



Détermination de K. — Ceci posé, cherchons h déterminer K. 
Pour les fonctions du premier genre, R est un polynôme 

d'ordre contenant h i coefficients homogènes. 

2 2 

Substituons dans l'équation différentielle, nous avons un poly- 
nôme de degré \- ly qui devra être identiquement nul. Le 

coefficient du terme de degré le plus élevé disparaît de lui-même 



à cause de l'hypothèse II = — n [n-\-i). Il reste alors 
coefficients à faire disparaître. 



n 



+ I 



n 



Les-^^^^ 1- I coefficients du polynôme doivent donc satisfaire 
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a \- i équations linéaires et homogènes, dont les coefficients 

dépendent linéairement de K; l'élimination des coefficients du 
polynôme donnera donc pour K une équation de degré 1- i . 

Si R est du second type, on aura à déterminer coefti- 

cients d'un polynôme de degré -. Le degré de l'équation 

en K sera donc, d'après un raisonnement analogue au précédent, 
n+ I 



Pour ceux du troisième, l'équation en K sera de degré — . 

Pour ceux du quatrième, l'équation en K sera de degré . 

Ces résultats figurent dans la sixième colonne du tableau pré- 
cédent. 

Il résulte de ceci qu'à chacune des racines des équations en K 
ainsi obtenues correspond un polynôme et un seul. 

Si n est pair, il n'y a de polynômes que de première et de 
troisième espèce ; leur nombre est 



n ^ n 

h I + J = 2/2 + I . 

2 2 



Si n est impair il n'y a de polynôme que de seconde et quatrième 
espèce ; leur nombre total 

n — I , „ /i 4- I 

' v) ■ = 2/1 + I . 



Donc il y a au plus 2/2+ i polynômes de Lamé. Je dis au plus, 
car les équations en K pourraient avoir des racines multiples ; 
nous verrons qu'il n'en est rien. 

Analogie avec les fonctions sphériques. — Supposons que 
x\ y, z soient infinis du premier ordre ; [jl et v restent essen- 
tiellement finis, mais p deviendra infini du premier ordre aussi. 

L'ensemble homogène des termes de degré le plus élevé dans 
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la fonction de Lamé, satisfait évidemment a l'équation de Laplace. 
C'est donc un polynôme sphérique que Ton peut écrire II (j:*,?^,.:;). 
Le terme de degré le plus élevé dansRMN en p sera ap"'MN. 
Passons aux coordonnées polaires 

'' X = r sin 8 cos '^, 
j y ^=^r sin 8 sin cp, 

\ z =7C0sQ, 

si p est très grand il difTère peu de r ; ul, v dépendront au contraire 
de c5 et de 8. 

En effet l'équation qui détermine |jl peut s'écrire 

' '' • — 1=0. 



■^ ^'^ ' ^' — b^ ' ^^x'—c' 



[JL" — a \x 



Si Xyi/,z sont très grands, l'unité est négligeable; il reste une 
équation du second degré ; \l et vies deux racines non infiniment 

grandes ne dépendent que du rapport-^ et —c'est-à-dire de cp 

X X 

et de 6. 
On peut alors écrire 

/•"Il (sin 8 cos cp, sin 8 sin cp, cos 8) = ap"MN, 

et comme r = p, on a 

n (sin 8 cos cp, sin 8 sin cp, cos 0) ^=aMN ; 

le produit MN apparaît donc comme une fonction sphérique, 
etllseradécomposable en un produit de facteurs quadratiques (*). 
M. Moutard a démontré que c'était le seul cas où un polynôme 
sphérique était décomposable en un produit de facteurs quadra- 
tiques et de facteurs linéaires. 

MN, considérée comme fonction de 8 et de cp, est donc une fonc- 
tion sphérique : il en résulte que toute fonction de [x et v peut se 
mettre sous la forme d'une somme 



2«MN, 



(^) Voir Journal de Liouville, t. XI, i8 îG, p. -278. 
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puisque toute fonction peut se mettre sous la forme d'une somme 
de fonctions sphériques; mais l'analogie des fonctions spliériques 
et des fonctions de Lamé, va apparaître davantage quand on étu- 
diera le cas où les surfaces homofoeales sont de révolution. 

Cas des ellipsoïdes de révolution. — Cela peut se présenter 
de deux manières différentes suivant que Ton di a = h, o\x h = c. 

1° Si a= hj on a des ellipsoïdes de révolution aplatis quand 
A- est supérieur à a^ et des hyperboloïdes à deux nappes si 
<^r<).^<6'^. Les hyperboloïdes à une nappe se réduisent à des 
plans passant par Taxe 0^. 

o,° Sib = 6*, on a des ellipsoïdes de révolutions allongés autour 
de O.r, pour ).^ > a\ Si a' > A^ > P on a des hyperboloïdes à une 
nappe, les hyperboloïdes à deux nappes se réduisant a des 
plans passant par Taxe Ojc, 

Dans ces deux cas, soit [ji, soit v ne peut servir de paramètre, on 
est obligé de prendre comme paramètre l'angle du plan avec un 
plan fixe, le plan xOz par exemple dans les deux cas. Nous 
avons : 

a: V (p^ — a') [a' — ^') {a' — V') ^ b' — c' ' 

Si b tend vers a, [jl tend aussi vers a, et l'on voit que Ton a : 



-^==limVS 5- 

Si cp est l'angle du plan passant par le point [x^ y, z) et Taxe O^ 
avec le plan xOz, on voit que jjl peut être considéré comme fonc- 
tion de îp, et l'on a 

\i} = a} cos^ cp -+- J"^ sin^ cp. 

On voit que l'on a 

^ja"-—)}} 



^^"Wt^ — k = 



., ,., — ces % 



CAS DES ellipsoïdes DE RÉVOLUTION ii^ 

Dans les mêmes conditions Téquatlon qui détermine v peut 
s'écrire en supprimant l'unité, négligeable devante, y etr, 

f!±il + -f!- = o 

V' — a' ^ ^' — c' 
Si on passe aux coordonnées polaires, on a 

r^ sin^ô r^ cos'^9 



+—> — :t- = o ; 



et ensuite 



2 2 ' 2 t 

V — a V — c 



in Q=\/ —:; T-» 



sin 



cos 6 



V cr—c- 



Dans ce cas on peut immédiatement former les fonctions de 
Lamé. 

Nous savons que MN est une fonction spliérique ; M ne dépend 
que de cp, N ne dépend que de ; MN est donc l'une des fonctions 
sphériques fondamentales considérées plus haut (p. 42). 

On peut alors écrire M = cosy;'^ ou sin jy^ \ et pour Non 
prendra 



en posant ; 



rj;(cosO)=F|;[v/^] 



On voit que les zéros de N se déterminent en calculant les 
zéros de F}[, abstraction faite des zéros ± i. Ces zéros sont tous 
réels et distincts ; car, comme on le voit par le théorème de 
RoUe, les zéros de Fj; séparent les zéros de Ff,~^ On a 






Donc si p est pair, FjJ ne changera pas de signe avec sin 9; au 
contraire sip est impair, F,,^ changera de signe avec sin 9. 

Si on remplace x^ijy i^par leurs valeurs en coordonnées polaires, 



126 



FONCTIONS DE LAME 



Q changera de signe s'il contient l'nn des facteurs x ou y ; 
mais s'il n'en contient aucun, ou s'il en contient deux, Une chan- 
gera pas de signe; p sera donc impair dans le premier cas, et 
pair dans le second. 

On voit aussi que y et sin /; '^ changent de signe en môme temps; 
donc si Q contient le facteur y, M devra être pris égal à sin/;'^, 
sinon M devra être pris égal a cosprf. 

De même si ù^ = c^, on prendra comme coordonnées polaires 

X = /'cos 0, 

y = r sin 6 cos cp, 

z = r »in9 sin es. 

La fonction N sera sinp'f ou cosy;'^, la fonction M sera 

FS (cos e) . 

On prendra y^ impair si Q contient l'un des facteurs y on z; 
p sera pair au contraire si Q ne contient aucun de ces facteurs ou 
s'il les contient tous deux. On prendra pour N, sin p ^ ou cos y? cp 
suivant que Q contient ou ne contient pas le facteur z. 

Résumé. — Le résultat de cette discussion est résumé dans le 
tableau suivant : la première colonne contient les 8 formes possi- 



V 



Pxj. . 

PXsj . . 

V X zx . . 

V xxy . . 
Vxxyz. . 



R 






Ry^p2__c2 

R^/(p2_c-^)(p2_^2^ 



a- = ^' 



o 



^') 



Cos 



Cos 



Sin 



Cos 



Sin 



Cos 



Sin 



Sîn 



/>2 m c2 



Cos 

Cos 

Cos 

Sin 

Sin 

Sin 

Cos 

Sin 
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bles de Q, la seconde contient les formes correspondantes de R, la 
troisième et la quatrième contiennent les cas limites pour d^ = i^; 
dans la troisième on a mis le signe o si p est pair, i si p est 
impair ; la quatrième indique si M se réduit à cos po ou sin po ; 
la cinquième et la sixième contiennent les indications corres- 
pondantes dans le cas où b'^ = c^. 

Les valeurs de K sont réelles. — Nous allons maintenant 
démontrer que les équations en K ont toutes leurs racines réelles. 

Supposons d'abord a' = P. On sait alors trouver les valeurs 
de K. Nous pouvons poser ^^ = a' cos- o + ù^ sin" cp, puisque jjl 
est compris entre a et 6 et lorsque a = b^ ceci est encore la 
valeur de M ; Téquation qui définit [jl 



(^-^)-t«(" + ')l-^-K]M = o, 



d^ \ d^f. 
peut s'écrire, en prenant comme variable o 



o. 



2 2N d'M , dU 

(u — c^) — r~i — \- 'Jr — a-) sin es cos 'o — ; — 

— [/2(/z + i)l^'— kIm =0 
et si on fait a^ = P on trouve l'équation 

[a^ — c^) -^ _r^(,e + iy _k1m=o. 

Cette équation doit avoir comme solutions cosp'f et sin/;'p, on 

a donc 

K — n(n-{-i) 2 

3-^^ Ty =P ' 

a" — c ^ 

L'équation en K a donc /e -|- r racines réelles, correspondant 
aux valeurs de /?, 

o, I, 2, . . . (/i — i), n. 

On voit que les équations en K ont toutes leurs racines réelles ; 
les valeurs paires de/7, séparent les valeurs impaires. 
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Lorsque b variera de a vers c les racines des équations corres- 
pondant aux deux groupes de fonctions de Lamé, contenant 
zéro et deux radicaux d'une part, sépareront les racines corres- 
pondant aux deux autres groupes contenant un ou trois radicaux ; 
cela ne cessera d'être, que si une des racines d'un groupe d'équa- 
tions, devenait égale à une racine cp d'un autre groupe. 

Soient, K cette racine, R et S les fonctions de Lamé corres- 
pondantes, appartenant par hypothèse à deux groupes diffé- 
rents : on a par hypothèse 

lV'_(IIp2+K)R=o, 
S'' — (IIp^ + K)R = o; 
on en conclut 

RS" — S''R=o. 

On peut intégrer immédiatement et écrire 

RS' — S'R = a. 

Si on suppose R et S exprimés en fonction elliptique du para- 
mètre u^ on sait que l'on pourra augmenter u de 2 o), de 2 (o^ ou 
de 2 (lig de façon que R par exemple ne change pas de valeur, et 
que S change de signe, et de même pour R' et S'; mais alors le 
produit RS' — SR' aurait changé de signe ; on a donc a == o ; 

on en conclut que -^ est une constante, ce qui est évidemment 

faux. 

Deux équations en K correspondant à une même valeur de n 
ne peuvent donc admettre de racine commune. 

Remarquons que dans l'équation en K le coefficient du terme 
de degré le plus élevé se réduit k l'unité; il est donc impossible 
que le degré de l'équation s'abaisse. 

Les fonctions de Lamé sont linéairement indépendantes. 
— Je vais démontrer maintenant qu'entre les p fonctions de 
Lamé, d'une même forme, et pour une valeur de n donnée il 

n'existe aucune relation linéaire à coefficients constants ; Ip est 

, , , , n n — i n -{- i n\ 

1 un des nombres 1- i , , ■ , — ) . 

2 2 22/ 
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En efTet : supposons que Ton puisse avoir une relation de la 
forme 

en portant R^+idansTéquatlon du second ordre, on aura à écrire 
que les p coefficients d'un polynôme seront nuls ; on aura pour 
déterminer a^ a^... a^. p -\- i équations du premier degré, dont la 
première sera satisfaite identiquement. Il restera donc p équa- 
tions homogènes à q inconnues. L^élimination de ces inconnues 
donnera un certain nombre d'équations en K de degré au plus 
égal à y et ces équations devront être satisfaites pour les p va- 
leurs de K ; ces équations se réduisent alors à des identités. On 
peut donc déterminer les coefficients a^, a^... a^ quel que soit 
K, et on aurait alors un polynôme comme solution de Téquation 
de Lamé, et cela quelque soit K; c'est une contradiction. 

Les polynômes de Lamé ont leurs racines réelles. — Je 
vais maintenant démontrer qu'un polynôme R quelconque a 
toutes ses racines réelles et comprises entre a^ et c^. 

On peut partager l'intervalle de — 00 à + 00 en quatre autres 
par les nombres a^, P et c^. 

Je vais d'abord démontrer que le nombre de racines compris 

dans chaque intervalle est constant lorsque i* varie de cr à c^. En 

efTet, Il ne peut en être autrement que si, pour une valeur de i, 

deux racines réelles devenaient imaginaires. Mais alors, pour 

cette valeur de è, ce polynôme aurait une racine double : R' serait 

nul. On a 

R''=FR = o, 

F étant (IPp +K); on a 

R''/ = FR'+F'R, 
R^v =FR''+ 2F'R'H- FR'' etc. ; 

toutes les dérivées seraient nulles, ce qui est impossible. 

Il pourrait encore arriver que Tune des racines devienne égale 
ù a^, i2 ou c^. Supposons par exemple qu'une racine devienne 
égale à a^ : comme on a posé 

PoiNCARÉ. — Figures d'équilibre. y 
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la valeur de u correspondante serait w^, or on ay/o>j = o, o)j est un 
zéro de p'tù^ ; donc ce serait un zéro de R' et l'on revient au cas 
précédent. 

Ceci posé soit e^ le nombre des racines imaginaires de Tensem- 
ble des polynômes correspondant à une valeur de n ; 

Sg, le nombre des racines réelles comprises entre a* et b^ ; 

Sg, » » réelles comprises entre h^ et c^ ; 

Êj^, » » réelles non comprises entre a^ et c^. 

Le nombre total des racines est égal à la somme des degrés des 
polynômes 



Tt l n \ . , 



Lorsque h^ est égal à a', toutes les racines t^ se réduisent h a^ 
et s'g se réduisent à a* ou à c^\ lorsque 6^ est égal à c^, toutes les 
racines £3 se réduisent a c^ et s^ se réduisent a a^ ou à c^. 

Mais si 0} = è^ les polynômes sont les dérivées successives 
de (v^ — i)** et ils ont toutes leurs racines réelles ; il y en a un qui 

a o racine, un qui en a une, un qui en a deux, ... un qui en a — 
abstraction faite des racines it i ; on a donc 



n [ n \ , 



de même on a 



^2 



n / n \ 



Nous avons donc l'égalité 



^ + ^2 + S3 + £, =-^ ^-^ + l\ + t, +£^+£'2+ s'aî 



d'où l'on conclut : 



et par suite 



e^ + £4 + e'2 + ^'3=o- 



ei=S,= £'3 = £',=0. 



Le théorème est donc démontré. 
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Pour savoir maintenant combien il y a de racines dans chacun 
des intervalles a*... i*, i*... c*^,on voit que dans le cas général il y 
en a autant que dans le cas particulier, où a^= i*. 

Et on voit qu'il y a : 

un polynôme ayant o racine entre a^ et i^, et — entre b^ et c*. 



)) 



)) 



I 


» 


n 

2 


» 


n 

2 


» 

1 





» 



Développement d'une fonction en une somme de fonctions 
de Lamé. — Soit un ellipsoïde de la famille, déterminé par le 
paramètre p ; si je pose 



/ = 



V/(?^-l--^)(?^-v^)' 

AMN.M,NjrfT = o. 

MN désignant un des produits déjà désignés, M^N^ en désignant 
un autre, et l'intégrale étant étendue à toute la surface de Tellip- 
soïde. 
Soit 

V = RMN, V, = R,M^N,. 
Puisque Ton a AV = o, la formule de Green donne 




V-^— V —]d^ = o 
dn ^ du I ' 



l'intégrale étant étendue a une surface fermée quelconque ; sup- 
posons que cette surface soit Tellipsoïde ; R et R^ sont alors des 

constantes. Calculons «^ et -7-^ : quand on se déplace sur la 
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normale à rdlipsoïde, jjl et v demeurent constants, p augmente 
de rfp ; et on a 



dN 
dn 


ÔV dp ÔV du 
dp dn du dp 

dV du 
àu dn ' 


rfp 



mais nous avons 

- — = mi^ --7— ; 
OU au 

rîntégrale devient donc 

MN. M,N, (RR^, — R,R') ^dfj = o; 

mais RR'^ — R^R' ne dépend que de p et comme c'est une cons- 
tante, on peut la faire sortir du radical. On a d'ailleurs 

RR', — R^R'z^o; 

il vient donc en changeant le signe 

— MNM^N, ^dfj=o. 



Tout revient h calculer 



du du dp 



dn dp dn ' 

nous avons vu que l'on avait, en posant p^ =pu-^ h 



du 
d'où 



do , 

au 



lorsque p tend vers l'infini, pu tend vers oo , u tend vers o, donc 
décroit ; on doit prendre le signe moins. 
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On a d'ailleurs, puisque l'élément linéaire est 

* 

et que p seul varie sur la normale : 

^^j Av/([JL^ — y-^) 

dn d Q 

On a donc, en posant 



dn ' 



V/l^^ — v^ 



ce qui est l'égalité annoncée. L'Intégrale, d'ailleurs, ne dépend 
pas de p. Car nous avons 

et par suite 

^~BCv/(p'^-|^o'^-v^) 

~ BC 

On a démontré que toute fonction déterminée sur l'elllp- 
soïde peut se mettre sous la forme d'une somme de fonctions 
de Lamé; 

4. =2| A* M* ^'*- 

Le théorème précédent permet facilement de trouver les coef- 
ficients. On a 

A* M A Na: rfcr = r V k^Mu N;t M, Nj rfa ; 



I 



e second membre se réduit i( Af^ l /M^Nfrfo" : ce qui détermine Aj 
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Les fonctions S. — Nous allons maintenant considérer une 
autre fonction de p satisfaisant aussi à l'équation de Lamé; nous 
la définirons ainsi qu'il suit. 

Soit R une fonction de Lamé de degré n définie par la relation 

R"+FR = o; 

la fonction S, associée a la fonction R, sera définie par la relation 

S'' + FS = o 
où F a la même valeur : on en conclut que 

S^'R — R"S = o ; 

en intégrant on trouve 

S'R— RS'=C% 

et nous prenons arbitrairement cette constante égale a 2/i -j- i ; 
on a donc 

S'R — RS' 2/1+1 



R^ R' 



Donc 



S I 2/1 + 1 - 



on peut choisir arbitrairement la limite inférieure d'intégration 
et on a 



Pour p^ très grand 

R = Ap», 

A étant une constante que nous pouvons prendre égale a Tunité ; 
Il est alors très petit : puisque 
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p^ est approximativement égal à — j 

I 

p » » — 

* Il 



On a donc 



R = z^~« , 



S = M-» I -V— du == 11''+^ ; 



S est donc sensiblement égal à ^i'* + * == ^^^ ; le produit SMN 
satisfait aussi à Téquation de Laplace. 



CHAPITRE VII 



ATTRACTION DES ELLIPSOÏDES 
FIGURES d'Équilibre 



Problème de Dirichlet sur r ellipsoïde. — Abordons maîii- 
tenant le problème de Dirlchlet pour un ellipsoïde, dans lequel 
on se donne, sur un ellipsoïde déterminé par le paramètre p^, 
les valeurs d'une fonction V, en proposant de trouver la fonction 
V satisfaisant a l'équation de Laplace et prenant sur relUpsoïde 
la série des valeurs données. 

Soit sur Tellipsoïde 

V==4.=^A,M,N,; 

on peut poser 

A,= a,R?S?, 

RJ étant la valeur pour p = p^ de la fonction R associée à M et N ^ 
et S^ » * » S,- » K,-- 

On a alors 



*=2a,R?S?M,N,. 



Ceci posé, considérons la fonction définie par l'expression 



S' 



S'RMN 



à l'intérieur de l'ellipsoïde, et 

aR'SMN 



S' 



à l'extérieur du même ellipsoïde. 
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Ces deux séries sont convergentes dans les espaces indiqués, 
la fonction est donc bien déterminée. 
Calculons les quantités 

OV ÔV 

et 



on a 

ôV vi ^« ^, ,, ,, du 



L=^a,RJS;M,N,. -£. 



= /y^a,R?S7M,N,. 
On a de m A me 



0/z. 



= — ZV a,S;R'?M,N,, 



V est une fonction continue quand on traverse la surface, mais 
ses dérivées ne le sont pas. V est donc le potentiel d'une masse 
superficielle, dont la densité o est déterminée par la relation 



Au lieu de supposer la densité variable, on peut supposer une 
distribution de densité égale à i, mais d'une épaisseur IJ, telle 
que Ton ait 

ôV dV _ ' 

'^ ^,—~^''^'^ 

équation qui détermine Ç, car on a 

4«J; = - l^/i (R? S'î - R'? s») M,N, 



= zVa(2rt+i)MN. 



Si l'épaisseur de la couche est dctermiuce par une fonction de 
la forme 



E 



Z3MX , 
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on en déduira réciproquement le potentiel sur la surface par la 
formule 



y 1 MN; 

jLJ 2/1+1 



d'où Ton déduit simplement le potentiel h Tintérieur et h Texlé- 
rieur de la surface. 

Application. — Voyons maintenant les plus simples des fonc- 
tions de Lamé. 

I® pour 71 = on a 

S= u. 
2° pour M= I on a trois fonctions R 

R. = n/P'^^ r, = v/p^^=^ïs R, = \/F=^- 



On a 



Z = //3R3M3N3, 



en posant 



h.= 



3° Pour /ï = 2, on a : 
d'abord deux polynômes du second degré 

p2 — h' 
et 
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que Ton peut calculer facilement; on sait que Ton a 
et ensuite les trois fonctions 

On a 

R^ = R2R31 

R5 = RiR3, 
Re = RiR2; 

le produit R^M^N^ est à un facteur constant près ?/c. 
» R^MjNg » » zx, 

» RjMjNj » » xij. 

Remarquons encore la signification des fonctions R^, R^, Rg, 
qui sont les demi-axes de l'ellipsoïde correspondant à p. Son 
volume est 

T = — t:R^R.^R3= — ttR^Rj^. 

Supposons qu'à la surface de l'ellipsoïde, soit répandue une 
couche de densité constante, et d'épaisseur 

le potentiel à la surface sera 

les indices inférieurs indiquant le numéro d'ordre de la fonction 
de Lamé, l'indice supérieur indiquant qu'on y fait p = p^. A l'in- 
térieur, le potentiel sera 

47:SJR„MA=4^'/o, 

le potentiel sera constant à l'intérieur; on peut donc consi- 
dérer la distribution de matière comme la distribution d'une 
couche électrique en équilibre à la surface de rellîpsoïdc sup- 
posé conducteur. 
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Le potentiel sera à l'extérieur 

Les surfaces équipoteiitielles seront donc les surfaces u = C*®, 
c'est à dire les surfaces p = constante; ce seront des ellipsoïdes 
homofocaux. 

La couche correspondante est celle comprise entre deux ellip- 
soïdes homofocaux comme le montre la valeur de J^. 

Attraction d'un ellipsoïde homogène. — Considérons (fig. 17) 
un ellipsoïde homogène dont le potentiel soit V. Déplaçons Tel- 

lipsoïde .parallèlement à Taxe 0^, d'une 
quantité c; le potentiel en un point (.r, y^ z) 




P" 
P' va devenir 



\' (^x — t,y, r.)=V {x,y,z)—t 



ÔV 



et l'ellipsoïde prendra la position E'. La 
différence de potentiel provient de ce qu'il 
y a en plus la région ombrée a droite, et 
Fig. 17. en moins la région ombrée à gauche. On 

peut supposer que c'est le potentiel d'une 

épaisseur Ç répandue à la surface de l'ellipsoïde, J^ étant positif 

à droite et négatif à gauche. 

Calculons J^ : si PP" est la normale commune aux deux ellip- 

soides, Ç = PP"; PP' = e, et le triangle PP'P" est rectangle 

en P'', on a donc 

PP"=ecos(PP',PP'0. 
Comme 

PF' cos (PP, PP'Q = projection de PP'' sur O.r, 
et, puisque PP" est normale a l'ellipsoïde p, on a 

Z cos P = -^ — rfp 

ô? ' 

— ' an 



dp dn 



^x dp y 
dp dn 
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de ces égalités on déduit 

^ dj: do Oj: du . d.r 

COS P = -r r- = -t:^ r— = — l 



dp dn du dn du 

Nous savons aussi que : 



du 
en posant 



=ar;m,\, 



* Op du 



— — X-L 



et comme 






on en conclut 



et l'on a 



IV, = R, 



'^•"^^ =AR,M,NV 



d// 



Il vient donc 



cosP = A/R;M,N,, 



On peut donc maintenant déterminer le potentiel V relatif h 
la couche d'épaisseur t^, en appliquant la formule connue : on a à 
la surface 

V'=— 4^ eARJSJRÎM.N.. 
A rintérieur de rellJpscïJe on a 
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et, puisque R^M^N^ =^, 



T=- 



r; ' 



on a 



donc 






r; 



A Textérieur on a 



ÔV 4it/i 



àx 



= ^^ rîr;s,m.n, 

"RT' 



ÙV 



A l'intérieur -r— est donc proportionnel à ^, à l'extérieur la 

relation est plus compliquée; -r^ représente la composante de 
l'attraction parallèle a Ox, 

m 

Théorème d'ITvory, — Considérons deux ellipsoïdes homo- 
focaux, correspondant, le premier E^ au paramètre p^, le second E^ 
au paramètre p^; p^ > p^ (fig. 18). Soient P^ et Pq deux points cor- 
respondants, c'est a dire ayant 
respectivement pour coordon- 
nées elliptiques (p„ ix, v) (p„ 
[ji, v) . Les coordonnées rectan- 
gulaires de deux points corres- 
pondants sont proportionnelles, 
c'est-à-dire que l'on a 




Fig. 18. 



rz 



1' 



p^ y, /• étant seulement fonctions de p^ et de p^. 

Le théorème d'Yvory consiste à comparer l'attraction de l'ellip- 
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soïde E^ sur P^ à Tattraction de Tellipsoïde E^ sur P^. La pre- 
mière a pour composantes suivant Ox 

47rARlSlR;M,N„ 

la seconde a pour composantes suivant Oj; 

47rARîR;S}M,N,, 

le quotient est donc 

R® 

Ri —^ ' 

on trouverait de même que les rapports des attractions sui- 
vant Oy et 0^ sont respectivement 

Rî ^, RJ 

et 



\\ï K 



Ellipsoïde de Mac Laurin, — Considérons un ellipsoïde 
homogène, que nous faisons tourner autour de Taxe Oo^", d'une 
rotation infiniment petite to ; le point de coordonnées x^ y, 3, 
aura pour nouvelles coordonnées 

son potentiel sera 

/ dV dV \ 

V(^, y + ziù, z—7jiù)=N{x, y, z)+iù\z~ y-^T)^ 

V — V est le potentiel dû à la couche superficielle comprise 
entre les deux ellipsoïdes. 

Si [x^ y^ z) représente un point intérieur, Faction en ce point 
a pour composantes 

T.rS; T//S5 T^SS 



Rî ' RS ' R5 ' 

elle dérive d'un potentiel 
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en appelant V^^ le potentiel au centre ; et cette formule est encore 
vraie h la surface : donc 

dV dV 

et 



sont proportionnelles à y et à ;;, et on a 

V — V = kiùyz 

En appliquant la formule démontrée précédemment (p. iS^), 
l'épaisseur de la couche en un point est 

5A(o/ M,N, 



47t s; 

Pour que l'ellipsoïde puisse être une figure d'équilibre, il faut 
qu'à la surface ou ait 



co^ 



Nous allons traiter un problème plus général ; de savoir si un 
ellipsoïde à trois axes peut être une figure d'équilibre, sachant 
qu'il est soumis, outre les attractions mutuelles de ses parties, a 
une force dérivant d'un potentiel 



a, p, Y étant des constantes. 

Pour que l'équilibre ait lieu, il faut que 

soit constant à la surface de Tellipsoïde dont l'équation est 

1 '- 1 I = o 

Kî ^ R^ ^ Ri 
Il faut donc que Ton ait 

(I) «RÎ-TS.R. = 13Ri-TSA=TRi-TS3R3. 



y. 



' .^ ■",. :• - 
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Dans le cas particulier qui nous occupe, c'est-à-dire quand 
Tellipsoïde est en rotation, on a, a = o, ^ = y == cj^, et nos 
formules deviennent 

TS,R, = w^R* — TS,R, = iù^Ki — TS3R3 ; 

on a donc 



/JJN —,=z^ ^fiz—^i^X __ Q RgSg — RjS 



1^1 



T ^ R^ ^ R| 

et d'ailleurs 

T= — RjR^Rg; 

T est une donnée du problème où les inconnues sont p^ — a*, 
— b\f — c\ 

On a une première condition, qui est que tù^ doit être positif; 
donc il faut que 

R,S, — R,S.>o, 
R3S3-R,S.>o. 

Ces inégalités sont toujours vérifiées. Démontrons par exemple 
que la première Test, c'est-a-dire que Ton a 

R.S^— R,S,>o. 

Remarquons que R^, R2, S^ et S.^ sont des quantités positives, 
puisque p^ > a^ ; il faut donc que Ton ait 



s, ^R. 



s. R. ' 



c'est-à-dire 




Rf 



PoiNCARÉ. — Figures d'équilibre. lo 
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du T^T 

Ri 



et on est ainsi ramené a un théorème d'arithmétique connu. Le 
premier rapport est compris entre la plus grande et la plus 
petite valeur du rapport des quantités placées sous le signe/: la 
plus grande de ces valeurs est i correspondant a w = o, p = oo, 

^2 2 

CL 

la plus petite est \ ji', le théorème est donc démontré. On 

en donnera plus loin une démonstration plus générale. 

On voit de suite que Taxe de rotation ne peut pas être Taxe 
0?/, car il faudrait que R2S2 soit la plus grande des trois quantités 
analogues. 

La relation (II) a une solution évidente, c'est R2 = R3, qui con- 
duit à un ellipsoïde de révolution autour de Ox : c'est l'ellip- 
soïde de Mac Laurin. 

Ellipsoïde de Jacobi. — Les relations (II) peuvent aussi être 
satisfaites pour un autre ellipsoïde découvert par Jacobi, mais 
nous devons, avant de le rechercher faire la remarque suivante : 

Si nous nous donnons un ellipsoïde homogène quelconque, 
nous pouvons toujours choisir les constantes a, p, y, de manière 
à satisfaire aux relations (I) c'est-à-dire de façon que cet ellip- 
soïde soit en équilibre soùs Faction de la gravitation et d'une 
force supplémentaire dont les composantes sont eux, P^, ys ; nous 
allons supposer a = o. 

La condition d'équilibre d'un ellipsoïde homogène E soumis 
à la gravitation et à une force perturbatrice dérivant d'un poten- 
tiel — ([3y^ -|- y^^) est que, à la surface 

soit une constante. 

Supposons que Telllpsoïde se déforme mais peu, de façon que 
la nouvelle surface 2 ne diffère pas beaucoup de l'ellipsoïde, et 
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cherchons à quelle condition S sera une figure d'équilibre, en 
continuant à supposer que le fluide reste soumis aux mêmes 
forces, c'est-à-dire à la gravitation et à la force ajr, ^y, yi:. 

La condition d'équilibre est que, en désignant par V + f» le 
potentiel dû à la figure S ; 

V + ^+— (Pi/' + T^^) 

soit une constante à la surface. 

Soit J^ la longueur de la normale comptée de E à S ; on pourra 
développer X^ en une somme de fonctions sphériques : 



i;=-^?mN; 



et l'on aura à la surface 

4«p 



=1: 



R^S^MN. 



2/1+1 

D'autre part, on aura, a des quantités près de l'ordre de Ç^ 
et, en appelant^ l'intensité de la pesanteur à la surface, il vient 

m 

d'où 

u=u„-g'j;. 

La condition d'équilibre est que U soit une constante et comme 
la condition d'équilibre est que 

^ 4- xrs 
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soit une constante à la surface c'est-à-dire que 

Je dis que gl est une constante sur l'ellipsoïde ; en effet, 
déplaçons Tellipsoïde parallèlement à Taxe Ox^ l'équilibre ne 
sera pas troublé, car par suite de la condition a = o, le travail 
de la force (olt, [3^, y^) est nul. Nous avons vu que Ton avait 
pour ÎJ la valeur 

K étant une constante ; Téquation précédente se réduit en ce 
cas a 

-^ r;s;m,n, _^;km,n, = c^ ; 

ceci ne sera constant que si 

gl = -^ RÎSî. 

L'équation (i) devient alors 

et l'équation d'équilibre devient, en supprimant l'indice zéro, 

R,S, R.S. 



2W + I 3 



0. 



Si la masse est en équilibre sous l'action de l'attraction et de 
la force centrifuge, elle le sera encore si on lui donne un mouve- 
ment de rotation. Cette rotation est une des transformations que 
l'on peut faire subir a un ellipsoïde de Jacobi, sans troubler 
l'équilibre. 

Nous avons vu que dans ce cas, on avait ii un facteur constant 
près (p. i44) 

J:-m,n,; 
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on en conclut que relHpsoïde de Jacobi doit satisfaire à la relation 



R»s, R.S, 



ô ô 

Réciproquement, supposons qu'un ellipsoïde soit en équilibre 
sous Taction de l'attraction et d'une force dérivant du potentiel 

je dis que si on peut le faire tourner d'un angle w sans troubler 
1 équilibre, c'est que P = y. 

En effet, à la surface on doit avoir, en appelant V le potentiel 
newtonien 

V + -f (?/ + T^^) = V, = C''; 

au sommet du petit axe V^V, ; si on fait tourner l'ellipsoïde 
d'un an^le (ù autour de O.r, le sommet du petit axe restera sur la 
surface de l'ellipsoïde, la constante V^ n'aura pas changé, V n'aura 
pas changé non plus. 
Or 

2 V + ?^ijot/ + yzo.z = oV, ; 
il faudra donc que sur la surface 

et en tenant compte des relations 

Oî/ = tôZ 

on devra avoir sur la surface 
donc 
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Discussion de l'équation des ellipsoïdes de Jacobi. — Nous 
avons donc a discuter l'équation des ellipsoïdes de Jacobi 



(,^ R.s, _ R.S 



l"^! 



5 3 

Cette équation ne convient pas aux ellipsoïdes de révolution, 
cela provient de la façon dont on a établi cette équation : dans le 
cas d'un ellipsoïde de révolution, on aurait en effet t^=o. 

Au lieu de discuter l'équation (ij,nous allons discuter l'équa- 
tion plus générale 

1< = ■ . — = o 

où l'on a p" > a^ ; m étant l'ordre de la fonction R/., n étant 
l'ordre de la fonction R,. 

F 

Fet-TTy^ sont de même signe, de plus R;^ ne peut pas s'annuler 

puisque les racines des polynômes R sont comprises entre a^ 
et c^ ; donc l'équation 

p _ F S, S, R? _ 

'- R-; (2/« + i)R, (2«+i)R, R| 

a les mêmes racines que F=:o. 

F 
Les racines de l'équation -j^j- = o sont séparées par celles 

Rit 
de la dérivée : 

^F, I d S, I RJ d S, 



du 2m-\-i du R/. 2/i+i R| du R, 

I S, rf / R; 



2/1+ I R, du 



m 



On a par définition 



d S;^ 2m-\-i 



du R, m ' 
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et par suite 

dv\ ^ 3 S, R, r;r,-r,r[ 

du 2«+ 1 R, ■ R,, Rf 

S..S, r;.r,.-r,r;. . 



u 



2/2+1 R 

(IV 
mais S. ne peut s'annuler : donc les racines de — j — - soiil les 

du 

racines de RjU^t — R,Hfc. Nous ne nous occupons évidemment 

que des racines supérieures à cr. 

Cherchons la dérivée de cette fonction par rapport à u, c'est : 

R;.'R,-R,R1'; 
et, d'après la définition de R, et R,., 

R;.' = (np^+K)R,, 

on a 

R,R, )[n (rt + I ) — m {m + , )] p^ + K, — K,) . 

Les deux premiers facteurs n'ont pas de racines en p" supé- 
rieures il a"^; le troisième ne peut évidemment en avoir qu'une 

dF 
seule : donc — - — ne peut en avoir plus de deux, F^ ou F ne peu- 
vent en avoir plus de trois. 

Supposons que -Tj^ soit constamment croissant, sa dérivée ne 
s'annulera pas ; donc Fi sera constamment décroissant et ne 

pourra avoir qu'une seule racine. Au contraire, si ■rr- est toujours 

Ra 

décroissant, Fi sera toujours croissant. D'ailleurs pour // = o, 

p =oc , -^ est à peu près égal à— ; donc F part de zéro et est 

toujours croissant ou toujours décroissant, donc n'a pas de 
racines, sauf la racine p = cx) qui n'est pas admissible. 

Dans le cas le plus général, F peut avoir, nous l'avons vu, 
trois racines : mais parmi celles-là se trouve la racine p =00 , 
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F a donc zéro, une ou deux racines supérieures à a^ \ s^"U^ ^^* 

toujours croissant ou toujours décroissant, elle n'en a aucune. 
Revenons alors à Téquation 

R^ R.S, 

O 2/1+1 

(R\ ^ 
-|~ j sera un polynôme 

dont toutes les racines sont réelles et comprises entre a^ et 
c^\ donc la dérivée ne peut s'annuler pour une valeur de p^ supé- 
rieure à a^. L'équation F n'a pas de racine supérieure à d^. Sup- 
posons maintenant que R, ne soit pas divisible par y p''* — a^ , et 
que n soit plus grand que i, je dis que l'équation a au moins une 
racine. En effet, substituons a^, le premier terme s'annule, le 
second est négatif : on a donc une somme négative ; remplaçons 
p par 00 ; le premier terme de l'équation est a peu près égal à 



p V ^ 2/1+1 / 



P 

le résultat est positif; il y a donc un nombre impair de racines 
entre a^ et + oo , il ne peut y en avoir plus de deux, il y en a donc 
une et une seule. 

Ce raisonnement n'est pas valable si /i= i, c'est-à-dire si 

ou bien 

dans ces deux cas, la quantité 

Rf _ ç^ — P 

et les quantités analogues sont toujours croissantes lorsque p^ 
varie de oo a a^. Donc l'équation correspondante n'a pas de 
racines : c'est un théorème déjà énoncé (p. i45). 
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Appliquons ce principe a Tellipsoïde de Jacobi. Il faut consi- 
dérer l'équation 



K ' 



3 5 

R^ n'est pas divisible par R^ : donc il y a une racine, et une seule. 
Donc parmi les ellipsoïdes homofocaux, il y a un ellipsoïde 
et un seul, qui soit une figure d'équilibre ; c'est l'ellipsoïde de 
Jacobi. 

Figures dérivées de l'ellipsoïde de Mac Laurin. — Suppo- 
sons que b^=:c^ et cherchons s'il y a des figures peu diDFérentes 
de l'ellipsoïde de Mac Laurin. Pour qu'il existe une figure d'équi- 
libre S peu difierente de cet ellipsoïde, il faut que l'on ait la 
condition 
, . R,S, RfS, 



3 2/1 -f-i 

Nous avons vu qu'en posant 



on avait 

Nj = cos jy'^, 



M,=AF 



(\/:-^)- 



h étant une constante. On a posé 

F (Ç)= (i — Ç')'^ D» + ^ (i — \'f, 

La condition nécessaire et suffisante pour que l'équation (i) 
ait une racine est que R,- ne soit pas divisible par y/p** — d^\ donc 
M; ne doit pas être divisible par y/a^ — |jl", donc n-^j) doit être 
pair, n etp étant déterminés, une figure d'équilibre voisine d'un 
ellipsoïde de révolution est déterminée par la relation : 

J^ r^ e /M,- cos p'f 

ou, d'une façon plus générale, 

t^ = rSli (e cos /r^ + e' sin pz>) . 
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On peut supposer £' = o : cela revient à faire tourner les axes 
de coordonnées d'un angle convenable autour de 0.^:. La figure 
que Ton obtient ainsi est une figure admettant Taxe Ox comme 
axe de symétrie d'ordre/^. F (;) ne change pas si Ton change ^ en 
— ^, donc le plan des yz est un plan de symétrie. S'ip = o, 2J ne 
dépend pas de cp, la figure est de révolution. 

On ne doit pas conserver parmi les équations obtenues celle 
oii n = Oy car le volume ne serait pas conservé; de même on doit 

exclure /i = i, car alors R,- serait divisible par y p^ — a^ : donc 
la valeur minimum de n est 2; à cette valeur de n correspondent 
deux valeurs dey^ : y^ = o, ety^=2. 

Quelles sont celles de ces figures qui sont ellipsoïdales : con- 
sidérons un ellipsoïde E^, et un autre E, qui en soit peu différent 
mais non liomofocal. Le potentiel de chacun d'eux à son inté- 
rieur est une fonction du second decrré des coordonnées carte- 
rs 

siennes ; la différence des potentiels est donc aussi un polynôme 
du second degré : c'est d'ailleurs le potentiel dû a la couche com- 
prise entre les deux ellipsoïdes. On peut par suite développer ce 
potentiel sous la forme 



VeRMN, 



R, M, N étant des fonctions sphériques d'ordre o, i ou 2. 
A la surface on a 



V=2 £R„MN, 



et a l'extérieur on a par conséquent 



-S-* 



SMN, 



l'épaisseur de la couche produisant ce potentiel est donc 

171 -\- i t 



^=s 



4« S, 



MN 
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Il n'y a que des termes d'ordre o, i ou 2, mais nous avons dit 
que M ne pouvait être d'ordre o ou 1 ; il reste donc n = i et 
p = o ovLp = 2. Siy; = o, la surface considérée est de révolution, 
c'est un ellipsoïde peudilFérent du premier; siy^ =2, on a un ellip- 
soïde de Jacobi. 

Variations de w^ avec l'aplatissement. — Si on fait varier 
p" de a^ à + 00, l'aplatissement des ellipsoïdes considérés varie 
de o cH I . 

Lorsque l'aplatissement est nul w" = o ; l'aplatissement crois- 
sant, (0" va aussi en croissant; mais il peut avoir des maximums et 
des minimums; lorsque l'aplatissement augmente, le moment 
d'inertie principal augmente aussi. 

On a vu que 



0/ 



TU, = II=\V + -::^ J 

2 



ne pouvait pas dépasser une certaine limite. On a de plus 



(0^ 



^11 = ^W + (ùJdi,) H dî: 

2 

dang une position d'équilibre on a 

^W + -^ rfJ=o; 

2 

il en résulte que pour une position d'équilibre on a aussi 

dll = (O îdtù. 

Lorsque l'aplatissement augmente, J augmente indéfiniment; 
donc Cl) doit diminuer et tendre vers zéro ;(0" part de zéro, revient 
a zéro, donc a un maximum dans l'intervalle. 

Remarque. — Faisons varier p- de 00 à a^. On obtient des 
ellipsoïdes de plus en plus aplatis : on en rencontrera pour les- 
quels on aura 

R,S. ^ R,S, 

3 2/2 + l ^ 
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je dis que la première figure d'équilibre que Ton rencontrera 
est celle d'un ellipsoïde de Jacobi. En effet, si on suppose que -^ 
est toujours croissant dans Tintervalle considéré, on voit que 

est toujours positil. 



5 271+1 

Pour p'^ = 30 on a 



R,S, ^ R,S, R,S, 



>— ^-^> 



O ô 2/4 + 1 

si p décroît, on finira par avoir 

RiSj R4S4 R,Sj 



3 5 2/1+ I ' 

p continuant a décroître, on aura alors 

R,S. ^ R,S, 

i 2/1+ I ' 

mais on voit que Ton rencontre d'abord un ellipsoïde de Jacobi. 

Figures dérivées de rellipsoïde de Jacobi. — Il reste a étu- 
dier les figures d'équilibre qui différent peu de l'ellipsoïde de 
Jacobi. Pour cet ellipsoïde on a 

3 ~ 5 ' 
pour une figure peu dififérente, on aura 

R.S, R.Si 



i 2« + I ' 

la hauteur de la surface au-dessus de l'ellipsoïde étant 

On a donc à considérer les deux équations 

R,S, ^ R,S, ^ R,S, 
3 5 2/1 + 1 * 
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On a trois înconnues p^ — a^, p^ — P, p^ — c^ ; et on a pour le 
volume qui est donné 

Pour que Téquation 

R,S. ^ R,S, 

3 2/1 + 1 

ait une racine, il ne faut pas que R,- soit divisible parRj=y p^ — d^; 
pour que l'équation 

R,S, ^ R,S, 

5 2/1 + 1 

Rf 
ait une racine il faudra que -^^ ne soit pas constamment croissant, 

Rf a donc comme seules formes possibles les formes suivantes, où 
P désigne un polynôme en p^. 

p, 

Pxy/ p' — 6% 
Pxy/ p^— 6% 

Pxv/{p' — ^')(p' — ^■'). 

Voyons si ces quatre formes sont possibles et rappelons que 
l'on a R, =v^(p^_^^)(p'_c'^). 

i" Type, Rj est un polynôme ; si a est la plus grande racine 
de R, on pourra poser 

R, = (p^-a) n,. 

a est compris entre «^ et c^ et on a 

R, V p^-/,^ • VT^r^ " 

n,. croît constamment, 






0^ — a A 

' croit constamment, 



p^— c^ 



J- 1- croît si a est inférieur a h^^ décroît dans le cas con- 

p^— ^^ ' 
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traire ; il ne faut pas que i / ~^^ soit décroissant : donc a < Z>- 

et toutes les racines de R, sont par suite comprises entre b' et c~. 
Parmi les fonctions de Lamé d'ordre /i, il y en a une et une 
seule ayant toutes ses racines comprises entre b~ et c^\ c'est celle- 
là que Ton doit prendre pour R,-. 

2® Type, Prenons R, = \/p^ — b^ P et appelons a la plus 
grande racine du polynôme P ; on a 



^=v.^\/^"- 



Les trois facteurs sont constamment croissants, donc toutes les 
fonctions de ce type sont h rejeter. 
3° Type. Posons 

R,=v7ir^(?'^-«)n. 

n étant un polynôme dont les racines sont inférieures à a; nous 
avons 



R. 



s/7^\/^^^' 



Pour que Péquation puisse avoir une solution, il est néces- 
saire que a< i^. Donc toutes les racines de Psont comprises entre 
b^ et 6'\ 

4® Type, On voit de môme que les fonctions du quatrième type 

sont telles que -~ soit un polynôme sans racines supérieures à 

a"^, donc il n'y a pas à considérer de fonctions du 4" type. 

Donc pour qu'une figure peu différente d'un ellipsoïde de Jacobi 
soit une figure d'équilibre, il est nécessaire que R, soit une des 
fonctions que nous avons appelée Ro„. Je dis que cette condition 
nécessaire est aussi suffisante : en effet, supposons que i" varie 
de c^ à a^; à chaque valeur de b' correspond une valeur de p" telle 
que 

3 5' 
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considérons alors hi fonction, 



5 271 -f- I ' 

pour P = 6'^, -5-^ est constamment croissant, donc F est positif; 
si P=za^^ la valeur correspondante de p^ tend vers a^, alors 

R, = o, 

R^.i^ o puisque R,- ne contient niy/p" — à' ni y/p^ — c^^ et que toutes 
ses racines sont comprises entre b^ et c- : donc 

F<o 

et il y a certainement un système de racines pour les équations 
considérées. 

Nous voyons donc qu'il y a des figures d'équilibre voisines des 
ellipsoïdes de Jacobi. 

Remarque. — Quelle est la première figure d'é([uilibre que 
Ton rencontre ? c'est-à-dire, pour quelle valeur de n la valeur 
correspondante de p est-elle maximum? 

Considérons deux fonctions Ro„j et Ry„, 

ju > Ji ; 
je dis que 



Hon 11 



est constamment croissant; nous allons montrer que la dérivée, 
prise par rapport à n^ est constamment positive; cette dérivée 
est à un facteur positif près 

R'/,R,-R,.R'., 

si R, n'est pas divisible par yp- — d^t R/ Test; et réciproquement. 

Mais ici R, et R,^. ne sont divisibles ni pary c^ — cr ni par y p" — 6"; 
il en résulte que 

R',R,-R,R': 

s'annule pour p^ = cr. 
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Si on exprime R, et R^t en fonction de l'argument //, on a 

R, = F(/»«), 

RÎ.=://«F'[/,(«)], 

et on voit que la dérivée n'a pas de racine supérieure à a^. En 
effet la dérivée de cette fonction (prise par rapport à u) est 

(Ri' R,-R,Ri') = (?* - 'f,)R.R* 

comme on Ta déjà vu; R»- et R/r. ne peuvent s'annuler, «p,- — '^^ ne 
peut s'annuler deux fois, puisque c'est un polynôme du premier 
degré en p^. La dérivée ne peut donc s'annuler plus de deux fois 
pour a^ ^ p^, c'est-a-dire pour o < u ^e^, 

La dérivée qui s'annule pour u = o, u = ^,, ne peut donc s'an- 
nuler dans l'intervalle; ■—- est donc constamment croissant lorsque 

Rf 

p*^ varie de a^ à -f- oo . Il en résulte que l'on a aussi 

R,S, ^ R,S, 



2/i+ I 2/72+1 ' 

lorsque p^ varie de oo à d^^ on rencontre donc d'abord une 
valeur satisfaisant a l'équation 

R,S, ^ R.S, 

5 271 + I ' 

puis ensuite la valeur de p^ satisfaisant à l'équation 

R,S, ^ R,S, 

5 2/2 + 1 

La première des valeurs de n que l'on trouve est donc /i = i , 
puis /i = 2 : on a déjà examiné ces valeurs; il faut maintenant 
examiner la surface voisine de Tellipsoïde de Jacobi et corres- 
pondant h Rç 3 ; nous allons étudier d'une façon générale la 
surface correspondant à R(,„. 

L'épaisseur de la couche sera 

i; = £/M N, ; 
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cette épaisseur sera nulle aux points où Ton a 

M,N,= o, 



i6i 



ou bien 



R.M.N. = o. 



Ceci est décomposable en facteurs quadratiques de la forme 



X' 






y 



(X^—b' 



CL 



I • 

2 Ji * > 





Fig. 19. 



les a sont d'ailleurs compris entre P et c^ ; les surfaces corres- 
pondantes sont donc des hyperboloïdes à deux nappes coupant 
l'ellipsoïde suivant des 
lignes de courbures 
d'un seul système. 

D'un côté d^une de 
ces lignes Ç sera positif; 
et de l'autre côté ^ sera 
négatif. On peut, avec 
ces données construire 
la figure dans l'hypo- 
thèse de /i = 3. On re- 
présente (fig. 19) en 

ponctué le contour apparent de l'ellipsoïde, en traits pleins l'in- 
tersection de la surface avec l'ellipsoïde ainsi que le contour 
apparent extérieur à l'ellipsoïde, en pointillé le contour apparent 
de la surface intérieure h l'ellipsoïde. L'une des lignes de cour- 
bure est d'ailleurs l'ellipse principale, comme on le voit en re- 
marquant que 

R0.3 = \/7^=^x P- 

On construira avec les mêmes conventions la figure d'équilibre 

correspondant à R,,,^ (fig. 20). 

Je remarque qu'on ne peut prendre pour les fonctions R,- de 
fonction du second ordre, car il faudrait prendre pour la fonction 
Rq2 un polynôme de la forme p^ — A^, ce qui donnerait un autre 
ellipsoïde et nous avons vu qu'il ne pouvait y avoir qu'un ellip- 
soïde de Jacobi pour une vitesse de rotation donnée. 

PoiNCARÉ. — Figures d'équilibre. 11 
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-C'-'' 

<.^* 




-----. 





Fi g. 20. 



En résumé, si on fait varier h^ de c^ à a*, on obtient des ellip- 
soïdes à chacun desquels correspond une valeur co et une seule ; 
(0 ne peut avoir de maximum ou de minimum, car s'il y en avait 

un cOp à une valeur voi- 
sine de cOj correspon- 
draient deux figures 
ellipsoïdales ce qui n'est 
pas. Donc co varie tou- 
jours dans le même 
sens. 

Supposons que b^ et 
6*^ tendent vers zéro, la 
forme de l'ellipsoïde 
tend alors vers une ai- 
guille allongée; co a d'ailleurs la même valeur pour deux ellip- 
soïdes semblables. Il suffit donc de voir ce qui se passe. dans le 
cas où 6 et c restant fixes, a augmenterait indéfiniment. On 
aurait alors un cylindre elliptique, la force résultante de l'at- 
traction du cylindre sur lui-même et de la force centrifuge 
doit être normale au cylindre, c'est-à-dire parallèle au plan yOz, 
Considérons un point d'une génératrice passant par l'axe de 
révolution : par raison de symétrie, l'attraction doit être parallèle 
au plan yOz^ il en résulte que la force centrifuge qui est per- 
pendiculaire a Oz ne peut être que nulle, et par suite co tend 
vers zéro lorsque a augmente indéfiniment. 

On peut d'ailleurs remarquer que le cylindre ne peut être 
que de révolution pour que l'équilibre soit possible. 



Stabilité des figures trouvées 



■\ 



Représentation graphique des résultats précédents, — 

Prenons deux axes de coordonnées (fig. 21) et portons en abs- 

R? R^ 

cisses les valeurs de -p|- et en ordonnées -p|^; ces quantités sont 

les rapports de l'axe moyen et du grand axe au petit. 

On a un premier point A correspondant au cas où ces rapports 
sont égaux à i, la figure correspondante est une sphère. La partie 



3/ 


! Vl' /^ 

! \ Vf 

! V' yé 




• • 


O 


JC 



Fig. 21, 
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de la bissectrice de l'angle xOy partant du point A correspond 
aux ellipsoïdes de révolution. Les ellipsoïdes de Jacobi sont 
représentés par une courbe CBC symétrique par rapport à cette 
bissectrice et la coupant en 
un point B. 

Pour les ellipsoïdes de 
Mac Laurin, co* va en crois- 
sant quand on part en crois- 
sant du point A, jusqu'à ce 
que l'on arrive en un point D, 
puis il décroît quand on s'é- 
loigne du point D jusqu'à 
l'infini où il arrive à la 
valeur zéro. Sur la courbe 
CBC deux points symétri- 
ques par rapport à AH re- 
présentent le même ellip- 
soïde ayant tourné de 90** : co^ décroit quand on part de B 
pour aller vers C ou vers C; les points C et C s'éloignant indé- 
finiment correspondent au cas du cylindre de révolution, la 
courbe CBC est asymptote aux droites a: = i , y = i. 

Nous marquons sur AH les points E, F, correspondant aux 
figures peu différentes de l'ellipsoïde de révolution, de même sur 
CBC, nous remarquons les points GKLM, G'K'UM', correspon- 
dant aux valeurs R^,,, donnant les figures peu différentes d'un 
ellipsoïde de Jacobi, que nous avons rencontrées précédemment. 

11 faut remarquer que les paramètres que nous avons pris ne 
correspondent pas aux vraies données du problème. Les para- 
mètres sont le volume T, la rotation *o et le moment principal 
d'inertie J. Si on pose [jl=(o^J, [jl et T sont des données du 
problème, puisque w^J doit être constant; [jl est du cinquième 
ordre par rapport aux longueurs, T est du troisième ; si on 
pose 

T 3 



1 » 



M sera une donnée du problème indépendante de l'unité de lon- 
gueur choisie. 
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Liouville a cherché à quelles valeurs de M correspondaient les 
divers points des lignes AH et CBC. Il a montré que lorsqu'on 
se déplace sur AH, M croît lorsqu'on va de A à D et décroît 
si on va de D vers H. Sur la ligne CBC, M décroît lorsqu'on 
s'éloigne du point B soit vers C, soit vers C. 

Donc si iù est très grand, il n'y a pas de figure d'équilibre 
possible; si co est compris entre cOg et co^, il y a deux figures pos- 
sibles qui sont des ellipsoïdes de révolution. 

Si (o est inférieur à co^, il y a quatre figures d'équilibre pos- 
sible, dont deux sont des ellipsoïdes de Mac Laurin, et deux des 
ellipsoïdes de Jacobi ; ces deux ellipsoïdes de Jacobi sont d'ailleurs 
égaux comme on l'a déjà remarqué. 

Courbes d'équilibre. — Considérons un système dépen- 
dant de n variables x^x.^ ... x„, soumis à un système de forces 
dépendant d'un potentiel F [x^x^ .,, x^t )^) ; la condition néces- 
saire et suffisante de l'équilibre est que F soit maximum par 
rapport aux jp, X étant donné ; on a ici 



2 2J 

La condition d'équilibre est que 

dF ÔF ÔF 



o; 



àx^ dx^ (^Xn 

pour que l'équilibre soit stable, il faut de plus que la forme 

soit une forme définie négative ; on peut la mettre sous forme 
d'une somme de n carrés : 



a> 



= 2j *^'^*'- 



Ces coefficients a^^. sont appelés coefficients de stabilité ; si l'un 
de ces coefficients s'annule, le discriminant delà forme <^ s'an- 
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nule, puisque cette forme est alors réductible à une somme de 
moins de n carrés ; ce déterminant est, en posant 



A = 



^Xj^Xk 


«1,^. 


^1,1 ^1,2 • • 


• • «l,n 


«2.1 «2.2 

• 
• 


«2,n 


• 
• 

«n,2 «n,2 


(^n,n 



Si deux coefficients clj^ sont nuls, tous les mineurs du premier 
ordre de A sont nuls ; nous supposerons qu'il n'en est pas ainsi. 
Si A est diflférent de zéro, on peut résoudre les équations 



ÔF 
^x. 



ôF 



o 



' d:r„ 



ôF 
= o, — , -^ — = 



dx 



n 



et on aura 



ôF 
Différentions la relation -^ — = o par rapport à )., on aura 



s 



^x^ 

d^F dxi ô^F 

+ 



hx^^Xi d\ ' dx^ô). 



ce qui peut s'écrire 



s 



a 



dxi 



Ô^F 



dx, ÔÀ 



Si A est diflférent de zéro, on peut trouver les valeurs de ,>' , 

puisque l'on a n équations linéaires à n inconnues et les x seront 
des fonctions uniformes de X, dans le voisinage des valeurs con- 
sidérées. 

Supposons A nul, sans que tous les mineurs du premier ordre 
le soient, supposons que l'on ait par exemple 



ÔA 



da 



^o, 



l,t 
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Nous pouvons alors résoudre les n — i dernières équations 
linéaires par rapport à 

dx^ dx^ dXfi 



d/x ' dX ' ' d\ 

ce qui montre que x^, x^^ ... Xn seront fonctions uniformes de x^ 
et de X. 

Substituons dans la fonction F les valeurs de x^^x^^ ... x^^ en 
fonction de x^ et X tirées des équations 

ÔF ôF 



= 0,...,-^ — = o : 



d^2 àxn 

on aura une fonction ^(x^^X), et on a 

_rf^_^ ôF_^ ôF dxn 

dx^ dx^ àx^ dx^ "* àxn dx^ 



mais par hypothèse 



ôF ôF 

rO,..,,-^ = O 



0^2 ^Xn 



quand on y remplace x^x^.,,Xn en fonction de x^ et de X ; on a 
donc 

di^ _ ôF 



dx^ hx^ 



La condition d'équilibre est alors que Ton ait 

4 _ 



dx^ 



C'est une relation entre x^ et \ qui définit une courbe, x^ sera 

une fonction uniforme de X si — r-^ est différent de zéro ; la 

dxi 

condition A = o équivaut donc à , ^ =o. 

ClX^ 

Supposons que cette condition soit remplie sans que l'on ait 

d^i^ 
en même temps --7 — ^=- = 0^ la courbe a alors une tangente ver- 

1 

ticale (fig. 22). 
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Quand X traverse la valeur de l'abscisse de cette tangente, il y a 
pour x^ deux valeurs réelles qui viennent se confondre, puis 
deviennent imaginaires ; de même 
pour x^y ... x^, 

—7 — ^ÏT^^^' point corres- 
pondant de la courbe est un point 
double, à tangentes distinctes ou 
confondues ; Tune des tangentes 
peut être verticale ou non, on a 
deux dispositions principales : 
celle qui correspond à deux tan- 
gentes distinctes et non verticales (fig. 23) et celle qui 
correspond à une tangente verticale et une non verticale 
(fig. 24). 




Fig. 22. 




wC 




Fig. a3 



Fig. 24. 



On a deux cas à distinguer : dans le premier cas, deux valeurs 
réelles de x se confondent et redeviennent réelles quand on tra- 
verse la valeur correspondant au point double. 

Dans le second cas, deux groupes de valeurs réelles se confon- 
dent et deviennent ensuite imaginaires lorque X traverse la valeur 
critique dans un sens convenable, une autre valeur de x est 
réelle avant et après la valeur critique. Les autres dispositions de 
figures peuvent se ramener à celles-là, du moins au point de vue 
qui nous occupe. 

Echange des stabilités. — H y a équilibre stable, avons- 
nous dit, lorsque F est maximum, x^x^ ... x^ étant variables. Il 
faut a fortiori que F soit maximum si x^ étant donné on fait 
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varier Xg, Xg, ...,-^n d'une façon arbitraire, il faudra donc que 

ôF ôF ôF 



ÔJC, 



ÔJC 



3 



ÔOTn 



O. 



On tirera oc^^x^^ •j-^n de ces équations dont le déterminant 
fonctionnel n'est pas nul par hypothèse, on portera ces valeurs 
dans F [x^x^ .-.^n)? c* on aura une fonction qui devra être maxi- 
mum pour qu'il y ait équilibre. Donc ^ {x^y X) doit être maximum; 

il y aura donc équilibre si -7-^ = o, mais il n'y aura équilibre 
stable que si l'on a en même temps i\ <o. 

UX \ 

Un arc de la courbe -7-^ = partage la région du plan où 
il se trouve en deux parties, 1 une ou jl > o, 1 autre ou ^J <o ; 



dx^ 



dx. 



si c'est un arc sans singularité, , \ ne sera pas nul et sera 

(IX A 

négatif si la région où , est positif est située au-dessous de 

l'arc de courbe ; et au contraire sera positif, si ^ est positif au- 
dessus de l'arc de courbe considéré. Le premier cas correspond à 
une position d'équilibre instable. 

En couvrant de hachures la région où /^ est positif, on voit 

que la figure 25 correspond au cas d'un équilibre stable et la 
figure 26 au cas d'un équilibre instable. 



3^ 








y 



JCy 








oD 



Fi g. 25. 



Fig. 26. 



Dans le cas où l'arc BAC a une tangente verticale en A, deux 
cas peuvent se présenter, suivant que c'est la partie convexe 
ou la partie concave qui est hachurée. 
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Dans la figure 27 : 

l'arc BA correspond à des positions d'équilibre instable, 
l'arc AC » » stable ; 

dans la figure â8 : 

l'arc BA correspond a des positions d'équilibre stable 
l'arc AC » » instable. 





Fig. 17. 



Fig. 18. 



Ce qu'il y a d'essentiel, c'est qu'en suivant l'arc et en passant 
par le point A, les deux systèmes de figures d'équilibre échangent 
leur stabilité. 

Dans le cas d'un point double sans tangente verticale, deux 
cas peuvent aussi se présenter; dans les figures 29 et 3o, les 
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Fig. ag. 

arcs BA et AC correspondent à des figures d'équilibre instable, 
les arcs D'A et AC correspondent à des figures d'équilibre stable. 
Remarquons encore qu'en suivant l'arc BAC ou l'arc B'AC les 
figures d'équilibre, de stables deviennent instables ou récipro- 
quement. 

Enfin dans le cas d'un point double avec une tangente verti- 
cale, on arriverait à la même conclusion, les figures d'équilibre 
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que l'on obtient en suivant un arc quelconque échangeant lear 
stabilité quand on passe par un point où ^;4-= o (fig. 3i et Sa). 




Fig. 3ï. 



St&bilité de l'équilibre des Bgures trouvées. — Appliquons 
^eci à notre problème et considérons la figure 21 ; ce que nous 
ivons appelé y. représentera ici ce que nous avons appelé ). et 
correspondra à x^. 

La fonction F, qu'il faudra rendre maximum est 



W- 



aJ 



M^ augmente lorsqu'on se déplace de A sur D, et décroit de 
D à H; sur AH, M va en croissant, et croit aussi de 6 à C ou 
à C. 

Il s'agit maintenant de traiter la question de stabilité. 

Parlons de A, où w ^ et où on a une sphère, nous savons 
que dans ce cas l'équilibre est stable : donc tous les coefiîcients 
de stabilité sont négatifs. Ils cesseront d'être tous négatifs quand 
on rencontrera une figure de bifurcation. 

La première bifurcation que l'on rencontrera en suivant AH 
sera au point B, où l'ellipsoïde de Mac Laurin diffère peu d'un 
ellipsoïde de Jacobi. Le coefficient de stabilité s'annulera donc 
en même temps que 

R. S. R, S, 
3 5 • 

Après cela il ne s'annulera plus ; quand on arrivera en E, un 
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autre coefficient s'annulera en changeant de signe, mais le 
premier coefficient restera positif. Donc de A à B nous avons un 
équilibre stable, au delà de B il est instable. 

Sur la branche BC, on a au contraire des ellipsoïdes qui sont 
stables jusqu'en G ; en G l'un des coefficients s'annule et au delà 
reste positif, donc à partir de G sur la ligne BC on aura des 
positions d'équilibre instable, mais à partir de G commenceront 
de nouvelles figures d'équilibre et ce sont celles dont il faudra 
étudier la stabilité ; il y a une difficulté, car, pour ces figures, 
M peut croître ou décroître. 

Il est aisé de voir que l'on a en tous cas un maximum ou un 
minimum ; on a en effet pour ces figures, si elles correspondent 
h des formes peu difTérentes de celle que l'on a en G 

en G même, e = o. Si on change e de signe, la figure a une 
forme symétrique de la première, M garde donc sa valeur ; il en 
résulte qu'en G la quantité M passe par un maximum ou un 
minimum. Si M<MoOn a après une seule figure d'équilibre; si, 
au contraire, M > M^ on aurait trois figures d'équilibre dont une est 
stable. On n'a pas encore fait le calcul complet, bien qu'il semble 
plus probable que l'on a un minimum. 

Conclusion, — Quoi qu'il en soit voyons les deux hypothèses : 
supposons une masse fluide homogène animée d'un mouvement 
de rotation et soumise à un refroidissement assez lent pour 
qu'il ait le temps de se propager dans la masse ; supposons de 
plus que le mouvement de rotation soit le même dans toute la 
masse, et cela en supposant le frottement suffisamment grand. 

La vitesse du mouvement de rotation n'est pas constante, ce 
qui est constant c'est iji = io'^J, en vertu du théorème des aires ; le 
volume T doit diminuer puisque la masse se refroidit, donc M 
augmente. On a des figures d'équilibre stable jusqu'en B ; après 
cela, ce qui devient stable ce sont les ellipsoïdes de Jacobi et cela 
jusqu'en G ; le refroidissement continuant, on ne peut plus avoir 
d'ellipsoïdes de Jacobi, puisqu'ils sont instables. Si, pour l'autre 
série, M^ est un maximum, cette seconde série ne donnera aucune 
figure, et le corps se dissipera dans l'espace; si, au contraire. 
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Mo est un minimum, on aura une sorte d'ovale étranglée que 
l'étranglement divise en deux parties inégales. 

On n'a pas fait le calcul jusqu'au bout pour voir ce qui 
arriverait ensuite : peut-être l'étranglement s'accentuera-t-il et 
la masse pourra se diviser en deux parties. Remarquons que ceci 
ne peut s'appliquer directement à l'hypothèse de Laplace, car 
la nébuleuse considérée par Laplace n'est pas homogène, elle 
doit au contraire être fortement condensée au centre. 

Quoi qu'il en soit, on ne sait pas jusqu'à présent si en G corres- 
pond pour M un maximum ou un minimum. 



CHAPITRE VIII 



ANNEAU DE SATURNE 



Etude des ûgures annulaires d'équilibre. — On peut faire 
trois hypothèses sur l'anneau de Saturne : on peut admettre qu'il 
est solide, qu'il est fluide ou qu'il est composé d'une masse de 
solides indépendants ; la troisième hypothèse connue sous le nom 
d'hypothèse de Cassini paraît être l'hypothèse véritable. 

HYPOTHÈSE DE l'aNNEAU SOLIDE 

Historique. — Laplace est le premier qui ait examiné cette 
hypothèse : et il a remarqué que si l'anneau était homogène son 
mouvement serait instable ; pour peu que son centre s'écartât 
du centre de gravité de la planète, Tanneau se précipiterait sur 
celle-ci ; l'attraction de Saturne ne tendrait qu'à accentuer l'ano- 
malie. Si donc l'anneau est solide il doit présenter des irrégula- 
rités ; Laplace ne chercha pas à déterminer l'ordre de grandeur 
de ces irrégularités ; Maxw^ell a fait le calcul et a reconnu 
qu'elles doivent être considérables ; ces irrégularités sont inad- 
missibles. 

Equations du mouvement. — Parmi les forces qui sont appli- 
quées à l'anneau il n'y a pas à s'occuper de l'attraction de l'anneau 
sur lui-même parce que l'anneau étant solide ses diverses parties 
n'éprouvent pas de déplacement relatif; il n'y a donc à s'occuper 
que de l'attraction de Saturne et de celle des autres corps 
célestes ; nous négligerons cette dernière action et nous allons 
nous occuper seulement de l'attraction de Saturne. 

Nous allons supposer que l'anneau a son centre de figure au 
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centre O de la planète : Tobservation montre qu'il en est ainsi; 
le centre de gravité G de Tanneau est supposé différent de O. 
Un point M de l'anneau, puisque l'on ne parle pas de l'attrac- 
tion de l'anneau sur lui-même, est soumis à l'action de l'attraction 
de Saturne ; l'ensemble des forces d'attraction fait donc équi- 
libre à la force d'inertie. Un point M de l'anneau est soumis à 
l'attraction de Saturne et à la force centrifuge : si l'anneau con- 
serve son centre de figure au même point cette force centrifuge 
et cette attraction sont constantes : M tournera alors avec une 
vitesse égale à celle d'un satellite placé à la même distance ; 
c'est la condition nécessaire et suffisante de l'équilibre. Si 
l'anneau est homogène, cette condition n'est pas nécessaire, car 
la vitesse de rotation peut être quelconque : chacun des points de 
l'anneau serait soumis à la même force, et l'équilibre aurait lieu 
quelle que soit la vitesse de rotation. Dans le cas général le 
centre de gravité G décrira donc d'un mouvement uniforme une 
circonférence autour de Saturne. 

Prenons comme unité de longueur le rayon de l'anneau, 
comme unité de masse la masse de l'anneau ; son moment d'inertie 
par rapport au centre de figure sera égal à i . Il reste à déter- 
miner l'unité de temps : on peut la prendre de façon que, en 
désignant par M la masse de Saturne et par fie coefficient de 
la gravitation, on ait /M = i . 

Etudions le mouvement : soient x et y les coordonnées du cen- 
tre de gravité, nous avons 



pds= I pcos5.rfs, 





fds =10 sin 5. ds^ 



en appelant p la densité de l'élément linéaire de l'arc ds de 
l'anneau supposé étroit. 



est la masse de Tanneau. 
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Les équations du mouvement du centre de gravité sont 




pds = I p cos s ds = 




d'où 



d^y 
dt^ 



==:X, 



= y' 



Les intégrales de ces équations sont, en prenant une origine, 
des temps convenable, 

j: = a cos ty 
y = p sin t. 

Le point G décrit donc une ellipse dans le cas général ; si 
c'est une circonférence, elle est décrite avec une vitesse angulaire 
égale h l'unité : c'est la vitesse d'un satellite placé k l'unité de 
distance et décrivant un cercle et dans ce cas, a = p. 

Stabilité du mouvement. — Il reste à voir si ce mouvement 
est stable. Supposant négligeable, par rapport à la masse de 
Saturne, la masse de 
l'anneau, nous pou- 
vons supposer Sa- 
turne fixe. 

Par son centre O, 
faisons passer deux 
axes fixes OX et OY 
(fig. 33) ; par le cen- 
tre de figure de l'an- 
neau, faisons passer 
deux axes ox^ oy in- 
variablement liés à l'anneau, leur position, et par suite celle 
de l'anneau, sera déterminée, si l'on se donne les coordonnées 
X et y de par rapport à xoy, et l'angle 9 de l'axe ox avec OX. 

Soit G le centre de gravité de l'anneau, qu'on peut supposer 




Fig. 33. 
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sur oj: ; on a oG = «, c'est une donnée du problème ; les coor- 
données X et Y du centre de gravité sont 

X = — y sin 9 + (« — x) cos 8, 
Y= ycos9 + (« — j:) sin 8. 

La force vive de Tanneau est égale à la force vive de trans- 

M 
lation (X'^ + Y'^) augmentée de la force vive de rotation J8'^ 

Le moment d'inertie par rapport au point o est donné par la 
relation I = J + a^. La force vive est donc une fonction de x, y, 8 
et de leurs dérivées et on a 

X'^ H- Y'^ = y^ + y ^+ 9'^ [y^ + (a - ^)^] - 29' [x'y H- (a -or) y']. 

L'énergie potentielle est 

U = — /MV, 

V désignant le potentiel de l'attraction de l'anneau au centre de 
gravité, et comme on a 

on a 

U=— V. 

Supposons que la densité de l'anneau soit une fonction déve- 
loppable en série de Fourier. Si di^ est l'angle sous lequel on 
voit du centre de l'anneau un élément linéaire, on aura, en appe- 
lant p la densité au point déterminé par l'angle ^, 

p = po| i + 2acos^H — ^ cos 2^j; + -J- sin 2^+ ... j; 

si ^ et 7) sont les coordonnées d'un point de l'anneau par rap- 
port aux axes mobiles, on a 



çds ^ I , 
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pT,rf5=0, 




p (?' + 71^) ds 




On trouve encore 

►2» 



p(Ç^ — 'r\^)ds= j pcos2Arfi]; = -^5 







In ds = I p sîn 2<J^rf^ = -^. 



Le potentiel en G de l'attraction de Saturne sur l'anneau 
est, comme on le voit en effectuant le calcul, 

^ V = — U 
Les équations de Lagrange donnent 

rf_ /dT\ _ ^ , au _ 

'dt [àq'J àq "*" dry ""^' 

mais comme U ne dépend pas de q' on a 

d d(T — U) d(T — U) 



dt àq' dq 

On peut écrire 



o. 



T U Ar+.B9'+C ,^.^ 

2 

oii Ton a posé pour abréger 

A = 3-\-7/-\-{a — j:)\ 
B = — x'y-i-ty{a—x), 
C = x" + y"; 

PoiNCARÉ. — Figures d'équilibre. n 
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remarquons que T — U ne dépend pas de 9 ; la première des 
équations de Lagrange donne donc : 



r a(T-U) i 

L 09' J °' 



et, comme U ne dépend pas de 0, il reste 



p étant une constante. Il est d'ailleurs facile de calculer cette 
constante, étant données les conditions initiales ; si on déve- 
loppe l'équation on a 

Ae'+B=/>. 

Supposons a:, = 0, y, = o, et par suite af^ et y', très petits, on 
trouve 

et par suite jo = 9'^; donc c'est la vitesse de Tanneau, supposé en 
équilibre : on a, comme on Ta démontré plus haut, 9'^^ = i ; 
donc p=i i l'équation obtenue est Téquation des aires. 

Transformons maintenant les équations de Lagrange : posons 

H=T — U— /?9'. 

Les équations de Lagrange conserveront la même forme, si on 
y remplace T — U par H ; car, si q est l'une des variables, on a 

ÔH _ Ô(T — U) ÔT 59^ 59^ 

dq ~ àq "^ 09' dq ^ àq ' 



et, comme on a 



ÔT 



on conclut : 



ÔH_ Ô(T— U) 
àq àq 
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De même on volt que Ton a 

ÔH __ d(T — U) 

On peut faire disparaître 9' de ces équations, on n'a en effet qu'à 
remplacer 0' par 

A ' 

et comme, ni T ni U ne contiennent Q, H ne sera plus fonction 
que des deux variables x et y, et de leurs dérivées. On trouve 

2A 2 

Supposons maintenant que le centre de l'anneau soit très près 
de Saturne, c'est-à-dire supposons j: et y petits, et voyons si dans 
le mouvement ultérieur x et y resteront petits. 

Les termes de degré zéro qui entrent dans H peuvent être 
négligés, car ils n'interviennent pas dans les équations de 
Lagrange, Les termes du premier degré en x et y doivent dispa- 
raître d'eux-mêmes : en effet, dans les équations de Lagrange ils 
donneraient des termes constants, et ces équations doivent 
admettre la solution j; = o, y = o. Quant aux termes du premier 
degré en x^ et y' ils disparaissent dans les équations de Lagrange. 
Il ne reste donc que les termes du second degré, que nous 
appellerons W, le reste étant négligeable. 

Développons H, en ne conservant que les termes du second 
degré ; on trouve pour W le développement suivant : 



x^^ . i— a^ 



Ex' + 2Cxy-^Fy'-]-^^ + — y"—x'y + (i— ^a') xy' 



en posant 



3 o . a 



E= -: 7.a 



2 



4 ^-^ ' 4' 
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Les équations de Lagrange, comme on l'a remarqué, peuvent 
se réduire à 

d dW ôW 
dt ^x' àx ' 

d ôW dW 

dt ^y' ôy ' • 

c'est-à-dire 

a:" — 2(1 — a^)y==2Ea: + 2 Cy, 
[i — a^) y" -{- ^ [i — a^)x' = ^Cx + iVy, 

Ce sont deux équations linéaires à coefficients constants ; pour 
intégrer on pose 

y=y^ e»"' . 

On aura quatre valeurs possibles de 10, et l'équation générale 
du mouvement sera 

les l et les m étant des constantes quelconques déterminées par les 
conditions initiales. Pour que le mouvement soit stable, quelles 
que soient ces conditions, il faut que les valeurs de co soient à 
partie réelle négative ou nulle. 

Pour déterminer 10, on remplace x et y dans les équations 
différentielles et on a 

{ù^x — 2(1 — a^)(i)y = 2 Ex + 2 Cy , 
(i — a^) 10^^+2(1 — a^) iùx=2Cx-{-2Fy^ 

on peut éliminer x et y entre ces deux équations, et on trouve 

(o*(i — a^) + to^ ( — iH 1 1 

+ -l(9-a^— P^— 24a' + 8a^a)=o. 

L'équation en (o^ doit avoir ses racines réelles : car si elles 
étaient imaginaires, les racines carrées de ces deux racines ne 
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seraient pas des imaginaires pures, et par suite, comme elles 
sont égales et de signes contraires, Tune serait à partie réelle 
positive. 

Or, pour un anneau homogène, où 

a =zaL= P = o, 
Téquation se réduit a 

(0* — lo^H — ^ = o; 

4 

les racines en sont imaginaires : si a, a, ^ sont voisins de 
zéro, elles seront encore imaginaires. Il faut donc que a, a et ^ 
soient grands, il y aurait donc de grandes irrégularités dans 
Tanneau : et sa régularité apparente serait inexplicable. 

Maxwell a fait le calcul dans l'hypothèse où l'anneau serait 
homogène, sauf en un point, où il y aurait une masse supplé- 
mentaire : et il a trouvé que cette niasse devrait être au moins 

les -^ de la masse totale de Tanneau, c'est-à-dire que l'on aurait 

a>o,8. 

M, Radau a examiné le cas où l'épaisseur de l'anneau serait 
variable d'un point à un autre, et il a trouvé que la densité 
devrait varier de 2,7 à o,o4, ce qui ne parait pas probable. 

De plus, la section de l'anneau étant mince, il lui faudrait 
une grande résistance pour résister à l'attraction des satellites : 
Hirn a calculé qu'il lui faudrait une résistance mille fois supé- 
rieure à celle de l'acier : donc Vanneau de Saturne n'est pas 
solide. 



HYPOTHÈSE DE L^ANNEAU LIQUIDE 

Potentiel d'une circonférence homogène. — Cherchons le 
potentiel d'une circonférence homogène, en un point P de l'es- 
pace. Soient PA = a, PB = b les distances minima et maxima 
du point P à la circonférence. La circonférence étant donnée 
de rayon R, « et è déterminent complètement le potentiel. 
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Soit <i) l'angle AOM (6g. 34) , on a 



2 2 




le potentiel 



Vp = 



= M 




Fig. 34. 



27rXPM 



Comme le potentiel ne dépend que de «, de b et de la masse 
de la circonférence, on peut ramener 
le calcul de V au potentiel d'une cir- 
conférence, en un point de son plan : 
on décrit une circonférence de dia- 
mètre [a -f- ^)î et on calcule le po- 
tentiel de cette circonférence en un 
point situé à une distance a d'une des 
extrémités du diamètre (fig. 35). 

Le potentiel en F est le même, 
on a 




Fig. 35. 



V=M 




diù 



27r 



K/ a^ cos^ — 



+ b^ sin^ — 

2 



mais on a aussi 



V=M 




4 
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VoA'cos^^l^+QP'sin^^ 

Q étant le point de la circonférence dont la projection sur le 
diamètre AB est en P, <J; étant l'angle AOM ; on a donc 

«p (a, b) = ï (OA, QP) 

=<p(-^'v/^)- 
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Si donc on forme une série de quantités 





a b 




• • 

• * m 

• 9 


telles que l'on ait 






^n-l + *«-l 




""" 2 




*n V^«n-1 *n-l> 


on aura 





et cela quel que soit l'entier n. On démontre facilement que «„ 
et bn tendent vers une même limite qu'on appelle la moyenne 
arithmético-géométrique de « et de b. 
Soit m cette limite, on a 




27t\ / m^ cos^ — 4" ^^2^ sin^ -^ . y o 



cp (a, b) est une fonction homogène de a et de è, de 
degré — i , c'est-à-dire que Ton a 

<p ().«, \b) = ) - ^ {p (a, è) ; 

on peut donc écrire 

?(«*)=-f?(|-, •) 

et, remplaçant «p (a, i) par cp (a), on a 



<p(a,i)=i.cp(|.) 



Que devient la formule si P est rapproché de la circonférence? 
a est alors une quantité très petite. On sait que l'on a, en dési 
gnant par L le logarithme népérien 

K 

V=2[x,L — +s, 

P 
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p désignant la distance du point P au point P^ le plus voisin du 
cercle, e étant une quantité tendant vers zéro avec p, iji^ et K étant 
des constantes ne dépendant pas de la position du point .P par 
rapport à P^ non plus que de la position du point P,, sur le cercle. 
En reprenant les notations précédentes on a 

et comme 

on a 

M ~ TîR a ' 

b diffère peu de 2R, à Tordre d'approximation adoptée ; on a 
d'ailleurs 

Cette relation permet de calculer la constante K : on a en efiet 
si on fait i = i on a 



f \ ^ / ^V ^ \ 



Il vient donc dans le cas qui nous occupe, en négligeant a 
devant Tunité 

<p(a)=2cp(2v/a); 



mais on a aussi 



ç (a) = L — ; 



a 



donc 



\ « / \ 2v/a / 

LK — La = 2 LK — 2 log 2 — La ; 
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d'où on conclut 






LK 2 La, 




K-4, 


et finalement 






M u \aj 



avec l'approximation donnée. 

On peut chercher une forme plus exacte pour '^ (a). Posons 

'f(«) = :^ log— ri + A,a + AX+...J 
en se servant de la relation 



. 2 / 2\/a \ 



( 



on trouve facilement que l'on a 

Aj = o, B„ = B. = o, 

4' ^' — 1 

de façon que l'on peut prendre pour 'f [a) la fonction 

4 



K = -r, B,=--^L2, 



I •- (i) 



et cela en négligeant seulement les termes en — r^ . 

Masse annulaire de révolution. — Soit G le centre de gra- 
vité de la section méridienne, menons Gy parallèle à Taxe de 
révolution et Gx dans le plan méridien, perpendiculaire à G// 
et tel que le sens positif soit dirigé de Taxe de révolution vers 
le point G, soit l la distance du point G à cet axe de révolution. 
Un élément de surface da' de la section méridienne, de densité p' 
et de coordonnées x^ et y^ engendre un anneau dont la masse est 

2TZp' [l -{- x') d<T^ . 
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Cherchons le potentiel de cet anneau en un point P du plan 
méridien, dont les coordonnées soient x et y ; il faut d'abord 
calculer les quantités a et b. Si on suppose le point P du même 
côté de l'axe de révolution que l'élément t?<r, on a 

« = v/(x-x')» + (y-Jy')% 

Supposons l'anneau mince et le point P situé près de G, de 
façon que x^ y, x'^ y' soient négligeables devant l; a est une 
quantité du même ordre ; on peut prendre pour valeur approchée 
de i 

car [y — y'Y est négligeal^le ;.il faut maintenant calculer cp (a, b). 
On a 

?(«»*)=t?(t) 

= — ^— — L-^. 

il -\-x-\-x' a * 

Le potentiel en P est donc 

b a 

on a 



il-\-x-\-x' 



comme par hypothèse x et x' sont petits, on peut écrire approxi- 
mativement 



l-^x' 



=î['+^] 



Il vient ensuite 



a a ^ IL 
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ce qui peut s'écrire, avec le degré d'approximation donné 

- b x'-\-x . 

L — r-= = — 1 

IL 21 

on peut alors écrire aussi 

Le potentiel de l'anneau sera donc 

les autres termes étant négligeables, et les intégrales étant 
étendues à la surface de la section. Nous allons supposer l'anneau 
homogène et nous écrirons p' = i. 

Les deux premiers termes sont des potentiels logarithmiques 
qui jouent, dans l'attraction des masses planes, le même rôle que 
le potentiel newtonien dans l'attraction des corps célestes ; les 
propriétés du potentiel logarithmique dans le plan sont anylo- 
gues à celles du potentiel newtonien dans l'espace : on démontre, 
par exemple, que le potentiel logarithmique d'un cercle homo- 
gène est le même que si toute sa masse était concentrée au centre. 

Si le point P est à la surface d'un tore, à une distance a^ du 
centre du cercle générateur, il suffît de considérer le potentiel 
logarithmique en P du cercle de rayon a^ ; sa surface est Tzal^ 
Son potentiel est donc 

■Rol L — ; 
la première intégrale est par suite 

271^0 L 

La seconde est la somme de deux autres : 



— rfd' — / -7- L — 
a fia 



a/-8Z,, i a: ^ SI j 
-y- L —— dd — I -7- L — - dd , 
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dont la seconde se calcule immédiatement, on trouve pour sa 

valeur 

X ,, 81 

la première va se calculer plus tard. 
On trouve de même 





rf(T'=4- I x^d:^=o 



puisque le centre de gravité est à Torigine des coordonnées. 
Il reste à calculer 




nous allons employer un procédé de calcul dont nous nous 
sommes déjà servi dans une autre question : supposons que la 
densité du cercle soit 



2 2 



le potentiel sera 




la densité ne dépendant que de la distance au centre, le potentiel 
logarithmique du cercle est le même que si toute la masse est 
concentrée au centre. Cette masse est facile a calculer, c'est : 
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= =4 "4 -8-J, 



TcaJ 



Le potentiel est donc, en un point A, situé à une distance a^ du 
centre du cercle a^ > o^, 

81 



Trao 



a, 



4 -. 

Supposons que le cercle C se déplace d'une quantité e parallè- 
lement à Taxe des x, de façon à 
occuper la position C^; le poten- 
tiel en A sera (fig. 36) 

dU 
la différence de potentiel sera donc 

àx ' Fig. 36. 

Nous allons étudier d'où provient cette différence : en un 
point de l'intérieur du cercle la densité était p', elle est devenue 

dp' 




p'— e 



dx' • 



Dans la région du cercle C extérieure au cercle C, la densité 
était nulle, elle est devenue égale à p' qui est petit puisque la 
densité p' est nulle sur la circonférence du cercle C, et l'aire de 
la région considérée est de l'ordre de e. Il y a encore à consi- 
dérer l'aire égale à la précédente, intérieure au cercle C et exté- 
rieure au cercle C, la variation de densité en sens contraire de 
la précédente est aussi de l'ordre de e au plus, puisqu'on un 
point de cette aire la densité p', petite sur le bord du cercle, est 
devenue nulle. 

La variation de potentiel due à ces deux aires est donc négli- 
geable, et on a 

ôU r dp' , 8^ 



àx 



àx' 



L 



d^\ 



a 
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Or 

il vîent donc, en supprimant le facteur e, 



,1 8i^, dU 
a ox 



L'intégrale que nous voulons calculer est donc, au facteur Z près, 

, , , ^U ôU - , ,, . . 

égale a -r — ; or — — est la composante de 1 attraction suivant 

.1 X 

Ox ; cette attraction est , sa composante suivant Ox est — ^ . 




a^ a{ 



On a donc 



x' 8f , TZa'^nX 

L dfs = 



l a 4^ ^1 



2 ' 



et, si on suppose le point situé sur la surface du tore, c'est-à-dire 
si Ton a a^ = a,,, il vient pour la valeur de l'intégrale 



2 ^ 



41' 

On a donc pour le potentiel en un point de la surface du tore 
V = 2'ezal L r- Tzal L -^ + -77" '^^If 



8f 


■ + 
al 


5x 

41 



en négligeant les quantités de l'ordre de -^ 

Section méridienne de l'anneau. — Si l'on appelle M la 
masse de Saturne, co la vitesse de rotation de l'anneau, une 
condition nécessaire de l'équilibre est que, à la surface du 

^2132 

tore, la fonction U -j soit une constante, U représentant 

le potentiel de l'anneau sur lui-même augmenté du potentiel de 
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Saturne, ce dernier peut être représenté, avec l'approximation 

M , 

employée , par — j— et Ton doit avoir à la surface de l'anneau : 



1»! 10 11 ^. 



En posant 



M 

l 


«o'R^ _ 

2 


R 


— l + X, 


R' 


— l^ + 2lx ; 



la distance d'un point de la surface au centre du tore est 



ce qui, avec 


l'approximation que Ton cherche, peut s'écrire 








?'-{i+xY, 


d'où 






p—l+x-, 


enfin on a 






I I I :r 




p l + x l l^ ' 


La condition d 


'équilibre peut s'écrire 




2TZa 


?L 


81 X .j SI 5x 




+ 


M 

l 


TT- H \riùHx^=0'\ 

L 2 



il faut donc que le coeflicient de x soit nul, ce qui donne 



4Z Z « l' 

Celte équation donne w'^, aux quantités du deuxième ordre 
près. 

Comme —j- est aussi petit, on a 
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La vitesse de rotation d'un tore engendré par un cercle, de rayon 
suffisamment petit par rapport à sa distance de Saturne, est 
donc la même que celle d'un satellite placé à la même distance. 
S'il n'y a pas de corps central, M = o : on voit alors qu'en 
l'absence de corps central il peut y avoir une figure annulaire 

d'équilibre : a^ est alors très petit, mais le rapport tend vers 

l'infini quand a^ tend vers zéro, car L — augmente alors 

indéfiniment. /'. . 

Si on voulait pousser plus loin l'approximation, il faudrait 
calculer les termes de la forme de 



r 



Q1 

^ a ' 



et 



x""' y"" Ç x"^' y"" d(s' , 



Supposons maintenant que le tore ait une section elliptique ; 
nous devons alors calculer le potentiel logarithmique d'une 
ellipse. Or le potentiel newtonien d'une droite indéfinie en un 
point, est le potentiel logarithmique en ce point de la projection 
du point sur la droite ; il en résulte que le potentiel logarith- 
mique d'une ellipse en un point de son plan est le potentiel 
newtonien d'un cylindre long par rapport à sa section droite, 
en un point situé a une petite distance de sa surface. On est 
donc amené à calculer le potentiel d'un ellipsoïde très allongé 
pour lequel on supposerait c^ infini. 

Le potentiel intérieur est donc une fonction du second degré 

dont les coefficients sont faciles à calculer (page i4o)- 
Pour le potentiel extérieur on a 

\ = k,+ k,x + k,x' + k,y'+ „, 
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On peut calculer A^, A^, k^ comme on Ta fait précédemment, 
mais les intégrales devront être étendues à la surface d^une 
ellipse dont l'équation est 

» •> 

a' b" 

Les quantités du troisième ordre sont négligeables. La fonction 
de force peut de même se développer en série de la forme 

La condition d'équilibre est que l'on ait 

(A. + 0-+(A,+C>^ + (A3 + C3)y^ = C 

a la surface de l'ellipse; il faut donc que l'on ait 

Aj + G, = 

(A,+ C,)a^ = (A3+C,)i^ 

Laplacc avait fait le calcul, et a discuté l'équation mais sans 
tenir compte des termes en A^ et A3. 

M™° de Kowalewski a été plus loin et a fait voir que la section 
méridienne du tore n'était pas symétrique par rapport à Taxe 
des y. 

Nous étudierons plus tard la question de la stabilité ; nous 
allons d'abord discuter l'hypothèse de Cassini. 



to 
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Vraisemblance de l'hypothèse- — Cette hypothèse est la seule 
compatible avec les travaux de Maxwell sur l'anneau de Saturne. 
L'observation confirme cette hypothèse, car l'anneau est trans- 
parent et la lumière le traverse sans trace de réfraction ; il faut 
donc supposer qu'il est composé de corpuscules solides (ou 
liquides), isolées les unes des autres. Les observations spectros- 
copiques montrent d'ailleurs que la vitesse d'une molécule dé 
l'anneau n'est pas la même sur le bord interne, ou sur le bord 
externe. 

Poii^CARÉ. — Figures d'équilibre. i3 
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Equations du mouvement — Nous allons commencer par 
étudier un cas très simple : nous allons considérer p satellites 
de même masse [jl, équidistants sur un cercle ; la distance 

de deux satellites sur ce cercle sera = 2O: le rapport -Ç- 

p *^ 6 

n'est pas infiniment grand. 

Un mouvement possible du système est celui oii chaque satel- 
lite parcourrait le cercle d'un mouvement uniforme, avec une 
vitesse co, déterminée par l'attraction de la planète, augmentée 
d'un terme où il faudrait tenir compte de la masse des autres 
satellites ; et nous allons chercher s'il y a un mouvement peu 
différent de celui-là. 

Soit Mk le satellite de rang k, son rayon vecteur sera i + p,,, 
son angle polaire sera co^ + ^^^ H~ '^k'y dans le mouvement non 
troublé on aurait 

<r/c = o ; 

si le mouvement est stable p^^ doit rester petit. 

La force vive T, et la fonction de forces U, sont fonctions de p,, 
et de (T^] les équations de Lagrange donnent 






ôT dU 

= 0, 



dt X^^'kl à^ji Oi,^ 

rf / ôT \ ôT dU 






dt \ dp',, / dp,, dp,, 



o 



en posant 



T— U= H, 



on a 



d ôH dH 

= o. 



dt dp',, dp 

Nous pouvons supposer que Saturne est fixe, c'est-h-dire que 
nous négligerons les perturbations apportées dans le mouvement 
de Saturne par celui des satellites ; la mise en équation du pro- 
blème est alors indépendante du mouvement de Saturne. 
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La force vive du système est 

et, en négligeant les infiniment petits du troisième ordre : 
La fonction de forces est 

R représentant la partie complémentaire du potentiel due à 
l'attraction des satellites les uns sur les autres ; si l'on pose : 



on a 



^ 2sinA ) 24 8 " ^ 



+ (p/.+P*) ^^S^ cotg.{<+ ^^*7"*^ cotg.{< 



«A — '* ..._j . ('* — «a) 



l«s termes non écrits étant négligeables. On peut écrire aussi 



I + ?fc 



Les termes les plus importants de R sont évidemment ceux 
pour lesquels ^ est petit, sin ^ étant un infiniment petit du pre- 
mier ordre, — . ", est un infiniment ffrand du second, et . ", ' 

sin ^ ^ sm ^ 

est un infiniment grand du troisième. 

U et T comme on le voit ne dépendent pas de t, et sont du second 

degré par rapporta p,p , 5- et <t ; — — , , - sont du premier 
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degré par rapport ii ces quantités; enfin — ^ — ( ■ , jestduprc- 

mier degré avec des dérivées du second ordre. Les équations de 
Lagrange sont donc des équations h coefficients constants de la 
forme 

A étant une constante, et B étant une fonction linéaire des argu- 
ments ; mais ces équations doivent admettre la solution p = o' = o; 
il faut donc que A disparaisse de lui-même ; les équations sont 
alors des équations linéaires et homogènes, qui s'intègrent par 
les exponentielles, et on a, 



p =2 ae>>^ , 



les ). étant des racines d'équations algébriques, et les a étant des 
constantes. 

Le mouvement sera stable si p ne croît pas au delà de toute 
limite ; en mettant en évidence la partie réelle et la partie imagi- 
naire de A, on a 

le mouvement sera stable si \ est négatif ou nul ; mais si nous 
démontrons que Téquation en X a ses racines deux à deux égales 
et de signes contraires, nous verrons que X^ne peut être que nul. 
Or si nous changeons t en — t les équations ne changent pas, 
puisque les termes du second ordre seuls interviennent dans les 
équations de Lagrange ; il en résulte que si l'équation en X 
admet une racine, elle admettra aussi la racine de signe con- 
traire, et par suite, la condition nécessaire de la stabilité de 
l'équilibre, est que les racines de l'équation en X soient purement 
imaginaires. La condition est d'ailleurs évidemment suffisante 
puisque, dans ces conditions, p s'exprime par des fonctions trigo- 
nométriques de la variable réelle t. 

Vous avons 2/? équations du second ordre, donc X satisfait h 
une équation algébrique d'ordre ^p ; c'est une équation en X^, de 
degré ip dont toutes les racines sont réelles et négatives. 
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Les premières des équations de Lagrange sont 

?//—<•>"( I +?/,) — 2(o^V+ I — 2p;i.= ;j. iç^j ; 

on peut écrire, en mettant en évidence dans R les termes de 
degrés o, i et 2 



et Ton a 



OR ôR, OR, 



^Mt, . . OR, 



estune constante, — ^-^ est une fonction linéaire 



Les secondes équations sont 

m, OR.. 

Dans la première équation les termes constants doivent dispa- 
raître : on a donc 

OR. 



tO" I = !JL 



Oo. 



4 Zji sin'i ^ 



'i prend les valeurs 0, 28, 39, ..., (/j — i) 8. Dans le cas ou u.=o 
on aurait (o* = i . 

Les équations sont donc 

^',; — 26)1,; —((O' +2)f;..= \J -T-^, 



<T 



'-»/ 



" .w..^/ 



OR 



c 



atoo,. = fjL 



/. — ^^^1^ — { 



^^/. 



Changement de variables. — Faisons le changement de 
variables, 

Pa — 7. ^v^ * ? i — I,. >.,...,/>>. 
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On a par exemple 

p, = Ç,e2'« H- Ç/« + 5,e«'» + . . . , 

Posons de même 



=2"^^ 



^A= 7 .VI. e'**-^*. 



Les Ç.^ d'une part et les v).^ de l'autre sont deux à deux imagi- 
naires conjuguées ; d'une façon plus précise je dis que Ç., et \p^^ 
sont imaginaires conjuguées. En effet, on peut écrire 

Cette seconde valeur de p^^ se déduit de la première en chan- 
geant i en — / et Ç^ en ^p^.^ : il faut donc que ip_.^ se déduise de E, 
par le changement de / en — i; donc Ç^ et Çp_.^ sont imaginaires 
conjugués. On démontrerait de même que 7)^ et y)p_^ sont imagi- 
naires conjugués. Si p est pair, ^p est réel, car il est à lui-même 

son imaginaire conjugué. 

^ ^ j sera un polynôme du second degré en Ç'v, qui ne ren- 
fermera d'ailleurs que les termes de la forme ^'.^ i^ï— ^ En effet, T 
et R sont symétriques par rapport aux p^ et aux C/^ ; si l'on per- 
mute circulairement les p^ d'une part et les c,^ d'autre part, c'est- 
à-dire si l'on remplace S^ par ^^ e^'ï^, Ç'^ ^^t multiplié par le même 
facteur. 

De même pour v)^ et 7)'^ . Par ce changement ni T ni R ne doivent 
être modifiés : il faut donc qu'un terme S'y SV ^^ transforme en 
lui-même ; donc il faut que le changement de variable multiplie 
ce facteur par e^^^P^ ; il faut donc que 

21 (y + y') = 2ihp ; 
donc 
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En appliquant toujours le même raisonnement, on voit qu'en 
définitive T et R ne doivent contenir que des termes de la forme 

S; Çp - Y Çy '-Op — Y 

Ç Y 'i/» — T ''Y '^P — Y 

Si maintenant on regarde de plus près la formation de la fonc- 
tion II, on voit que Ton peut poser 

II = II,+ H,+ ... + IIy + .. + H^ 
en écrivant IL, à la place de 

i'-. \'p _ ..+2(0 (ç., r/p _ ..+Ç^ _ .;/i'^) + iù% \p _ .^+71 '.^ Yi'p _ ^+ 2i(ip _ . 
+ a I — U; .,;p _ y+^^ (^y'^P - y+*^ y^P - y) + ^y"^ y'^p - y j* 

L, M, N sont donnés par les formules 

V^ / sin^Y^]; cos'i]; cos*^y^ \ 

'^ ^^ \ 4 sin"^^ 2 sin ^ / ' 

\M V^ sin 2y4' cos 4» 

VI / sin-yi|ycos"^ siu^Y^I; \ 

^ ^j \ 2 sin'*^}^ 4 siu ^ / * 



Les équations de Lagrange seront alors 



/ ôII, \ ÔH.^ 



d /ôII. 

Tt 



et on aura en supprimant les indices 

5" — ^ (u,^ + 2) — 2<o-/i'= (X [Li + M/.], 

V'+2u)i' =iA[Mi + Nr.]. 

L'équation en X est 

().i _ (o-^ _ 2 — L^) (À- + Na) — {^^W — 4w*V) =^ o. 
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Discussion des équations. — Quel que soit y, les racines de 
cette équation en A* doivent être réelles et négatives ; on voit 
que cette condition sera remplie si (jl est suffisamment petit ; 

2,3 

Maxwell a démontré qu'il fallait que l'on eût [jl< \ — . Dans 

ces conditions jjl/j tend vers zéro quand /? augmente indéfiniment; 
c'est là une difliculté, car la masse de Tanneau précisément égale 
à |x/;, serait négligeable ; mais cette difficulté disparait d'elle- 
même si on suppose les satellites répartis non sur une circon- 
férence, mais dans un certain volume. 

Si les racines de l'équation en A^ sont négatives, on pourra 
poser 

n étant une quantité réelle, et en prenant 

A = A^e»"% 

on aura 

i=A,e"''-^'>, 

On aura une solution particulière si à ;^ y,^ on ajoute les solu- 
tions, 

7,0 = o. 

et on aura alors une solution particulière pour les p^^ en posant 

à cette solution on peut ajouter la solution imaginaire conjuguée; 
on obtiendra ainsi des solutions de la forme 

PA.= Ap cos (/2; + 'f + aA-yB), 
0-^= Bq cos (w; + ç + aA-yO). 

Si on suppose un système d'axes entraînés constamment avec le 
satellite M;^ dans son mouvement non troublé, et tel que l'un d'eux 
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coïncide avec le rayon OM^, p^^. et Œ/t. seront les coordonnées du 
satellite, par rapport à ces axes, et on voit que le satellite décrit 
par rapport a eux une petite ellipse. 

On passe de la trajectoire d'un satellite à celle du suivant en 
changeant k en A* + i ; le second satellite est donc en retard sur 
le premier, et ainsi de suite ; il y a ce que Ton appelle un déca- 
lage égal il yO. Si AyO est un multiple de 'itz le satellite de rang 
A+i sera dans la même position que le premier satellite; celui 
du rang A + 2 dans la môme position que le second et ainsi de 
suite. 

Le mouvement sera comme s'il se transmettait d'un satellite à 
Tautre avec un certain retard; on peut aussi dire que le mouve- 
ment se transmet comme une onde, dont la longueur serait — '—, 

î 

Les ondes qui sont les plus dangereuses pour la stabilité 

sont les ondes courtes; les ondes qui compromettent le plus la sta- 
bilité sont donc celles qui correspondent aux grandes valeurs 
de y. Considérons en effet une valeur de n, suffisamment grande; 
pour une grande valeur de y, l'onde serait très courte, deux satel- 
lites voisins pourraient se rapprocher d'une manière sensible, et 
leur action mutuelle ne serait plus négligeable par rapport à 
l'action de Saturne. 

La masse totale [x/; doit être inférieure à -^; si on suppose 

p = 100, on voit que la masse de l'anneau devra être inférieure à 

2,3 ,. . 

~ , c'est un inconvénient de la théorie de Maxwell ; mais 

I o 000 

c'est un inconvénient artificiel, car l'hypothèse que nous avons 
faite est trop simple et il faut supposer une distribution de satel- 
lites occupant un certain volume de l'espace. 

La perturbation apportée au mouvement d'un satellite par l'ac- 

tion d'un autre est de l'ordre de —7^^ , en appelant la dislance 

des satellites. 

Si les satellites sont répartis sur une circonférence de 

longueur/, on a 0= : le mouvement sera stable si la per- 
turbation a une valeur moindre qu'un nombre donné A. Dans 
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l'hypothèse de satellites répartis sur une circonférence, il faut 
donc que l'on ait 



ou bien 





<A 


tx< 


..3 



Si les satellites sont répartis sur une couronne comprise entre 
deux cercles de distance comparable à la longueur de l'un deux, 
le nombre des satellites sera comparable à y/. S'ils sont dans un 
anneau de volume fini, leur nombre sera comparable à y/, la 
masse de l'anneau sera égale à [kp^ et la condition de stabilité 
sera que l'on ait 

la masse de l'anneau devra donc être inférieure k un nombre qui 
ne dépend pas du nombre des satellites. 

Stabilité de Vanneau supposé liquide. — Nous allons main- 
tenant, d'après Maxwell, examiner de plus près la stabilité du 
mouvement, dans le cas où l'on suppose l'anneau, soit fluide, 
soit composé de corpuscules indépendants. Dans cette partie du 
travail de Maxwell, les raisonnements ne sont pas tout à fait 
rigoureux, ni même très clairs, mais on doit en accepter la con- 
clusion; il est probable qu'un calcul rigoureux conduirait à des 
résultats peu différents de ceux auxquels nous arriverons. 

Soient x^ y, z les coordonnées d'un point dans le mouvement 
normal; n, v, (v, les composantes de sa vitesse ; 

•^' + i^ y + 'f^^^ y + ç, 

d\ d'f\ dt 

"^HT' ^+-jr' "'+^' 

les mêmes éléments dans le mouvement troublé. 
La force vive du système est 

^di 




(-#)V(..'+^rj^ 



STABILITÉ DE IJ ANNEAU SUPPOSÉ UQUIDE 2o3 

réncrgie potentielle — U est égale à 



(VH -V')l? + ?0 ^^^ 



en appelant V et p le potentiel et la densité dans le mouvement 
normal, V + V, ^-\-^'y ce que deviennent ces quantités dans le 
mouvement troublé. V est une somme de quatre intégrales : 







'A 



La première intégrale est une constante ne dépendant pas 
de i, T^y î^ ; la seconde et la troisième sont égales, comme on l'a 
vu au commencement du cours. Elles représentent la moitié de 
l'énergie potentielle due à l'attraction de la matière répartie, 
comme dans le mouvement normal, sur l'excès de matière qu'il y 
a en chaque point dans le mouvement troublé. 

On peut poser en négligeant le cube de 5, Tj, Ç, 



V + V = V + ;-T \-r, -. h^ 









La somme de la seconde et de la troisième intégrale peut donc 
s'écrire 



U faut maintenant calculer la dernière intégrale; je peux supr 
poser que le déplacement ^, r,, tj est la somme de deux déplace- 
ments, l'un dans le sens de la propagation de ce mouvement, 
l'autre transversalement à ce sens : on posera ainsi : 



y y , f 
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On sait d'ailleurs, d'après les formules connues des vibra- 
tions longitudinales et transversales que 

est une différentielle exacte et que l'on a 

dm d{?r,,) rf(?g ^.^_ 

dx dij dz 

On a d'ailleurs la relation 

^ dx d]i dz 

Dans ces conditions on peut écrire 



V'^^=-gb- / p''^-^[?^+•v^+!;M• 



Prenons maintenant des axes mobiles tournant avec une 
vitesse to autour de l'axe de Tanneau, l'expression de V ne 
change pas. T devient dans ces conditions 




[w+^T- -^^ (.'/+ ■'')]"+[''+ tÎt + '' (•'■+^)] 



+(..■+§)■'.. 



La condition d'équilibre est que 




— lî)rf/ 



soit maximum; si l'on tient compte de la relation 
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on trouve les équations 






V'/ 



O/i ^ d// ' 
(">P . dO 



H- 



P est un polynôme du second degré en (ç, yj, J^), est une 
variable auxiliaire liée aux trois premières par la relation. 

u.r 0// Oi 

Conclusion. — Maxwell est arrivé au résultat suivant : pour 

(tue l'équilibre soit stable, on doit avoir : p<-n — . 
* * ^ ooo 

(]e résultat est vrai pour un amas de poussières cosmiques 
comme pour un anneau liquide ; mais il y a une condition de plus 
pour un anneau liquide : la pression a la surface doit être dirigée 
vers rintérieur ; or, Laplace a trouvé pour l'anneau à section 

elliptique que la densité doit être supérieure a —, pour l'anneau 

extérieur et à 2 pour l'anneau intérieur. 

On peut alors se demander s'il y a d'autres figures d'équilibre. 
Rappelons que l'on a trouvé qu'il fallait que 

Dans cette formule, w doit être pris égal à la vitesse angulaire 
d'un satellite dont l'orbite coïnciderait avec l'anneau. On a donc, 
en désignant par M la masse de Saturne et par L le rayon de 
Tanneau. 

M 



(I)" 



\J 



Si Pq est la densité de Saturne, R son rayon, on a 

M = -|- ttR^o,. 



ao6 AIVyEAU DE SAITHNE 

On trouve donc 



4 W 



p >4 



Po 



m- 



Le second membre est supérieur à — ^ : si ou rapproche ce ré- 
sultat du résultat précédent, on voit que Tanneau de Saturne ne 
peut être liquide. On a d'ailleurs vu plus haut que les observa- 
tions montraient qu'il était probablement composé d'une masse 
de petits corpuscules. 
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